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Le coefficient de dilatation 


Résumé. — Cette seconde étude a pour objet le calcul du coefficient de dilatation et de sa dérivée 
par rapport à la température. 

Le coefficient de dilatation est rapporté directement à l’énergie potentielle et auxtranslations 
élémentaires du milieu cristallin. On trouve uu’il croît avec la température, en gros, comme la 
chaleur spécifique à volume constant. Il s’annule au zéro absolu, et quand la température approche 
du point de fusion, il augmente avec la température absolue à peu près selon la loi parabolique. 


Abstract. — The object of this second study is the calculation of the linear expansion coefficient 
and of its derivative with reference to temperature. 
- The expansion coefficient is directly referred to the potential energy and to the elementary 
vibrations of the crystalline medium. We find that it increases with temperature, roughly as the 
specific heat at constant volume. It is annuled at absolute zero and when the temperature comes 
near to the melting-point, it increases with the absolute temperature nearly according to the para- 
bolic law. 


_ A. Le coefficient de dilatation. — Quand la tem- D'autre part, selon la théorie de l’élasticité, la 
pérature s’élève de dO@, la tension thermique dilatation « dO développe une tension 
T (12) (2) croît de 


; Ca 
Prdr nat = AT: ee AT = (ch + 2 0°?) a«d0 = = d6. (34) 
ar = (TS + + + …) 46. (32) 
et, d’après les expressions de T, (13), de T, (17) Cette tension est égale (en valeur absolue) à 
et de T, (18) l'accroissement (32-33) de la tension thermique. 
’ 46 Compte tenu des fonctions 6 (4) et & (5) et des 
47, _dT,do , 1 Le (5) relations (11), l'égalité des expressions (32) et (34) 
d® dv d® N Fr Lo W \0O / donne : 
æ Le do | 27: de) en 
do du ‘ og du) do fs- 3 ! ; DRE To 
CES JA MEN LUE ST RON PE US7S 
7 DT do, 26 y, AW sic + les Rp fé (urgies Rs) 
d® dv d® pv ” do 
d7,: dT;do | 8e wW2 AWm, (33) + Wa (b + ai |] 
do du d®@ uv ?” dO 
. L TRES 
(2) Les formules dont les numéros vont de 1 à 31, et les eee FF ] W. &. (2b + 2cW | :* Pas 
références bibliographiques dont les numéros vont de [1] mn x 2 lee * m (26 + 3CWm)] ; (35) 
à [5] se trouvent dans la première partie du mémoire « Les ps . De 
facteurs de la dilatation J. Physique Rad. Juillet 1961, 22, d Log.o _ D dw, d Log.g _vdg 
451 ». d Log.v odv’ dLogv gdu' 


LE JOURNAL DE PHYSIQUE ET LE RADIUM.—T. 22. — N° 8-9, AOUT-SEPTEMBRE 1961. 31 


466 


6m est la fonction d’Einstein (4) qui se rapporte à 
à une oscillation harmonique de fréquence 


vm (SE) (16). 


Posons : 
ALI SN FAI De (36) 
PE x 2 ER 
TE Wan 6m (20 + 3cWhn). (37) 
xl £ : US) = d Ke] 
m= ip [£[r Poe AE Log .o] | 
(38) 
K = W2, (b + cWm). (39) 


Comme le facteur Q (14), et pour la même raison, 
à température constante, les facteurs G (36) et 
H (38) sont constants, caractéristiques de l’espèce 
cristalline. Les facteurs J (37) et K (39), et le 
coefficient C (8) ont la même propriété. Et de 
l’égalité (35) nous déduisons le coefficient d’élas- 
ticité : 
: RIGHTS. 
7 ue + 3(H + K) 
Les expressions de G (36) et de J (37) contien- 
nent en facteur la fonction d’Einstein & (4) qui 
s’annule quand la température tend vers le zéro 
absolu. Donc le coefficient de dilatation « s’annule 
aussi quand la température tombe au zéro absolu. 
Dès que les fonctions d’Einstein 6, (4) attei- 
gnent, dans l’ensemble, sensiblement l’unité, ce 
qui se produit déjà à la température ordinaire pour 
la plupart des cristaux envisagés, nous avons gros- 
sièrement : 


(40) 


1e 3kÈ 
. 
D'où il vient : 
x pese (41) 
3k 


En portant dans cette expression les coefficients 
d’élasticité mesurés, et compte tenu des grandeurs 
négligées, J, H et K, on trouve (1) : 

(2) Al  Feæ{() Cu 
€ X 102? gramme 0952072160 2 
(1) Ces facteurs Ë se rapportent aux coefficient € qui sont 


(42) 


déterminés par les modules de compressibilité x (< = 5) 


du fer x à 130 °C. (xre — 0,664 X 10—12) et du cuivre à 
20 0C (xcu — 0,75 X 10—12) mesurés par Grüneisen [6], et 
par le module de compressibilité de l’aluminium à 20 °C, 
mesuré par Bridgmann [7] (xa1 = 1,37 X 16—12), 

(2) Le facteur Ë est homogène à une masse. 

(3) Compte tenu des pulsations «s- déterminées par 
H. Curien [8] on retombe, par l’expression (9), à peu près 
sur la grandeur 1,6 X 10—?? pour le facteur 6 du fer &,si l’on 
prend les coefficients DSP suivants (en unités C. G. $.) : 


DE = 1,6 X 101? pour « — 1 V3 (premiers voisins). 
De$ 1,2 X 101? pour a — 21 (seconds voisins). 
De DR M0 Dour 27 V2 (troisièmes voisins). 


Ces coefficients sont bien compatibles avec les constantes 
de rappel déterminées par H. Curien [8]. 
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Présentement, l’énergie potentielle du milieu 
cristallin n’est pas définie avec une précision suffi- 
sante pour qu’on puisse en déduire, même approxi- 
mativement, ses dérivées cinquièmes et septièmes 
par rapport aux translations qui transportent les 
positions moyennes des atomes les unes sur les 
autres. De ce fait, les coefficients b (28) et c (29), 
et partant les facteurs J (37) et X (39) ne peuvent 
être calculés. D’autre part, les coefficients d’élas- 
ticité mesurés permettent seulement d’atteindre 
l’ordre de grandeur de ces deux facteurs, soit la 
part que les tensions secondaires T, (28) et T3 (29) . 
prennent ensemble à la dilatation, et non la part 
qui revient à chacune d’elles. Le coefficient b oscille 
probablement entre 10-° et 10-10 (unités C. G. S.). 


‘Quant au coefficient c son ordre de grandeur même 


reste inconnu. Afin de développer le calcul de la 
dilatation, je supposerai que ce coefficient s’élève 
aux environs de 10? unités C. G. S. (J’admets, de 
la sorte, que le rapport c/b est grossièrement égal 
au rapport £/b.) 


B. La dérivée 2 du coefflcient de dilatation. — 


dO 
De lexpression (40), on déduit : 
de _ 
d® 
.{4G : dJ = dlog:C. dif dE 
. Es + 2 ee (ue à Log.D a | me 


uC + 3(H + K) 


Le facteur Æ (38) croît avec la température : 


a4 | PH 
do Go), | 
; dH\ dLog.o, ù Log.H 
VE di (o Es d Log.v ‘ en 3; 


Un calcul détaillé révèle que le second terme, 


34 (o (2%) d Log.o | H Ÿ re) 
dw/06 d Log.v 


à Log.v : 
FREE 6 À 
de la dérivée 10 (44) reste, à toute température, 


Le NO 
très petit par rapport au premier 30 ’ et par rapport 


à 
d Log.C dk 
Xæ——— € ù 
d Log.v da® 


k fe L de) e 

3a (ao a} À 
Nous négligerons donc ce second terme dans le 

calcul de la dérivée, c’est-à-dire nous poserons : 


do 
{4G  dJ ; d Log © dH aK 
; Le dE a) — É [uex d Log.v a 5). ) 


DOPESSUHEEUK) (25) 
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Les expressions des facteurs G (36), J (37), dK 
H (88) et K (39) donnent : | de de Po ME EAP a 


LES ENT 
pi o? ol. 


&6\ d Log. & d Log.g 
C4 46 
+ deg [eo ose (a Te Log.v Did rl FF (a) 
aJ . d & 
ao = À En (20 + 66 Wm) + Won © (2? + 3c Wyn). (47) 
D). = 1 6g/: d _. 
Fe = rÈ [re (1 Ta Fl 
gd Es (Sr) | 18 
o d Log.v (90/45 ol 
O,—=10 
(x = 0) 
G 0 
J 0 
11 ] > ['£ (: d Log-g , d ee] 
N £+ le d Log.v  d Log.v 
2 k° &w2 ño 
nie en + 
L (s Fe 2 
dG 
«© À 
J 0 
d® 
dH 
LEA 0 
(ie). 
dK 
= D) 
FAO) 


et, compte tenu de l'expression qui définit la fonc- 
tion F (5), il vient : 


ga d Log. 2) = 
o d Log.v (0 


__ kg d Log.wzx° es (ee + 1) 
© » wd Log.v Log.v (et — 1) 


(50) 


C. Les grandeurs qui déterminent le coefficient 

de dilatation. — Les facteurs G, J, H, K et leurs 
G dJ SE ) dK 

40’ 49’ (36 /,,'d® 

et à température élevée, les grandeurs suivantes : 


dérivées —— ——ont, au zéro absolu, 


OSLE "# 


(z : 
3È 
kO (2b + 3ckO) 
d Log.g 


1 gkO d Log. 
roles ie Log .v 24 a le 


œ 3Ék® (d’après la relation (26)) 


© k262(b + ck@) 


e DEA ORUNONE ( Log.g 9 d rxe)] 
N £ L6k 05 d Log.v d Log.v/Jsr 
° : 
en. . SE £gsr (d’après la relation (26)). 
ST 


2k(b + 3ckO) 

k LS (1 “dl DES +2 Ee)] 
N F7 Low? d Log. d Log.v/]s7 
= 36k (d’après la relation (26)). 

k?O0(2b + 3ckO) = kJ. 


Quand la température s’éiève dès le zéro absolu, 
le facteur A (38) diminue d’abord légèrement, 
puis augmente et finit par croître comme 3Ë40. 
D'autre part, le facteur X (39) augmente cons- 
tamment à température croissante. Mais, même 
-orsque la température est à peine inférieure au 
point de fusion, la somme 3 (4 + K) n’atteint que 
le dixième, environ, du produit €. Et, comme le 
coefficient © (8) diminue rapidement lorsque la 
dilatation grandit — en gros comme une fonction 
exponentielle décroissante de la dilatation —, dans 
les formules (40) et (45), qui AN le coefficient 
de dilatation «& et sa dérivée — Te le dénominateur 
uC + 3(1 + K) diminue sans cesse à température 
croissante. 


la dérivée TE (46) croît 


d’abord lorsque la température s’élève ; elle passe 
par un maximum à une température y plus basse 


Nulle au zéro absolu, 


que la température où la dérivée LE Cal) (for- 
mule 6, fig. 1) s’annule, c’est-à-dire : 
hvm (Ou) 
On < ——— RTE. (51) 


Yu(Ow) étant la fréquence moyenne (16) des oscilla- 
tions fondamentales à la température Oy (1) ; puis la 


(*) Quand la température s’abaisse, la fréquence moyenne 
Vmn (27 ) s’élève donc 


M (Om) > vm EI 


A ue 
212 2R 
RME 2: 


en || 
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k uand la température s’élève, la fonction 
Trés décroît, et, quand —""@» {ombe au- Q P è 
dO kO 3 1 P>2W? HAE ee w = ) 
dessous de l’unité, elle diminue sensiblement ce (sor ). Ed A ps 5 he 
comme 1/3. k : 
nulle au zéro absolu, commence par croître, passe 
dG ,402 par un maximum de 1,45 pour x égal à 3,9, puis 
Le L diminue, et, lorsque x devient inférieur à l'unité, 


elle tend vers l’unité (courbe de la figure 4). 


F1G:03: 


La courbe de la figure 3 représente, en fonction 


de ©, la variation de la dérivée TS pour le fer «. La fonction 

Elle se rapporte à la densité p(v), déterminée par : ( 5) Sr (2x +W ($) } 

H. Curien [8], son maximum se situe à 58 0K envi- ROSE RE 29/w 

ron (?). 4 s’annule aussi quand la température tend vers le 

D'autre part, les expressions des tensions T, (28) zéro absolu ; et lorsque la température s’élève elle 

et 73 (29) donnent : augmente constamment, sensiblement comme © 
aJ Ar, d?W2 da? W£, dès que x décroît au-dessous de l’unité. Toutefois, 
cie Al MOST Ï quand elle devient égale à 0,82 environ, pour 


x — 3,7, sa dérivée seconde s’annule (courbe de la 
L'énergie moyenne W, (3) est celle d’un oscil- figure 5). 
ÉCHAE Or En conséquence, lorsque la température s’élève à 
lateur linéaire dont la fréquence vn (5°) (27) est 


À Se : 
L partir du zéro absolu, la dérivée —, nulle à cette 


. , L 
constante. Les dérivées ù se 
température, commence par croitre, passe par un 
d'w? n d'Wy L LATE 
ACTE 0? ? maximum pOur 76 — légèrement plus grand 
varient donc en fonction de la température comme que 3,9, donc à une température plus élevée que la 
les dérivées partielles température Om (51), puis décroît légèrement, 
0 V 
D? W? 22 W3 tombe à un minimum pour -—— aux environs de 1, 
CONS EN 10? Ja © 
[0] 


enfin croît en gros comme ©. Mais, compte tenu des 
données expérimentales, le maximum de la déri- 


(HsPour le lens à va (Ou) 20240 
ee __ = 66 0K. vée 5 reste petit par rapport au maximum de la 
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dérivée _. (probablement inférieur au dixième de 


ce dernier). Donc, aux basses températures, la va- 


riation de la dérivée Se en fonction de ©, est déter- 


en >243 
ETATS EN 
3 xlircost?x) -sinh x 
x 


| td G 
minée principalement par celle de la dérivée 


Et, tant que la température @ reste inférieure 
à Om, nous avons approximativement : 


dG 
dax do 
OA UC 


Aux température élevées, quand la différence 


relative LS tombe au-dessous du dixième 
F 
Or 
lin (5) LARÉCTANINRERESS 
0 > , la dérivée —— devient nég'igea- 


k dO 


ble parrapport à la dérivée 10: Donc, aux tempéra- 


tures élevées, nous avons de même approxima- 
tivement : 


du k2 
© _uC + 3{4 + À) 


(52) 


[2 (b + 3 ckO) + o (O)], 
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__ 3 [3E + kO (26 + 30]? 
7 uC+3(4 +4) 


1 Trees! 
x b JEU KI [à Logo [ 


ue 


o(@) 


Notons que pour tous les cristaux cubiques envi- 
ar d 00 
sagés le module de la dérivée TE excède forte- 
: d Log v 
ment l’unité (pour les cristaux de cuivre il s’élève 
aux environs de 5. 


D. La variation du coefficient de dilatation. — 
Quand la température tend vers le zéro absolu, le 
coefficient de dilatation x (40) et les dérivées 
dœ 
dO 
tombe aussi à zéro (courbe OD, fig. 6). Lorsque la 
température s’élève à partir du zéro absolu, la déri- 


RÉ TR A AA ; s PRET 
10 (46) et 10 4) s’annulent ; donc la dérivée 


demo" 


Cuivre 


FIG. 6. 


La courbe OD concerne la dilatation des cristaux de cuivre. 
Elle se rapporte aux grandeurs suivantes : 

1) densité des fréquences p(v) à 290 0K déterminée 
par E. Jacobsen [9], et compte tenu de mesures plus 
récentes, par la diffusion des neutrons, faites par D. Cri- 
bier (non publiées) ; 

Y étant la fréquence au zéro absolu, on admet que : 


(0) 
v(@) — v, e—aë@) ; he PARA 2 oO = (O) «© ; 
(@) o P 


E— 92 x 10—2? (25): b—3 x 10—1 (28); 
Ci AO = 025 1 0722 nnités CGR 
Co étant la grandeur du coefficient € (8) au zéro absolu, 

on a pris : 
C(O) = €, eve) ; 


) Ge erilimeiC CG CAE 
PA 
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do 


dO 
qui, compte tenu du facteur J (37-47), se situe à 
une température Oy légèrement plus haute que la 
température Oy (51), mais inférieure, de peu, à la 
température  . (Yu(Ou) étant la fré- 
quence moyenne des oscillations fondamentales 
dans le cristal à la température Ow) ; puis la déri- 
de 


10 
constante aux températures moyennes (branche 
BC) ; enfin, quand la température atteint le voi- 
sinage du point de fusion, elle augmente de 
nouveau avec © sensiblement selon l’expres- 
sion (52), donc plus vite qu’elle croîtrait selon la 
loi linéaire (branche CD). 

En conséquence, quand la température s’élève, 
dès le zéro absolu, le coefficient de dilatation «, 
nul à l’origine (courbe OD, fig. 7), commence par 
croître lentement (branche OA), puis la courbe qui 
décrit sa variation en fonction de @ présente un 


vée —= croît (branche OA), passe par un maximum 


décroît (branche AB), devient à peu près 


x x10° 


Cuivre 


F1G. 7. 

Coefficient de dilatation linéaire du cuivre calculé à partir 
des grandeurs prises en compte pour évaluer la dérivée 
dœ RE , : : : 

Jo: Grüneisen a mesuré la dilatation du cuivre entre 


186 °K et 873 °K, Dans cet intervalle de température, 
les coefficients de dilatation calculés se confondent à 
16—? près avecles coefficients mesurés par Grüneisen [10]. 
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point d’inflexion Z à la température On 
(branche AB); ensuite, aux températures 
moyennes, le coefficient croît approximativement 
selon la loi linéaire (branche BC) ; enfin, aux tem- : 
pératures élevées, proches du point de fusion, il 
augmente sensiblement comme 


k [3€ + A0(2b + 3ck®)] 
uC+3H+X) 
donc plus vite que selon la loi parabolique (53) . 
(branche CD) ; l’écart entre sa loi de croissance et . 
la loi parabolique étant dû principalement à la di- 
minution du dénominateur uC + 3(4 + K). | 
J’admets qu’à toute température chaque oscil- 
lation fondamentale conserve une énergie moyenne 
sensiblement égale à celle d’un oscillateur harmo- 
nique ayant la même fréquence (3). Cela n’est pas 
démontré. Si l’énergie moyenne véritable, W;, d’une 
oscillation fondamentale ne se confond pas avec 
l’énergie moyenne W (3) d’un oscillateur harmo- 
nique, l’écart relatif [W; — W]/W ne peut être 
notable qu'aux températures élevées. Toutefois, le 
coefficient de dilatation est principalement fonc- 
tion, non des énergies W,4., mais des dérivécs ” 
ee et un faible écart |[W; — W}\5./Ws. n’im- 
plique pas que la différence relative 
APP AE 
d®O  dO s-/ |«® 


soit minime. Si, aux températures élevées, les diffé- 
rences relatives 


, 


EE 


— 1] x 


deviennent considérables, le coefficient de dila- 
tation est, à ces températures, plus grand ou plus 
petit que le coefficient prévu, selon que les écarts 


FA 
TO 


k, 


sont positifs ou négatifs. 

Notons que, même si l’énergie moyenne W; d’une 
oscillation fondamentale ne se confond pas avec 
celle d’un oscillateur harmonique, le coefficient 
d’élasticité « reste déterminé, comme je l’ai établi, 
par les fonctions (4) et (5) et les dérivées (6) et (7). 
Mais ces fonctions et ces dérivées se rapportent à 
l’énergie moyenne véritable W,. 

La dilatation thermique des cristaux métalliques 
a été l’objet de nombreuses mesures. Mais presque 
toutes ont été faites sur des cristaux à des tempé- 
ratures moyennes (branche BC de la courbe OD, 
fig. 7) et les intervalles de températures qu’elles 
couvrent sont le plus souvent fort restreints (quel- 
ques dizaines de degrés centésimaux). Il importe de 
connaître la dilatation dans deux zones de tempé- 
ratures restées Jusqu’ici inexplorées. 
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19 Aux basses températures, afin de vérifier, 
d’une part, si le coeflicient de dilatation et sa déri- 
vée s’annulent au zéro absolu, d'autre part si la 
courbe (0D, fig. 7) qui décrit la variation de ce 
coefficient en fonction de la température porte un 
point d’inflexion. Si ce point d’inflexion existe, il 
est possible de repérer la température O} où il se 
situe, et par cette température de déterminer appro- 
ximativement la fréquence moyenne vy(@y) des 
oscillations fondamentales qui ont lieu dans le 
cristal à la température @y (1). 


= 4,5kO; 
YM(O@m — h Mote 


(*) La fréquence moyenne vu (@m) trouvée se rapporte 
à la température Om, et non à la température an et 


comme les fréquences des oscillations fondamentales 
s’abaissent quand la température s’élève, on a 


re (Oo > in (%) 
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20 Mais c’est surtout aux températures élevées, 
immédiatement au-dessous du point de fusion, 
qu’il‘convient de mesurer la dilatation avec une 
haute précision. Cette mesure révèlerait l’impor- 
tance des tensions thermiques secondaires, et sinon 
la grandeur exacte, tout au moins l’ordre de gran- 
deur des dérivées cinquièmes et septièmes de l’éner- 
gie potentielle, et peut-être même la variations des 
facteurs Esesr (30) et Ésssrs (31) en fonction 
de la température. 

Remarquons enfin que le coefficient de dilatation 
est essentiellement fonction des rapports de gran- 
deurs entre les dérivées successives de l’énergie po- 
tentielle et non de la grandeur absolue de ces déri- 
vées. 


Manuscrit reçu le 11 avril 1961. 
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L'EFFET PHOTOÉLECTRIQUE FONDAMENTAL DES SEMI-CONDUCTEURS ET DU Cs;,Sb 


Par A. MEESSEN, 


Centre de Physique Nucléaire, Université de Louvain, Belgique. 


Résumé, — Une théorie de l’effet photoélectrique fondamental est proposée, pour le cas de 
transitions directes, verticales et pour Je cas de transitions indirectes entre la bande de valence 
et la bande de conduction. La confrontation et la discussion des résultats expérimentaux pour 
le Cs,Sb permet de détacher l’effet photoélectrique fondamental d’effets secondaires et de donner 
quelques orientations nouvelles pour d’autres recherches expérimentales. La profondeur d'émission 
des photoélectrons diminue quand on s’approche du seuil photoélectrique, ce qui est dû aux diffu- 
sions sur les centres d’impuretés. 


Abstract. — A theory of the fundamental photoelectric effect is proposed for both direct 
and indirect transitions from the valence band to the conduction band. Comparison and discus- 
sion of the numerous experimental results for Cs,Sb allows one to separate the fundamental photo- 
electric effect from secondary effects and to suggest new lines for other experimental research. 
The photoelectric emission depth becomes smaller near to the photoelectric threshold of the elec- 


tron scattering by impurity centers. 


introduction. — L'effet photoélectrique des mé- 
taux a été décrit théoriquement comme un effet de 
surface. Ce n’est que récemment qu’on a été con- 
duit à l'introduction d’une nouvelle théorie basée 
sur l’idée d’un effet de volume [1, 2]. 

L'effet photoélectrique des semiconducteurs par 
contre fut toujours considéré comme un effet de 
volume, surtout parce que la distribution spectrale 
de l’efficience photo-électrique reflète celle de l’ab- 
sorption optique qui est un effet de volume, mais 


B. de cond 


électrique correspondant à une absorption aux 


centres d’impuretés ou de déformation du réseau ou : 


encore à une absorption avec formation d’excitons 
Mais celui-ci sera une perturbation de l'effet photo- 
électrique fondamental. 

Nous ne considèrerons donc que la bande de 
valence et la bande de conduction, séparés par un 
intervalle d'énergie AE et une affinité électro- 
nique W pour les électrons de la bande de conduc- 
tion (fig. 1A). Cette théorie s’appliquera donc aussi 


B. de val. 


F1G. 1, — A: Modèle. — B : Transilions directes. — C. Transitions indirectes. 


aussi parce que l'efficience photoélectrique est bien 
plus grande que dans le cas des métaux. D'ailleurs 
Burton [3] a estimé la profondeur d'émission des 
photoélectrons dans le cas du semiconducteur 
Cs,Sb, à environ 250 À. 

Les semiconducteurs présentent en général un 
effet photoélectrique dont le caractère est assez 
complexe et variable d’un cas à l’autre. Nous ne 
considèrerons 1c1 que l’effet photoélectrique fonda- 
mental, correspondant à l’absorption fondamentale. 
On peut évidemment observer aussi un effet photo- 


bien aux isolants qu'aux semiconducteurs, mais le 
seuil photoélectrique des isolants sera reculé plus 
loin dans l’ultraviolet. 


Théorie de l'effet photoélectrique fondamental. 
— &) LA FORMULE DE APKER, TAFT ET DICKEY. — 
Apker, Taft et Dickey [4] sont les seuls auteurs qui 
ont proposé une formule générale pour la distri- 
bution spectrale de leffet photoélectrique fonda- 
mental des semiconducteurs. Ils utilisent le modèle 
sulvant, 
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— La bande de valence correspond à un dia- 
gramme e(%) parabolique, tel que la distribution 
d'énergie de ces électrons soit donnée par 

n(e) f(e, T) + ef 
pour ÀT < AE. 

— On fait abstraction de toute règle de sélection 
en X. Des photons d’énergie Av pourront donc 
exciter des électrons de la bande de valence de 
n'importe quelle énergie e. Le seuil photoélectrique 
sera Av, — AE + W et l’énergie des photoélectrons 
émis sera comprise entre £ — O0 et 


La probabilité de photoexcitation d’un électron de 
la bande de valence d’énergie € par un photon 
d'énergie Av n’est pas explicitée. On suppose sim- 
-plement qu’elle peut s’écrire 


se, v) = g{v) en—5/?, 


où g(v) varie lentement. 

— On suppose alors que la probabilité de sortie 
des électrons excités soit p — En. Ceci peut être 
rendu plausible par le raisonnement suivant : 
l’énergie cinétique de l’électron excité dans la bande 
de conduction est e — E + W. Mais cet électron ne 
pourra sortir que si la composante normale de son 
énergie cinétique est supérieure à la hauteur de la 
barrière d’énergie à la surface W. En supposant 
donc que la réflexion de la surface soit nulle dans 
ce cas et qu’il n’y ait pas d'absorption des électrons 
sur leur parcours vers la surface, on aura pour la 
probabilité de sortie 


p={112}(41—VWIJe) ElsW 


près du seuil photoélectrique où 0 < ÆE < Es. 
Avec ces hypothèses, le courant photoélectrique 
sera finalement 


iv) — se n(e) f(e, T)s(e, v) pde 
 q(v) (hv — hv,jn 


Ces auteurs reportent pour le Te et le Ge une 
expression de ce genre avec m = 4. 

Ce modèle aboutit donc à une formule simple et 
apparamment assez générale. Mais en fait il ne 
pourra être valable que dans le cas des transitions 
indirectes, comme nous le montrerons plus loin. On 
sait par contre que l’absorption fondamentale des 
semiconducteurs et isolants est déterminée surtout 
par les transitions directes, verticales avec une règle 


de sélection très stricte : 4° = k. 


b) L’EFFET PHOTOÉLECTRIQUE FONDAMENTAL 
DANS LE CAS DES TRANSITIONS DIRECTES. — L'étude 
de l’absorption optique des solides [5] aboutit aux 
conclusions suivantes : parmi les transitions di- 
rectes, c’est-à-dire sans intervention des vibrations 
du réseau, on peut distinguer les transitions verti- 
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cales (% — %) et les transitions obliques 


(k! — k + G). Ces dernières résultent d’une réflec- 
tion de Bragg des électrons sur le réseau, simul- 
tanée avec la photoexcitation. Les transitions verti- 
cales, donnent lieu à un coefficient d'absorption 
plus important. Mais dans les métaux (alcalins) 
seul les transitions obliques interviennent parce 
que les transitions verticales sont interdites. Cela 
n’est pas le cas pour les semiconducteurs et isolants 
où les transitions verticales sont prédominantes. 

Dans les transitions directes, verticales l’énergie 
des électrons excités est bien déterminée. On a donc 
toute raison d’appliquer l’approximation mono- 
cinétique qui permet de calculer le courant photo- 
électrique pour un film d’épaisseur d, caractérisé 
par une profondeur d'émission des photoélec- 
trons ZL. Comme pour le cas des métaux [2], on 
trouve : : 

iv) = ea gl, pLii —e—dL) = Map 


I, est l'intensité lumineuse incidente et g7, est 
l'intensité lumineuse dans le film, déterminée par 
les propriétés optiques classiques de la surface du 
film. g et L peuvent dépendre de la fréquence, 
donc M aussi. Mais il est probable que M varie len- 
‘tement près du seuil photoélectrique et que ap 
donne le comportement caractéristique. 

Le coefficient d'absorption « est donné par [5, 6] 


a(v) = (AJhv) (kv — AE)? 
pour des transitions permises (par ex. : s — p) et 
a(v) = (BJhy) (hv — AE)? 


pour des transitions non permises (par ex. : s —s). 

Ce résultat peut être retrouvé très simplement 
puisque la probabilité de photoexcitation pour un 
électron d'énergie e est s(e, v) — Cte ou « et que 
a(v) + n(e)s(e, v) avec n (e)  el”?, Si de plus les 
sommets des paraboles e(k) correspondent à la 
même valeur de k (par ex. : fig. 1B)on a 


E (Av — AE). 


Mais dans ce cas on a aussi £ + (hv — hv,) et 
comme p + Æ près du seuil photoélectrique, on 
aura dans cette région 


uv) æ (hv — hy,) (hv — AE)V? où 3/2, 


Suivant que les deux paraboles e(k) sont tournées 
en sens opposés (fig. 1B) ou dans le même sens, 
on aura respectivement pour le seuil photoélec- 
trique : Av, — AE + W(1 + y) ou encore 
hv, = AE + 2W et hv, = AE puisque y est le 
rapport des masses effectives. 

Donc dans le dernier cas seulement on peut avoir 


i(v) © (hy — AE)? où 5/2, 
Notons que dans le cas de films minces (d & L) 


on aurait encore pour l'efficience quantique par 
photon absorbé, comme dans le cas des métaux 
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Y ap. L’eflicience photoélectrique des semi- 
conducteurs du type Cs;Sb n’est pas seulement 
plus grande que celle des métaux à cause de la 
plus grande valeur de «, mais aussi à cause de l’aug- 
mentation plus rapide de p, suite à la valeur plus 
petite de l’affinité électronique W. Pour les métaux 
W = 5 à 15 eV et pour les semiconducteurs 
We tev (fs. 3). 

L'analyse de l’effet photoélectrique fondamen- 
tal dans le cas des transitions directes, verticales, 
peut se faire de la façon suivante : 

1. Comparer le spectre d'absorption déterminé 
expérimentalement, avec une des formules théo- 
riques de œ{v). 

2. Déterminer AE et Av, par le spectre d’absorp- 
tion et par le spectre de l’efficience photoélectrique. 
Ceci permettra de calculer 

1 AN, — AE) 
P = À (: = VE — a) 
si les sommets des paraboles e(k) correspondent à 
la même valeur de #. 

3. Diviser 2(v) par le produit «p pour obtenir une 

idée de la variation de W(v) et de L{(v). 


€) L'EFFET PHOTOÉLECTRIQUE FONDAMENTAL 
DANS LE CAS DES TRANSITIONS INDIRECTES. — Les 
transitions indirectes impliquent une émission ou 
absorption simultanée d’un phonon d’énergie Æ, 
avec l’absorption du photon d’énergie Av. Les tran- 
sitions indirectes constituent donc un effet du 
second ordre et sont moins probables que les tran- 
sitions directes. 

Les transitions indirectes ne sont à considérer 
que dans le cas où les transitions directes sont 
interdites. Ce sera le cas par exemple quand les 
paraboles e(k) n’ont pas leur sommets pour une 
même valeur de 4 (fig. 1C) et que le seuil d’absorp- 
tion par transitions directes verticales est supérieur 
à AF. Mais pour retomber sur le modèle de Apker, 
Taft et Dickey, il faut encore que la bande de 
conduction (ou la bande de valence) présente un 
diagramme e(k) assez aplati, c’est-à-dire ne cou- 
vrant qu’un intervalle d’énergie de l’ordre de 2 XT. 
C’est le cas du Ge auquel la théorie des transitions 
indirectes peut très bien s'appliquer. En effet, 
l’énergie de l’électron excité sera 

e —=hy— AE —sT+tE, 
respectivement suivant qu'il y a eu absorption ou 
émission d’un phonon. 

Comme pratiquement Æ£, < kT, on couvrira une 
bande d’énergie de l’ordre de 2 AT. 

Admettons donc que hvy — AE = kT. Alors 
l'énergie des électrons qui peuvent être excités par 
des photons Av sera0 Le <ey =hv—AÂE+HE, 
et l’énergie des électrons qui peuvent être excités et 
émis sera 0 < € < Em = Av — hv, avec 


AVS SAE L'W. 


comme seuil photoélectrique. 
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Le coefficient d'absorption a été calculé par 
Bardeen [6] et peut être calculé plus simplement. 
d’après [7] 
ab) = [M ne) ste, v) ne de + (av — AE à Epjreu 


Si nous utilisons encore l’approximation mono- 
cinétique pour le calcul du courant photoélectrique, 


nous devons faire intervenir par contre un autre. 


coefficient d'absorption, correspondant aux élec- 
trons qui peuvent être émis : 5 


iv) = Ma'p avec 


RASE Je n(e) s(e, v) n(e’) de  (hv — hv,)? oui, 


Ici nous utilisons pour la probabilité de photo- 
excitation s(e, v) — Cte ou €. Pour la probabilité 
de sortie nous admettons encore p + E,. Ainsi nous 
obtenons près du seuil photoélectrique 


iv)  (hv — hv,)ou4, 


C’est le résultat trouvé par Apker, Taft et Dickey 
pour le Te et le Ge avec l’exposant 4, ce qui indi- 


querait que la transition à eu lieu entre états des 


même parité. 


Données expérimentales sur le Cs;,Sb. — Les … 


photocathodes au Cs-Sb, dont l'efficience pour la 


lumière visible est de 20 à 30 0), furent découvertes 
par Gôrlich [8]. Il est bien démontré actuellement 
que ce composé intermétallique forme un semi- 
conducteur pouvant être du type p, puisqu'il se 
forme un composé Cs,;Sb avec éventuellement un 
léger excès de Cs [9, 10]. 

Nous avons choisi le Cs;Sb parce que son effet 
photoélectrique, étudié par beaucoup d’auteurs à 
cause de sa grandeur et de son importance pratique, 
est probablement un des meilleurs cas d’effet photo- 
électrique fondamental connu actuellement. En 
rassemblant les résultats expérimentaux on 
observe de grandes divergences qui montrent com- 
bien il est dangereux de se baser uniquement sur 
des résultats partiels d’un seul auteur, à cause de la 
densité variable des centres d’impuretés ou des 
imperfections dans les différents échantillons. On 
en tirera aussi quelques orientations nouvelles pour 
d’autres recherches expérimentales. 


&) LE SPECTRE D’ABSORPTION OPTIQUE. — Le 
spectre d'absorption optique du Cs,Sb a été mesuré 
par Burton [3]; Wallis [9], Spicer, Morgulis ainsi 
que par Taft et Philipp [11] dont les résultats sont … 
comparés à la figure 2. La courbe en trait plein 
correspond à la formule théorique pour des transi- 
tions directes, verticales, non permises : 


av) = (1,1 405 Jhv) (hv — 1,6)#2 cm—1, 


Il en résulte que AE — 1,6 eV. L'écart pour 
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._ kv > 3,5 eV est dû à l’approximation de la formule 


qui n’est valable que près du seuil photoélectrique. 
Les autres déviations peuvent être considérées 
comme des perturbations de l’absorption fonda- 
mentale. 


ÆD=0-0— 0 
ui 0-0 O0-0-9. 
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1. Pour hv < 1,6 eV on observe encore une 
faible absorption, que Wallis [11] à interprété 
comme un effet des centres d’impuretés dont la den- 
sité serait d'environ 101? cm—#, 

2. Autour de hv — 2,4 eV, on observe dans cer- 
tains cas une autre absorption supplémentaire quise 
compose probablement de plusieurs maxima deve- 
nant plus visibles à basse température. Il s’agit soit 
d’une absorption avec formation d’excitons, soit 
d’une absorption aux centres d’impuretés, soit 
encore d’une superposition des deux. 


b) LA DISTRIBUTION SPECTRALE DE L'EFFICIENCE 
PHOTOÉLECTRIQUE. — Les résultats expérimentaux 
de Burton [11] ; Zworykin et Ramberg [9], Sakata ; 
Miyazawa ; Apker, Taft et Dickey ; Eckart [12] 
ainsi que Spicer ; Taft et Philipp [11] sont com- 
parés à la figure 3. Le parallélisme entre ces courbes 
et celles de l’absorption optique est évident. Il 
s’agit bien d’un effet photoélectrique fondamental 
auquel se superposent encore des effets secondaires. 

4. Pour hky = 1,2 eV Tañft et Philipp [11] trou- 
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vent un maximum à 90 OK, si la surface a été 
irradiée préalablement par une lumière de fré- 
quence plus élevée. Mais ce maximum s’efface pro- 
gressivement sous l’action d’une irradiation à 
1,2 eV. 
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Schaetti et Baumgartner [13] ont fait une obser- 
vation semblable. Ils ont montré aussi que le cou- 
rant noir augmente progressivement avec la durée 
d'exposition préliminaire, mais qu'une saturation 
est atteinte après 5s pour le Cs;Sb et après 1205 pour 
le Li,Sb. Le courant noir dépend également de la 
fréquence utilisée dans lillumination préliminaire 
et présente une dépendance spectrale semblable à 
la distribution spectrale de l'effet photoélectrique, 
avec une influence plus prononcée de la lumière 
rouge (1,6 eV < hv < 2,5 eV). 

Nous en proposons l'interprétation suivante : 
l'excitation préliminaire a pour effet d’exciter des 
électrons de la bande de valence à la bande de 
conduction. Mais un certain nombre de ces élec- 
trons sera retenu dans la bande de conduction ou 
sera capté par les centres d’impuretés ou retourne à 
la bande de valence — cela à cause des diffusions 
sur le parcours vers la surface et à cause de la 
barrière de potentiel W à la surface. L2 courant 
noir sera produit par une émission thermique dos 
électrons qui ne sont pas encore retournés à la 
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bande de valence. Le seuil de cet effet doit être 
évidemment AE — 1,6 eV. L'émission photoélec- 
trique pour Av — 1,2 eV sera due à l’excitation 
d'électrons captés dans les centres d’impuretés. 

2. Vers Av — 2,4 eV on observe un autre maxi- 
mum secondaire, qui se détache d'avantage à 90 °K 
et se localise alors à 2,2 eV d’après les mesures de 
Shalabutov et Maslenikova [15]; Taft et 
Philipp [11] ainsi que Miyazawa [12]. La courbe en 
trait plein semble bien représenter l’effet photo- 
électrique fondamental. 

La figure 4 donne la comparaison des courbes 
d’efficience photoélectrique pour différents com- 
posés du type M,Sb, mesurées par Spicer [11]. On 
remarque l’analogie entre les courbes du Cs,Sb et 
Rb,Sb d’une part et des courbes du K,Sb et Na,Sb 
d'autre part. 

Sommer [10] a montré que les deux derniers sont 
normalement du type nr et que les deux premiers 
sont normalement du type p. On peut donc 
admettre une émission d’impuretés p autour de 
2,4 eV pour les deux premiers, absente pour les 
deux derniers. Ce point de vue est confirmé par 
l'apparition de nouveaux maxima dans la même 
région du spectre d’absorption après addition de K 
au K,Sb et de Na au Na,Sb si on admet qu’on pro- 
duit ainsi des centres d’impuretés p. Il semble donc 
bien s’agir d’un effet d’impuretés. Mais l’existence 
de plusieurs maxima dans le spectre d'absorption 
du K,Sb, observés par Spicer ainsi que par Taft et 
Philipp [11] pourrait indiquer qu’il existe aussi un 
effet d’excitons, plus visible dans le K,Sb et le 
Na,Sb qui ne sont pas du type p. Il serait intéres- 
sant de vérifier si l’émission photoélectrique autour 
de hv — 2,4 eV dépend également de l’excitation 
préliminaire, puisqu'elle devient plus prononcée à 
basse température et puisqu'elle n’est pas égale- 
ment importante chez différents auteurs. 
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Remarquons la grande efficience photoélectrique 
des photocathodes multialcalines (fig. 4) et la dimi- 
nution du travail d'extraction photoélectrique par 
Pabsorption de Rb et de Cs sur la surface. Ce sont 
les meilleures photocathodes connues actuellement. 

L'influence de la température sur l'efficience 
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photoélectrique du Cs,Sb a déjà fait l’objet d’une. 
autre étude [14]. Notons ici simplement que le | 
déplacement de 3 à 5 10—4eV/degré du seuil photo- | 
électrique Av, avec la température s’explique direc- 
tement par le fait que 


AE = AE, — cT avec c = 10—* eV/degré. 


Il n’est donc pas nécessaire de faire appel à un 
effet d’impuretés, comme Spicer [11] la supposé. 


C) LA DISTRIBUTION D'ÉNERGIE DES PHOTOÉLEC- 
TRONS. — Les déterminations expérimentales de la . 
distribution d'énergie des photoélectrons du Cs,Sb 
pour différentes valeurs de Av sont rassemblées à la 


r 


Fe 
hv(ev) 


res, 


figure 5. Ces courbes sont normalisées par leur 
maximum et déplacées suivant la valeur de Av. Ces 
mesures ont été faites par Sakata ; Miyazawa ; 
Apker, Taft, Dickey [12]; Shalabutov, Masleni- 
kova ; Pyatnitskii ; Zheludeva, Akhmatova [15]; 
Taft, Philipp [11]. 

Le trait essentiel de ces courbes est l’apparition 
d’un second maximum à basse énergie pour 
hy > & eV. Comme dans le cas analogue des mé- 
taux [1, 2] on peut considérer ces électrons lents 
comme des électrons secondaires. Il faut alors que 
hy D Dys —= 2(A Eh + W). 
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Notons que Apker, Taft et Dickey ainsi que Taft 
et Philipp [12] observent plusieurs maxima dans 
les courbes de distribution d’énergie pour 
hky > 4 eV. Shalabutov et Maslenikova [15] par 
contre observent un maximum secondaire pour 
hv — 2,34 eV. 


d) LES EFFETS THERMIQUES. — Le Cs,Sb est un 
semiconducteur de type p. Cela résulte des mesures 
du coefficient de température de la force thermo- 
électromotrice et de la mesure de l'effet Hall par 
Sakata [16]. Cela est confirmé aussi par les mesures 
de l'influence de la composition stéchiométrique 
sur la conductibilité électrique par Sommer et 

Borzyak [101]. 

_ La dépendance de la température de la conduc- 
tbilité électrique, permet de déterminer l’énergie 
d'activation ÆA1 des centres p, ou éventuellement 
l’énergie AEm dans le cas de conductibilité intrin- 
sèque. Les résultats de Lukjanov, Mosover ; Juze, 

 Mochan, Ribken ; Schaetti, Baumgartner ; Suhr- 
man, Kangro; Harper, Choyke [17]; Mor- 
gulis [11] ; Sakata [12, 16] ; Wallis [10] sont repor- 
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tés sur la figure 6 en fonction de l’intervalle de 
température dans lequel ils ont été mesurés. 

On observe : 

1. Une grande fluctuation dans les valeurs de 
l'énergie d'activation. Ceci pourrait être (au moins 
partiellement) expliqué par une variation de la den- 
sité VA des centres d’impuretés. Pearson etBar- 
deen [18] ont montré en effet que l’énergie d’acti- 
vation des centres n et p diminue d’environ 0,1 eV 
quand V, augmente d’un facteur 100, par suite de 
l'interaction entre ces centres. Les cassures des 
courbes In o/T" (r = 3/4 ou 3/2) en fonction de 1/7 
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permettent de distinguer trois groupes d’énergies 
d'activation : 0,2 eV æ 06 eV, = 1,2 eV. 

On constate aussi : 

2. Une grande fluctuation de la valeur de la 
conductibilité électrique, qui provient d’une fluc- 
tuation des valeurs de la densité VA des centres 
d’impuretés ou de la mobilité des charges pou n 
(ce qui peut encore être au moins en partie une 
suite des variations de V;). La dépendance de la 
conductibilité électrique de la densité V4 semble 
bien résulter des mesures de Sommer ; Borzyak ; 
Wallis [9, 10]. 

De la mesure de Widmaier et Engstrom [19] on 
peut conclure que la conductibilité électrique est 
largement influencée par l’illumination de la photo- 
cathode Cs,Sb. Il s’agit probablement encore d’élec- 
trons captés dans ces centres p ou au fond de la 
bande de conduction, suite à l’illumination. 

Le travail d'extraction thermoionique ® mesuré 
par Wright [20] ; Sakata ; Eckart [12]; Wallis [9] 
et d’autres [15], [16], [17] fluctue également entre 
1,10 eV et 1,50 eV, mais chaque mesure donne 
une valeur unique de ®m. Pour un semiconducteur 
de type p on aurait Dm — AEm + W = E/2 et 
pour un semiconducteur intrinsèque 


DPrn — (AEm + W)/2. 


Les fluctuations des valeurs Om proviennent 
probablement encore des fluctuations de VA puisque 
Sakata trouve une augmentation de Dm quand Æ1 
diminue et que Pyatnitski [15] trouve 
Pr = 1,15 eV pour un rapport «optimum » de 
Cs/Sb. 

La confrontation de ces résultats expérimentaux 
montre surtout l’insuflisance des données actuelles 
pour une compréhension complète des phénomènes 
électroniques dans le Cs;Sb. Il serait particuliè- 
rement intéressant d'analyser l'influence de la den- 
sité VA et de l’illumination et de réaliser des 
mesures simultanées des différentes grandeurs me- 
surables. 


Modèles électroniques proposés pour le Cs,Sb. — 
a) LE MODÈLE DE WRIGHT ET BURTON. — Le: 
modèle proposé par Wright [20] et Burton [11] est 
schématisé à la figure 7a. On mesure hv,, AE, Dm, 
E\ et on suppose AE — AEm. Les désaccords 
éventuels entre ces valeurs sont expliqués par un 
réarrangement après photoexcitation (Wright). 
L'apparition des électrons lents pour Av > 4eV est 
expliquée comme une émission secondaire par 
Apker, Taîft et Dickey [12]. 


b) LE MODÈLE DE Miyazawa. — Le modèle de 
Miyazawa [12] fait intervenir une double branche 
pour le diagramme e(k) de la bande de valence, 
comme lindique la figure 7b. La photoexcitation 
se fait par transitions verticales. [1 ÿ aurait un 
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effet photoélectrique de surface avec le seuil Avset 
un effet photoélectrique de volume avec le seuil Av, 


Les électrons lents pour Av > 4 eV seraient dus à 
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un second effet photoélectrique de volume avec le : 
seuil Av; pour la deuxième branche de la bande de 
valence. 


Pie. 9, —@ Wright-Burton. — b. Miyazawa. — c. Piatniskii. — d. Apker, Taft, Dickey-Spicer. 


c) LE MODÈLE DE PYATNITSsKII [15] propose une 
courbure vers le bas de la bande de conduction et 
de la bande de valence suite aux états de surface du 
semiconducteur pur et aux atomes Cs absorbés à la 
surface dans le semiconducteur de type p. La zone 
de distorsion s’étendrait à environ 1 000 A. Le 
seuil photoélectrique est hv, mais l’effet photoélec- 
trique ne devient important qu’à partir de Av’ 
puisque la profondeur de sortie des électrons 
L & 250 À, d’après Burton [11]. 


d) LE MODÈLE DE APKER, TAFT et DICKEY ET LE 
MODÈLE DE SPICER. — Le modèle de Apker, Taft et 
Dickey [4] basé en fait sur l'hypothèse de transi- 
tions indirectes, a déjà été discuté dans la théorie 
générale. 

Le modèle de Spicer [11] est basé sur les mêmes 
hypothèses, mais introduit de plus l'influence de la 
profondeur d’émission photoélectrique 


L = Ja + B, 
déterminée par la profondeur de pénétration de la 
lumière (1/x) et par la profondeur de sortie des 
photoélectrons (1/8). Sa formule correspond à une 
approximation monocinétique dans un film d’épais- 
seur infinie : 


ilv) = ex GLL. 
Mais il introduit alors les simplifications sui- 
vantes : 
L = 1]x + Bæ1/x + Cte 
œ  (hv — hv,)*? au lieu de &° — (hv — hv,)2 ou 8 
et G = Cte au lieu de G — gp où p — E varie 
rapidement avec hv. Ceci explique que Spicer ne 


pouvait comparer sa théorie avec les résultats 
expérimentaux pour le K,Sb et Na,Sb et que la 


comparaison n’est possible pour le Cs,Sb, le Rb,Sb 
et le (NaK),Sb, que si on admet un effet photo- 
électrique d’impuretés, très important près du 
seuil. 

La diminution de l’efficience avec la température 
dans cette région n’est d’ailleurs pas une preuve 
d’un effet photoélectrique d’impuretés, puisque la 
dépendance de la température de AE en donne une 
explication suffisante. 

On peut proposer d’autres modèles encore. 
Nottingham [21] propose par exemple une jonc- 
tion p-n à la surface. Si nous admettons un modèle 
semblable à celui de Wright et Burton en suppo- 
sant pour la bande de conduction un diagramme 
e(k) parabolique et en portant l’énergie E,, en 
tenant compte de lénergie kv correspondante, on 
peut estimer la forme du diagramme e(k) pour la 
bande de valence. Ce procédé montre que la bande 
de valence n’a pas de branche fortement remon- … 
tante, comme le suppose Miyazawa. 


L'effet photoélectrique fondamental et la pro- 
fondeur d'émission du Cs;Sb. — Sakata [12] a 
décrit la distribution spectrale de l’efficence Y 


pour leffet photoélectrique fondamental par la for- 
mule 


Y = 3,5 10—2{hv — 1,7)3 


qui correspondrait à des transitions indirectes per- 
mises. Mais comme la courbe de l’absorption fonda- 
mentale coïneide mal avec les formules de l’absorp- 
tion par transitions indirectes, mais bien avec les 
formules de l’absorption par transitions directes, 
verticales, non permises, nous admettrons celles-ci, 
d'autant plus que la courbe Y peut être influencée 
par une dépendance de la fréquence de la pro- 
fondeur d'émission L. et que l’intervalle spectral 


“ 
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est trop large pour admettre des transitions indi- 
- rectes. 

Nous connaissons AE = 1,6 eV et hv, = 1,8 eV, 
ce qui permet de calculer p en supposant que les 
courbes e(k) de la bande de conduction et de la 
bande de valence ont leurs maxima pour la même 
valeur de À. p représente l'efficience maximum à 
laquelle on peut s'attendre et détermine en fait les 
ordres de grandeur (fig. 3 et 8). « est connu égale- 
ment (fig. 2). L’efficience photoélectrique est 


= ilel, = apgL(i — e—d/L) 


- ou simplement Y — opgL si d > L ou Y — xpgd 
si d & L. Pour des films très minces on aurait donc 
Y — ap. La distribution spectrale de Y mesurée 
expérimentalement s’écarte de «p et permet d’esti- 
mer la dépendance de la fréquence de gZL et donc 
de Z, si on suppose des films épais (fig. 8). 
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On trouve que L = 1 À près du seuil photoélec- 
trique et que L = 10? À plus loin du seuil. Comme 
la profondeur de pénétration de la lumière est 
{fa = 10% a 104 À et varie en sens contraire, la 
variation de L est bien un effet des diffusions des 
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électrons dans leur mouvement vers la surface. 
Dans les métaux, la profondeur d'émission L varie 
en sens contraire : L = 10% À près du seuil et 
L = 1 À plus loin [1, 2]. Mais dans les métaux la 
profondeur d'émission est déterminée par la diffu- 
sion (interaction individuelle) de l’électron excité 
avec les autres électrons de la bande de conduction 
ou par une excitation d’ondes de plasma (inter- 
action collective) pour des énergies plus élevées. 
Dans les semiconducteurs et isolants les inter- 
actions électron-électron ne sont possibles que pour 
hy = ? À En. 

Pour hv < 2 A Em seul les diffusions de l’électron 
excité avec des phonons ou avec les centres d’impu- 
retés sont possibles. Mais la diffusion par phonons 
ne change pratiquement pas l’énergie des électrons 
et le libre parcours moyen est indépendant de son 
énergie [7]. Par contre, la diffusion par les centres 
d’impuretés est une diffusion coulombienne, pour 
laquelle le libre parcours moyen croît fortement 
quand l’énergie augmente, puisque do = dQ/Q? E 
si E est l'énergie de l’électron et Q sa perte d’énergie. 
Comme les faibles pertes d’énergie sont les plus 
probables, les photoélectrons subiront un grand 
nombre de diffusions, avant de ne plus pouvoir 
sortir du solide. Notons que les-électrons très lents 
peuvent aussi être capturés dans les centres d’impu- 
retés. , 

Pour Av > 2 AEm des électrons secondaires 
peuvent être excités de la bande de valence à la 
bande de conduction et pour 


IN OIN ES EMNErCE 


on pourra observer une émission secondaire. La 
moitié des électrons lents observés pour Av > 4eV 
dans le Cs,Sb seront des électrons secondaires. 
L'émission secondaire explique aussi pourquoi 
Y = Cte pour hv > 4 eV, bien que «p © Cte et 
que Z diminue, comme dans le cas des métaux. 
L’excitation de plasmons demandera sans doute 
une énergie beaucoup plus élevée (kv > 14 eV). 

Quant à la dépendance de la température de 
l'efficience photoélectrique, on remarque que près 
du seuil Y © a(hv — hy;)" et que 


hvy AE = AE — cT, 


dont il résulte un coefficient de température 
K — dYJdT posiüf et tendant vers l'infini si 
hyv tend vers hv,. Par contre pour kv > 4 eV, 
l’émission secondaire devient prédominante et 
celle-ci a un coefficient de température négatif 
[14]. 

Nous remercions M. le Professeur de Hemptinne 
pour l'intérêt qu’il a pris à ce travail. 


Manuscrit reçu le 10 février 1961. 
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LE SPECTRE D'ÉMISSION D'UN ÉLECTRON NON-RELATIVISTE LIBRE DANS UN GAZ IONISÉ 


Par Lours GALATRY, 


Laboratoire des Hautes Pressions, C. N. R. S., Bellevue (Seine-et-Oise). 


Résumé. — Une théorie spectrale par intégrale de Fourier du rayonnement d’un électron libre 
est proposée. Après des applications au rayonnement de freinage et au rayonnement cyclotronique, 


on examine le cas où les effets de friction dynamique et 


enroulement des trajectoires autour des 


lignes de force du champ magnétique ne sont plus négligeables. 


Abstract. — À Fourier integral theory of the radiation by a free electron is proposed. After 
applications to bremsstrahlung and cyclotron radiation, cases are studied where dynamic 


friction effects and rotation around the lines of force of a magnetic field are not negligible. 


Introduction. — On connaît deux modes prin- 
cipaux d'émission par un électron libre, l’un est dû 
au freinage de l’électron au moment d’un choc avec 
des ions ou des atomes neutres ; l’autre se mani- 
feste lorsque, dans un champ magnétique H, l’élec- 
tron est accéléré par la force e/c (v X H) (+ est 
la vitesse de l’électron). Le premier type d’émis- 
sion (par transition hyperbolique) a été étudié en 
détail dans le cas où l’électron émetteur ne subit 
que des chocs isolés. Denisse [1], puis Oster [2] (*) 
ont donné des expressions de la répartition spec- 
trale émise dans le domaine des grandes longueurs 
d'onde. 

Dans le premier cas le traitement consiste à ane- 
lyser le rayonnement émis lorsqu'un électron décrit, 


Fic. 1. — Trajectoire d’un électron dans le champ d’un ion. 


(*) On trouvera dans la référence [2] une bibliographie 
plus détaillée. 


par rapport à un ion supposé infiniment lourd, une 
trajectoire hyperbolique (fig. 1) ; si on suppose que, 
lors d’un tel choc, tout le processus d’émission est 
concentré dans le temps p/v que passe l’électron au 
voisinage de l’ion déflecteur et si, de plus, la statis- 
tique sur les paramètres de choc p est faite en 
supposant qu’au delà de p = prindépendant de », 
le choc n’est plus optiquement actif, on obtient 
alors une distribution spectrale essentiellement 
constante dans les grandes longueurs d’onde. 
Denisse [1] a discuté la valeur de pr suivant les 
conditions physiques et Oster [2] à examiné les 
modifications à apporter à l’intensité du fond 
continu quand le processus d’émission n’a plus une 
durée courte devant la période du rayonnement 
émis. 

Le second type d'émission envisagé plus haut 
fournit, dans l’approximation non relativiste, une 
raie spectrale à la fréquence angulaire w. — eH [me 
qui est la fréquence de giration de l’électron autour 
des lignes de force du champ uniforme H [31 
Lorsque l’électron est isolé, cette raie possède la 
seule largeur naturelle. Si, par contre, on fait inter- 
venir pour calculer son profil l’action perturbatrice 
des autres particules du plasma, on obtient une lar- 
ceur supplémentaire. Pour des actions interparti- 
culaires schématisées par des chocs de durée infi- 
niment courte qui provoquent seulement de 
brusques changements dans les caractéristiques de 
la trajectoire hélicoïdale de l’électron, la raie pos- 
sède un profil de Lorentz de demi largeur égale à la 
fréquence moyenne des chocs +57. De plus l’aniso- 
tropie du mouvement des électrons entraîne une 
variation de la forme de la raie suivant l’angle 
d'observation par rapport à la direction du champ 
magnétique [4]. 

Les deux effets présentés ci-dessus seront traités 
ici en utilisant la théorie de l’intégrale de Fourier 
des raies spectrales qui est d’usage courant pour 
l’étude d’un électron lié et qui sera étendue au cas 
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d’un électron libre (première partie). Dans la 
seconde partie, on examine comment cette formu- 
lation permet d’aborder l’étude du spectre d’un 
électron libre dans un plasma dense, éventuel- 
lement soumis à l’action d’un champ magnétique 
uniforme. 


PREMIÈRE PARTIE 


Théorie classique de l’intégrale de Fourier pour 
un électron libre ; applications simples. — 
a) L'énergie rayonnée par unité de temps dans l’in- 
tervalle de fréquence dv — d(w/2x) par une 
charge — e décrivant une trajectoire r — r({) est 
donnée par [5] : 

+T .. 2 
die “lim el f de cie) | EU Ta) 
avec x — 2e2/3c3 ; rest le vecteur accélération de 
la particule, égal à m1 F, où F est la force 
s’exerçant sur elle à l’instant {et m sa masse. 

Par transformation du carré du module en une 

intégrale double, on obtient à partir de (1) : 
1 


F4 3 
es ex 


+T +T 
f- ue d'E(t).F(r)eiet—0 | av. (2) 


Dans la plupart des cas la force F possède un 
caractère purement aléatoire ; il est impossible de 
connaître sa forme analytique, donc de calculer les 
intégrales ci-dessus. Mais on peut remarquer que 


seul nous intéresse, pour comparaison avec l’expé- 


rience, le rayonnement moyen par rapport à toutes 
les trajectoires électroniques possibles entre — T 
et + T. En notant cette moyenne par < >, la 
relation (2) nous donne : 


< def: AT 
da > = À im 


4 1 di 
Es af, dé < F(D.F() > ciot—nl av. (3) 


Si l’on remarque qu’à cause du caractère aléatoire 
de la variation de F, la quantité < F(t).F(1') > ne 
dépend que de { — t = +ret non pas de 4, il vient 
immédiatement, après intégration sur L : 


LENS 
a fe < F(0).F(r) > cior drdv. (4) 


m? 
La répartition spectrale est done proportionnelle 
à la transformée de Fourier de la fonction 


of) = < F(0).F(r) > (5) 
dite fonction de corrélation de la force F (pour 
obtenir le coefficient d'émission par unité de 


volume, 1l faudra multiplier @ par le nombre n, 
d'électrons par unité de volume). 
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La relation (4) n’est autre que la relation fonda- 
mentale de l’analyse harmonique généralisée (théo- 
rème de Wiener-Khintchine [6]), appliquée au cas 
où la fonction aléatoire à analyser y(t) est le vecteur 
force subie par un électron. au cours de son mou- 
vement. La définition de la fonction « densité spec- 
trale » : 


G(f = Jim 


imp lAUNE: ta = [UE art est 


qui, d’après la théorie est la transformée de Fourier 
de la fonction de corrélation de y(t), coïncide préci- 
sément, dans sa forme analytique, avec la rela- 
tion (1) reliant l’intensité émise par unité de temps 
dans l’intervalle dv à l’accélération r de la particule 
émettrice. 
b) L’exemple le plus simple d'application de la 
relation (4) est relatif à un électron tournant autour 
d’un noyau d’un mouvement circulaire et uniforme 
de vitesse angulaire «,. On le suppose de plus. 
soumis dans son mouvement à des changements de | 
phase brusques et aléatoires représentant l’in-… 
fluence des chocs avec les atomes perturbateurs. 


Soit ®(t) la somme de ces déphasages au temps £. 


La force mise en jeu dans un tel mouvement est. 
lPattraction du noyau — Ze? r /r% qui est égale, pour 
une trajectoire stable, à la force centrifuge « r. 
On a donc 


F(0).F(<) — offr/? cos (wo, + + D(r) — D(0)). 

D'où 

< E{0).F{r) >-— 04 r2 
{ cos wÿ + < cos (D(r) — D(0)) > 
—sin 64 + < sin (D{r) — D(0)) > }. 

Sachant que, si les chocs détruisent complè- 

tement la corrélation de la phase, on a [7]: 
< cos (D{r) — D(0)) > = exp (— +frr) ; 
< sin (D(+) — D(0)}) > = 0 

il vient, après introduction de (5) dans (4) : 
2% 


2 


4 »2 
O5 ? 


CAUSE mn? 
_—1 —1 
TL, I TL 
: A = —2 
(© + @0)* + 7x 


dv (6) 


la) Pre 


qui est la forme spectrale de Lorentz bien connue. 
c) La relation (4) permet également de retrouver 


lexpression du coefficient d'émission par transi- 


tion hyperbolique. A cet effet on posera 
FE D Fr; 
ï 
ou Fi est la force exercée sur l’électron par le i-ème 
ion: d'où: 


plr) = 2 < F;(0).F;(r) >. 


dj 
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Si l’on néglige les corrélations entre les positions 
des différentes particules les termes où à 7 sont 
nuls et il reste : 


o(T) = nr Es F;(0) F;(r) > (7) 


ni étant le nombre d’ions par unité de volume. Or 
la force F, ne s’exerce que pendant le temps passé 
par l’électron au voisinage de la particule 1. Done, 
contrairement au cas de l’électron lié, le produit 
F,(0).F,(7) est ici nécessairement nul pour + supé- 
rieur à une certaine valeur +, assimilable au temps 
nécessaire pour parcourir la portion AB de la tra- 
Jectoire (fig. 1) ; soit +, — p/v, où p est le para- 
mètre de choc et v la vitesse relative à l’infini. Par 
ailleurs la partie réelle du facteur exp (tw+T) appa- 
rassant dans la transformée de Fourier (4) a son 
premier zéro pour + = ty — 7/20. Nous avons 
donc 


7mlcu = (plv) Lofr) — po(mIkT)"/?. 


Même dans le domaine des microondes 
(à æ 0,1 cm) et pour des valeurs raisonnables de T 
et de p (105 0K et 10-56 cm), on a tu/tmy > 6.102. 
La fonction œ(r) est donc nulle pour des valeurs 
de + dix fois plus faibles que le premier zéro de 
exp (twr) ; on approchera donc (+) dans le calcul 
de la répartition spectrale, par une fonction de 
Dirac : 


o(r) = n; a Ô(r) (8) 
où a est l’aire totale de la fonction de corrélation 


exacte ; soit, en faisant en sens inverse les opéra- 
tions qui conduisent à l’équation (3) : 


= +00 +00 3 
a sf dr Elto).Flto + 7) dte- 
—CO —0CO 

La barre signifie une moyenne sur toutes les con- 
ditions initiales du choc (p, v). Le spectre d’émis- 
sion devient alors : 


2% EURE COr 
dP{v) ane ma [°° eiot ÿ(r) dr dv 


2% a 
— — nr madv (9) 
me? 


L'évaluation de l’intégrale a peut se faire par les 
méthodes usuelles et conduit pour un plasma d’hy- 
drogène au résultat suivant [1], [2] : 


Ghrnre m 22) 
V) — x p \ (D 
AP(Y = mem DE os À) RE a dv. (0) 


(ep est la distance de Debye et p, — e?/2kT). On 
peut remarquer que, dans ce formalisme, le résul- 
tat qui consiste en un fond continu d’intensité 
constante est directement lié à l’hypothèse d’une 
fonction de corrélation représentée par une fonc- 
tion de Dirac. Mais comme le rapport tn/ty est 
assez proche de l’unité, il peut arriver, si certains 
effets sont capables d'augmenter le temps de corré- 
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lation, que cette approximation ne soit plus valable. 
Ceci entraînerait une variation de lintensité du 
fond continu avec la longueur d’onde. 

d) Un dernier exemple d’application simple de la 
relation (4) est celui du rayonnement cyclotronique. 
Considérons un électron en mouvement dans un 
champ magnétique uniforme ÆH et subissant de 
plus des chocs avec d’autres particules ; ces chocs 
sont supposés instantanés et tels que la direction de 
la vitesse de l’électron après un choc soit statisti- 
quement indépendante de la vitesse immédia- 
tement avant le même choc. On a alors : 
< F(0).F(r) > = . < (e(0) x H) (6x) X H)>. (1) 

La moyenne dans le second membre peut se dé- 
composer de la manière suivante : 


SF (0). Fi) >= LT X H}.(< ü(r) > x H) 


où << v(r) > est la moyenne sur tous les v(r) 
possibles pour un t(0) fixé et où la barre indique 
une moyenne sur ®(0). En vertu de lhypothèse 
faite sur la nature des chocs, < v(+) > — 0 dès 
que lintervalle (0, +) contient au moins une colli- 
sion. Ainsi la fonction de corrélation cherchée sera 
égale à : 
NCAA) — 


Solr) (v(0) X H).(v(T) X H) 


où ®(r) est la vitesse de l’électron dans l'hypothèse 
ou il n’y a aucun choc entre 0 et + ; elle se déduit 
de (0) par une rotation d’angle «, + autour de la 
direction du champ magnétique («x — eH mc). De 
plus #,(r) est la probabilité qu'il n’y ait aucun 
choc entre les instants 0 et +. En supposant valable 
ici, même en présence d’un champ magnétique, une 
distribution de Poisson pour les chocs, on a: 
Pot) entr Ilmientalorss 


< F{0).F{r) > — = v£(0) H? cos over 


12 


où 0, — 0 sin 0 est la projection de ®(0) sur un 
plan perpendiculaire à H. Or e2H?/[m°?c? = «£ ; on 
obtient donc la répartition spectrale suivante : 


a _— CA 
dP(v) = 2rxo8 02 R eiOT e—T/Tr, COS we + dr dv. 
0 


Ou encore, en prenant une distribution maxwel- 
lienne des vitesses : 


2kT 


dP(v) = 2kw8 = 


71" ne 
| en 5 bu ZE =) GI 2) 
(OMR MO ot Es 


Ce résultat est identique à celui obtenu par utili- 
sation du théorème de Parseval [4]. 
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DEUXIÈME PARTIE 


Comme on l’a vu précédemment, la relation (4) 
permet de retrouver rapidement les résultats essen- 
tiels de la théorie de l’émission cyclotronique et de 
l’émission par transition hyperbolique. On peut 
espérer que son caractère de généralité permette 
d'obtenir des indications valables lorsque la théorie 
élémentaire de ces deux effets ne l’est plus. 

Considérons un électron en mouvement dans un 
plasma, soumis de plus à l’action d’un champ ma- 
gnétique extérieur uniforme H. La force qui 
s’exerce sur cet électron est donc : 


Po Fe 


où F, est la force coulombienne due aux charges 
des autres particules et Fu la force due au champ 
magnétique (e/c(t X H)). La fonction de corré- 
lation du problème s’écrit : 


* o(r) = < (Fol0) + Fx(0)).(Folr) + Fa(r)) >. (13) 


La relation (4) indique alors que, d’après (13), le 
spectre de l’énergie émise par unité de temps est la 
superposition de deux spectres distincts obtenus 
par transformation de Fourier des deux fonctions 
de corrélation suivantes : 


1 = < Fç(0).Folr) > 


De = < F0) .Fa(r) >. M) 


Ces fonctions se rapportent respectivement au 
rayonnement de freinage (modifié par la présence 
de H) et au rayonnement cyclotronique (modifié 
par la présence du champ coulombien). Les termes 
rectangles, qui n’apportent pas de contribution 
appréciable à l’intensité totale, n’ont pas été men- 
tionnés. 


I. Érupe DE @,(r). — On sait que dans les mi- 
eux ionisés l’influence des interactions multiples, 
à longue distance, est au moins aussi importante 
que celle des chocs binaires, à courte distance. 
Deux points de vue peuvent donc être adoptés pour 
une étude qualitative de +, ; l’un utilise la notion 
de choc, donc surtout valable pour les faibles 
valeurs de +, telles, par exempie, que rv << p, où 
Po = Ze?[mv? est le paramètre de choc conduisant, 
pour un choc binaire électron-ion, à une déviation 
de x/2 ; l’autre utilise la notion de champ interne ; 
on en tiendra compte pour des valeurs de + supé- 
rieures aux précédentes. 

a) Examinons d’abord linfluence des « chocs », 
c’est-à-dire des approches électron-ion telles que, 
par exemple, la distance d'approche devienne infé- 
rieure à PQ, ; On supposera de plus que, si à l’ins- 
tant O l'ion le plus proche de lélectron est l’ion 
numéroté 1, alors l’électron subira l’influence de 
nombreux autres ions avant de sortir de la sphère 
d'influence de cet ion n° 1. Ceci implique que 
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Po <-pD Où on —6,9 (Tn.)” est le rayon del 
Debye caractérisant le milieu ionisé. On se placera. 
ensuite dans l’hypothèse d’un plasma d'hydrogène 
considéré comme un gaz de Lorentz (ions immobiles 
de masse infinie, charge des électrons de champ 
uniformément répartie dans l’espace). On prendra 
comme origine des coordonnées le centre de l’ion 1. 
A l'instant 0 on peut considérer que seul cet 1on 
agit sur l’électron étudié puisque celui-ci se trouve 
près de la singularité du potentiel e/r, (autrement 
dit le fait prépondérant à l’instant O0 est le choc 
avec l’ion n° 1). On a donc : 


Fc(0) — F;(0). 


De plus la seule fraction de F.(r) pouvant avoir 
une corrélation avec F,(0) est constituée par la 
force F;(x) et par la somme des forces F(r) dues 
aux particules se trouvant dans la sphère de Debye 
centrée en 1. On a ainsi : 


PralT) = < F;(0) ») Fr) > = < F;(0) Fip(r) 2 (15) 


où la barre indique une moyenne sur toutes les 
positions possibles des ions dans la sphère de Debye 
de 1. L’expression et F,,(r) est le champ de 
Debye-Hückel : 


F1 = e?rifr$ exp (— rilep)} 4 + ri1lep). 


La moyenne figurant dans le membre de droite 
de l’équation (15) ne porte plus que sur les positions 
possibles de l’électron émetteur (électron-test) à 
l'instant +, soit r, + s, par rapport à sa position r; 
à l’instant 0. Si on introduit la probabilité 

P(S ; vo, 7) ds 
pour que l’électron, animé à l’instant 0 d’une 
vitesse (0) subisse un déplacement compris entre s 
ct s + ds pendant la durée + et la probabilité 
P(r;) dr, pour qu'il se trouve repéré à l’instant 0 
par un vecteur compris entre r,;etr, + dr,,on a : 


gas) = f Pins) dis P((0)) de(o) 

fPis: vb), à Fr) Fiptr, + s)és (16h 
P{&(0)) est la fonction de distribution de la vitesse 
de l’électron à l'instant 0. On a de plus : 

P{r,) = ne nj exp (— V,JkT) 
où V, est l’énergie potentielle de l’électron placé 
cn +,. On la négligera et on posera : 

P(r;) = ne ni, 
Pour les petites valeurs de + (c’est-à-dire telles 

que + < p,{v) on négligera la dispersion de r, +s 
par rapport à la valeur calculée r, + (0) + en 


l’absence d'ions. D’où P(s ; v(0), +) = à(s — rv(0)). 
Dans ces conditions il vient : 


ua = ne m6 | PO) de(o) [Fi] Fins + re(0)) as 
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Le calcul de cette intégrale peut être assimilé 
(sachant que [F2] < [F;] et moyennant les mêmes 
considérations sur le domaine spectral étudié : 

. © < v/po) à celui de la fonction de corrélation (8). 
On a donc : 


Dia = Ne M GOT). (17) 


Pour les valeurs de + plus grandes, mais néan- 
moins suflisamment faibles pour que l’électron 
n’ait pas de fortes chances d’être sorti de la sphère 
de Debye au temps + (+ < p»n/v), on peut faire 
_ intervenir, dans (16), par l'intermédiaire de 
… P(s ; v(0), +), le freinage du mouvement de l’élec- 
tron par les particules situées à l’intérieur de la 
sphère de Debye de l’ion 1 ; ceci peut se faire 
approximativement en posant 


P(S ; v(0), <) = (4x Dr)" exp (— [s/*J4D+) (18) 


où D est le coefficient de diffusion de l’électron au 
sein du plasma. D sera pris égal à AT ]mB [8], 6-1 
_ étant le temps de relaxation généralisant pour le 
cas de forces à longue distance la notion de fré- 
quence moyenne des chocs. En fait (18) n’est qu’une 
approximation pour la densité de probabilité 
- Wa(s ; vo, +) d’un déplacement s au temps 7, 
connaissant la vitesse v(0) au temps 0. L'influence 
de ce processus ne se fera sentir de manière appré- 
 ciable que si, du fait du freinage, l’électron reste 
sous l’influence de l’ion 1 pendant une durée de 
l’ordre de grandeur, ou supérieure, à t/2w, c’est-à- 
dire si la distance moyenne / parcourue pendant le 
temps 7/20 est de l’ordre de, ou plus petite que 
le rayon de Debye. Soit, d’après le résultat classique 
d’Einstein : 


1 =V2Dr 2 pn = 6,9 (T{ne)/? 


En posant ci-dessus + — 7/20, on obtient la 
condition suivante, donnant la limite inférieure du 
domaine spectral dans lequel cet effet pourrait être 
observé : 


NS 100 cT (Dn,)\r 


(à titre d'exemple, on aurait, pour T = 1040K et 
me — 1012 em—$, À < 10=% cm). On constate que 
. pour les conditions physiques choisies, cette Inéga- 
lité exclut le domaine des ondes centimétriques et 
millimétriques. Néanmoins, quand elle est réalisée, 
la fonction de corrélation comporte, en plus de la 
fonction de Dirac (17), le terme supplémentaire : 


sos) = neni f Pots, =) &s | Fi(r:) 
L'intégrale en r, sera calculée en assimilant, 
comme plus haut F,, à F,. On obtient alors [9] 


[rit 


).F,p(r, + S) dr;. (20) 


).Fi(r, + 8) dr = 4re [sf 


d’où Ê 
P1d = Ne M AV D Avr, (21) 
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D’après l’équation (4) la contribution de 6,4 à 
la distribution spectrale s’écrit alors : 


DA F 
aP(y) = 5% [7 cios(ose + qu) dr dv 


te Hein tt ve dv (22) 


ou Zx est l’intensité du fond continu donnée par 
l’équation (10). Le fond continu tend donc, d’après 
cette correction, à diminuer d'intensité vers les 
grandes énergies. 

Dans le cas où il existe un champ magnétique 
uniforme EH, les électrons s’enroulent autour des 
lignes de force de H et la diffusion électronique 
dans la direction perpendiculaire à Æ est très di- 
minuée [10], [11]. L'effet de cette nouvelle entrave 
apportée au mouvement des électrons ne sera 
évidemment sensible dans le spectre que si le rayon 
de gyration bG est plus faible que le rayon de 
Debye op, au delà duquel les corrélations dyna- 
miques s’éteignent ; d’où la condition : 


H> (6,9)-1 2% (S} 


Ne 


(23) 


(à titre d'exemple on peut voir que si nr, — 101? 
em, ceci exige À > 2 500 gauss). Si la condi- 
tion (23) est remplie et si, à l’instant 0, la vitesse 
v(0) de l’électron est sensiblement perpendiculaire 
à H, celui-ci restera dans le champ des ions situés 
dans la sphère de rayon p» + p& et centrée sur sa 
position moyenne pendant un temps beaucoup plus 
grand que p,fv. Dans ce cas aussi l’approximation 
de la fonction de corrélation par une fonction de 
Dirac n’est plus valable, et le problème est plus 
complexe car, d’une part lémission du rayon- 
nement devient anisotrope et, d’autre part le mou- 
vement se fait, avant que l’électron ne s’échappe de 
la sphère d’action de lion {, sur une trajectoire 
approximativement hélicoïdale non centrée sur 
lion n° 1. On prendra comme exemple le cas d’un 
plasma d'hydrogène totalement ionisé soumis à un 
champ magnétique constant Æ —3 000 gauss et où 
le rayon de Debye est 

pp = l0 Sen. (ler OSOR ane 10 en 
Dans ce cas un électron de vitesse 3.107 cm.s-1 
perpendiculaire à H décrit un cercle de rayon 
ea —D,7.10-# cm. Cet électron n’aura une proba- 
bilité notable d’avoir parcouru par diffusion la dis- 
tance on qu'après le temps & — 10-6/4D. En esti- 
mant d’après [10] D à 0,6.10-* cm? s-1 pour la 
direction perpendiculaire au champ magnétique, on 
obtient t —0,4.10-1$, Dans le domaine des ondes 
centimétriques, d'autre part, l’approximation de 
Dirac n’est plus valable dès que 

i 


EE 


C’est donc bien le cas pour les conditions expéri- 


DR 
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mentales ci-dessus et cet exemple montre qu'il 
serait alors nécessaire, pour le calcul du spectre, de 
tenir compte de l’effet de confinement des électrons 
par un champ magnétique dans la sphère de Debye 
d’un ion donné. 

b) Passons maintenant à l’étude, dans le cas 
I — 0, de l’influence exercée sur @(r) par la pré- 
sence des particules ionisées du plasma considérées 
maintenant non comme susceptibles d’entrer en 
collision avec l’électron-test, mais au contraire 
comme génératrices d’un champ interne (somme 
de champs individuels) s’exerçant sur de grandes 
distances (de l’ordre de la distance de Debye) et de 
gradient relativement faible. On sait que ce champ 
interne est capable de dévier la direction de la 
vitesse de l’électron. La variation par unité de 
temps des composantes moyennes de ® est connue 
sous le nom de coefficients de diffusion [3] ; la varia- 
tion de © permet de plus de définir un temps de 
relaxation +, au bout duquel le vecteur vitesse, 
régi à l’instant O par une distribution micro- 
canonique, peut être considérée comme statis- 
tiquement indépendant de (0) et régi alors par une 
distribution maxwellienne. On a [3] pour un gaz de 
Lorentz : 


m°? 3 É =) ve 


ZE —————————— ———— Se À —= -e 
 — srme log À ; 2° 


T le 

Puisque l’électron subit ainsi une altération dans 
sa vitesse, on peut considérer qu’il est soumis de la 
part des autres particules (la notion de choc étant 


ES 

exclue) à une force stochastique (+) dont on 
ignore l’évolution dans le temps, mais dont on 
connaît la valeur moyenne. Cette dernière peut se 
déduire du travail de Chandrasekhar [8] selon 
lequel l’accélération (ou plutôt la décélération) 
moyenne — Av de l’électron pendant un temps At 
crand devant la durée d’un choc (At > p,fv) mais 
plus petit que le temps de variation appréciable des 
grandeurs macroscopiques caractérisant le milieu 
est relié à sa vitesse par la relation : 

20 go 

Le coefficient B, dit coefficient de friction dyna- 
mique, peut être considéré [8] comme l’inverse du 
temps de relaxation 7, de la direction de la vitesse 
de l’électron émetteur. La force moyenne subie par 
électron du fait des actions à grande distance 
s’écrira donc : 

FU) — — mp(v(t)) v(t). 

On conçoit que cette force, dépendant du vecteur 
vitesse de l’électron, ait des propriétés de corré- 
lation différentes de celles des forces coulombiennes 
examinées plus haut. Ces propriétés seront rapi- 
dement étudiées par l’intermédiaire de celles de la 
fonction : 


Te 
Quel) = <F(0).F(r) > 


= m2 < B{v(0)) B(v(r)) v(0).v(r) >. 


(25) 
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On supposera que [v| varie peu avec £ et on. 


: v(t anse 
adoptera pour la corrélation de a une décrois- 
sance exponentielle conditionnée par le temps de 
corrélation +. D’ou : 


quels) — m? B2{0(0)) J0(0)F er. 


La moyenne indiquée par une barre ne porte 
plus que sur (0). En portant dans (4) on obtient la 
distribution spectrale suivante : - 


b 
CES 


C’est une moyenne sur des distributions de type 
Débye dont le maximum (pour & — 0) a la valeur : ” 


dP(v) = 2x B#{v(0)) [v(0)[? 


£r n; et log A 


dP(0) = 2xB(v(0))[v(0)|?dv = 2xr me v(0) dy 
8ret ]m 
re 24e 4 MERTL NELCA dy (26) 


De même que la fonction ow, la fonction x ne 


peut pas être assimilée à une fonction de Dirac car 


le temps de corrélation mis en jeu r — 8-1 n’est 
pas nécessairement petit devant la période de la © 
radiation émise. Par exemple pour : 


Ne = Ni —= 1077 CE 


T — 104 0K et dans la région des ondes centi- 
métriques on a d’après (23) % = 0,5.10$8 s et 


r/20—10-11s. Remarquons que, contrairement à … | 
ce qui a été vu plus haut cet autre effet de corré- n : 


lation n’est pas dû à l’extension spatiale du do- 
maine d’action des forces en présence. En effet les 
ions agissant sur l’électron à l’instant 0 sont diffé- 
rents de ceux agissant à l’instant + dès que l’élec- 
tron a franchi une distance 26, c’est-à-dire dès que 
2op [v < +, par exemple pour les mêmes conditions 
physiques que ci-dessus dès que + > 4.10-11s, 
c’est-à-dire dès le millième du temps de corré- 
lation +,. Le présent effet, au contraire, est dû à la 
persistence d’une grandeur physique attachée au 
systéme émetteur lui-même, à savoir la direction 
de la vitesse de l’électron. 


IT. ÉTUDE DE ®,(+). — D’après (14) on a : 


paf) = < (d(0) X H).(v&(r) x H) >, (27) 

Alors que la forme de w,(r) dépendait de la corré- 
lation de la force de caractère électrostatique agis- 
sant sur l’électron (donc, partiellement au moins, 
de sa position), on constate que celle de œ,(t) ne 
dépend que de la corrélation de sa vitesse entre les 
instants 0 et +. On adoptera les coordonnées de la 
figure 2, le champ magnétique étant dirigé suivant 
l'axe z et la vitesse ®(1) de l’électron étant repérée 
par les deux angles d’Euler 0(t) et œ(t). En l’absence 
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de toute particule perturbatrice, la fonction de 
corrélation est proportionnelle à : 


(Lx — Ur10) COS we T) D(Py — Py'0) sin «we +) 03 — 03/0). 
(28) 


Fic. 2 — Évolution du vecteur vitesse d’un électron 
par rapport à un champ magnétique uniforme. 


On a étudié dans la première partie la forme de 
o,(T) dans le cas où l’on suppose : 

a) que les chocs sont instantanés (rayon de gyra- 
tion très grand devant le rayon d’action des forces 
ioniques) ; 

b) que les chocs sont tels qu'il n’y a pas de corré- 
lation entre © avant le choc et v après le choc. 

L'étude du cas général où les hypothèses a) et b) 
ne serajent pas faites nécessiterait la résolution des 
équations du mouvement d’un électron dans le 
champ d’un ion, l’ensemble étant plongé dans un 
champ magnétique constant et uniforme d’inten- 
sité quelconque. On se contentera ici de lever la 
restriction b). 

Dans les coordonnéss de la figure 2 on a : 


@2(7) = v(0) sin O(U) < v{r) sin O(+) cos o(r) > H? (29) 


(la barre indiquant une pondération sur les condi- 
tions initiales). On se placera dans le cas d’un gaz 
de Lorentz et on admettra que [v(r)| varie peu au 
moment des chocs : )| 

calculer : 


» << sin O(r) cos œ(r) > 


CO CO 
— D; Fnlr) < sin 0 cos p >n — Ÿ Fi(r) at) 
0 0 


où #,(7) est la probabilité pour que lélectron 


subisse r chocs dans l’intervalle de temps (0, +) et 
< >" la moyenne de la quantité sin 0 cos © 
lorsque la durée + comporte n chocs. On a : 
n—1 
an) — Reiort < IT ciAv; sin 0, (cos Av, +isin Am) > (?) 
il 


où Ao; est la variation de l’angle ® lors du 1-ème 
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[o(r)| = |v(0)|. Il faut donc. 
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choc et 0, l'angle 0 après ce 1-ème choc. Si on dé- 
signe par < >>, une moyenne sur toutes les possi- 
bilités au n-ième choc (c’est-à-dire, puisque |v| est 
constant, une moyenne sur le paramètre de choc p), 
on obtient : 
n—1 
at = < IT eiâe,;(< sin 0, cos Aer >» 
£ + i < sin 06, sin Apy a) >. 


Les chocs étant supposés de durée très courte 
devant la période de gyration, on raisonera sur le 
choc en ignorant la présence du champ magnétique. 
Les déflections A, et — Aw, étant alors également 
probables on aura : 

< sin 0,,sin A on > = 0/* 


De plus on a la relation suivante entre les angles 
de la figure 3 : 


sin 0»— Sin 0, cos Aw, = CoS 4n — cos 0,4, cos Û» 


Où Xn est la déflection des tangentes à la trajectoire 
lors du n-ième choc. D’où : 


sin Op < sin 04 cos Am n 
C0 pin 00 0,22,.7c0s 0,124 0180) 


Fic. 3. — Déflection d’un électron par un ion supposé fixe 
{on n’a représenté que les tangentes aux trajectoires). 


Le dernier crochet peut être obtenu à partir de la 
relation : 


COS 0, — cos y. Cos0, "+ Sin, sin 0,2, cos LE, 


où [”, est l’angle dièdre entre le plan contenant H° 
et la tangente avant le choc et le plan analogue 
après le choc. Puisque l’azimut de l’ion déflecteur 
a une probabilité constante, par rapport à la tra- 
jectoire avant le choc, on aura : 


< sin {y COS Ty >n = 0 
d’où : 


NOUS CON COS Ye 
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et, en portant cette égalité dans (30), il vient : 
<'sin 04 c08 Les n — Sin Ones < COS Yn n: 


Nous avons donc : 


n—1 


an — < COS Yn n < IT aa) sin 0,1 >(#—1) 
1 


soit, par itération : 
an) = < COS {n nn < COS Xn—1 m1 + 

< cos x1 > 1 sin 0(0) 
et, puisque les chocs sont indépendants : 


an) — sin 0(0).(< cos x >}2. 
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Bien que % dépende en toute rigueur de H, on 
adoptera la distribution de Poisson 
Far) = (nt (r/rnt exp (— rÎrr) 


où +72 est la fréquence moyenne des chocs. Il vient 
alors : 


C9 
2 
ge = = 10(0)1* sin? 0(0) cos we + 2 M < cos x 7 H° 


ei T 
BE 2 —; (1—<cosXx>) pJ2 
na (0) cos Ter, H 


où o1 est la composante de v sur un plan perpen- 
diculaire à H. Le spectre correspondant à cette 
fonction de corrélation est alors, par transior- 
mation de Fourier : 


d 


2 — £ os y > 
ER NT ms 


m? c? TL 


1 — < cos x > 


a 


(oo =")? + 
Calcul de < cos 4 >.— Ona 


où Pn(v) est le paramètre de choc au delà duquel on 
considère que la notion de choc perd son sens pour 
un électron de vitesse v. Connaissant l’expression 
de cos 7 dans un choc électron-ion [3], on obtient : 


1— < cos x > = 2{efmvt)? pre log (1 + (mwu?fe?) p3). 


Si on définit, à titre d'exemple, la notion de choc 
à partir du moment ou l’électron subit une déflec- 
tion de x/2 (Ppn — e?[mv?), il vient alors : 


PA coS A0 op 


et, puisque cette quantité est indépendante de », 
la relation (31) se simplifie et donne, compte tenu 


112 COS y ef dy. 
TL 


ni 


du fait que, pour une distribution maxwellienne des 
vitesses #1? — 2AT ]m : 


(o + wc)? + 


2x e? ,,2kT 2 log 2 
Ce RE +) D aise ê 
dP{v) = ns RE 
1 1 1 
Lo — «ç)? + ( | (Oo Eee er SE | fa 
FL 7x, 


C’est donc seulement à partir de la définition un 
peu arbitraire de choc qui a été prise ci-dessus que 
l’on retrouve (au facteur 2 log 2 près) la forme spec- 
trale obtenue par la théorie plus élémentaire de la 
première partie (équation (12)). 
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MESURE ABSOLUE DES SECTIONS EFFICACES DES RÉACTIONS (n, p) ET (n, d) 
AVEC UN TÉLESCOPE A COINCIDENCES 


Par MM. P. LÉVY et BONNEL, 
Laboratoire de Synthèse Atomique et d’Optique Protonique du C. N. R. $. 


Résumé. — Afin d'étudier sélectivement les réactions (n, p) et (n, d) avec divers éléments, nous 
avons effectué un montage comportant un télescope de coïncidences conduisant à une différen- 
tiation des particules par l’observation des coordonnées d’un point lumineux sur un oscilloscope. 
Associée à une mesure absolue du flux de neutrons, cette installation permet d’atteindre lessections 
efficaces totales et différentielles de réactions nucléaires. La réaction F{n, d)'80 et la réaction 
de rupture du deuton sont données en exemple et les résultats comparés à des résultats antérieurs. 


Abstract. — In order to study selectively, (np) and (nd) reactions for various elements, we 
built an arrangement compuring ‘a coincidence telescope allowing differentiation of particles by 
observation of the coordinates of a spot on an oscilloscope. Together with absolute measures of 
neutrons flux, this installation allows us to measure total and differential nuclear reaction cross 
sections, #F{n, d)180, and deuteron break up reactions are given as illustrative examples, and our 
results are compared with earlier results. 


A l’époque où nous avons entrepris ce travail, un 
INTRODUCTION certain nombre de résultats avaient été obtenus 
sur les répartitions angulaires de protons soit par 
L’étude de la répartition angulaire des particules des caméras à plaques (avec inconvénient de mau- 
émises lors des réactions avec des neutrons de  vaises statistiques) soit avec un télescope à proton ; 
14 MeV jointe à une mesure précise du flux de ces expériences ne permettaient pas de déterminer 
neutrons incidents donne des renseignements inté- la nature des produits de la réaction et leurs résul- 
ressants. tats laissant supposer l’existence de réaction (n, d). 
La connaissance de 36/00 (E, 0) permet de dis- Nous avons donc essayé de réaliser une instal- 
tinguer la proportion d'interaction avec formation lation permettant d’analyser sélectivement les pro- 
de noyau composé de la proportion d’interaction tons ou les deutons dans une très large gamme 
d’énergie tout en mesurant de façon précise le flux 
de neutrons incidents. La photo n° 1 montre la dis- 
position des appareils de détection. 


I. Mesure du fiux de neutrons. — a) Nous avons 
choisi la méthode de la particule associée [1] qui 
consiste à compter les x de la réaction d(T, «)n — 
dans un angie solide donné. La direction a été 
prise à 900 pour laquelle la répartition en énergie 
des « est la plus étroite. Une légère correction est 
appliqué: pour tenir compte du fait que les neu- 
trons utilisés ne sont pas pris dans une direction 
d'émission associée à la direction des &« mesurés. 
Un certain nombre de précautions doivent être 
prises pour obtenir une grande précision et une 
bonne stabilité des mesures, en effet : 

a) La largeur du pic varie avec l’énergie des 
deutons incidents. 

Nous avons utilisé successivement le tube 300 kV 
construit par l’un de nous (PL) et le générateur 
directe, d'atteindre les températures nucléaires, et 900 KV du Laboratoire de Synthèse Atomique. 
dans certains cas d'interaction directe de déter- b) Le bruit de fond est fonction de la tension de 
miner l’énergie et le moment orbital des premiers fonctionnement et du réglage de la focalisation de 
niveaux du noyau résiduel, l'accélérateur, 
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c) Les produits des réactions parasites 
D(D, n)Hr,  D(D, pr 


qui ont lieu soit dans la cible elle-même soit dans 
son support sont présents dans le spectre ; leur 
nombre augmente en fonction du temps et varie en 
fonction des conditions de marche du tube accélé- 


rateur. (L'observation du spectre fournit d’ailleurs 


une méthode de réglage rapide du centrage et de la. 


focalisation du faisceau.) Il faut donc isoler le pie 
au maximum en réduisant le bruit de fond, en 
interposant le minimum de matière sur le trajet 


des « et encadrer ce pic par la fenêtre d’un sélecteur. 
La figure 1 bis montre la disposition des diffé- 


2 


Fi. 4 bis. — Coupe schématique du moniteur. 


1. Pièce porte cible. — 2. Circuit de réfrigération. — 3. Écran d’aluminium 
amovible (un sur cinq représenté). — 4. Photomultiplicateur. — 5. Blin- 
dage u-métal. — 6. Cristal. — 7. Écran d’aluminium (200 ug/em?). 


— 8. Masque amovible. 


Nb. coups /canal 


7500 
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Fie. 2. — Spectres obtenus avec le moniteur &. La courbe IT 
montre les pics par.sites dus à la réaction d(dp)T. 


rents éléments. La feuille de 200 ug/em? d’alu 
sert à éliminer les photons lumineux produits au 
voisinage de la cible. Aucun diaphragme ne se 
trouve sur le trajet du faisceau et l’angle solide est 
défini par la surface du cristal (9 — 4 mm, 
e — 1/10 mm) ce qui évite les diffusions parasites 
des & (pour la même raison les parois des pièces ont 
été éloignées du faisceau utile). 

L'interposition d’écrans d’aluminium calibrés 
permet d'identifier la nature des particules obser- 
vées par le déplacement des pics sur le spectre. La 
plaque effaçable située devant le porte cible sert de 
butée pour le positionnement de la cible, de masque 
permettant les régiages du tube sans émission de 
neutrons et de support de source pour les réglages 
du moniteur. | 

La figure 2, montre deux spectres types corres- 
pondant l’un à une cible neuve non chargée en 
deutons, l’autre à une vieille cible chargée en deu- 
térium. 

Si l’on utilise un cristal épais (2 mm), il n’y a 
plus de séparation nette entre le pic et le bruit de 


N° 8-9. 
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fond, et il apparaît d'autre part un pie corres- 


pondant au protons de 13,8 MeV de la réaction 
SHe(dp)*He ($He provenant de la désintégration du 
tritium). Une étude séparée en comparant les taux 
de comptage «&-protons simultanément mais avec 
des angles solides différents, a montré que le 
nombre d’'x produit par cette réaction parasite et 


SECTIONS EFFICACGES DES RÉACTIONS (n, p) ET (n, d) 


4 


491 


comptés simultanément sont suflisamment faibles 
pour ne pas mettre en cause la précision de la 
mesure, Ce phénomène serait toutefois à considérer 
pour des sources stockées un certain temps. 
Compte tenu des mesures précises de la géo- 
métrie du système à mieux que 0,5/00 on peut 
estimer la précision du moniteur à + 2,5 %. 


FRS 


mi 


ÈS 


q RSS) 
Hé jury] 


\ 
F1. 3. — Coupe du télescope. — I-II. Zones utiles des compteurs proportionnels. — 
III. Cristal d’ICs. — IV. Radiateur. 


1. Tube de champ. 


— 2. Photomultiplicateur. 


3. Boîtier du préamplificateur PM. 


— 4. Billage. — 5. Roue porte radiateurs. — 6. Boîtier. 


II. Construction du télescope. — Le télescop 
réalisé est dérivé du télescope construit par 


Johnson et Trail [2] — avec un certain nombre de 


modifications mécaniques facilitant le démontage. 
Le choix de la géométrie est assez arbitraire étant 
donné le nombre de paramètres entrant en jeu ; un 
compromis est nécessare entre l’efficacité et la 
dispersion angulaire. A efficacité et résolution angu- 
laire égales, nous avons préféré un télescope de 
petites dimensions qui facilite la préparation des 
radiateurs notamment dans le cas d’isotopes sépa- 
rés et diminue les fluctuations des impulsions dans 
le temps (voir $ 3), nous fixant comme objectif 
la meilleure séparation possible des protons et des 
deutons. Nous nous sommes orientés vers un téles- 
cope utilisant pour la coïncidence deux compteurs 
proportionnels coaxiaux dotés de tubes de champ 
fournissant simultanément l’information (voir $ 3) 
avec la meilleure définition possible. Par contre 
l'expérience a prouvée que la présence du tube de 
champ introduisant au bruit de fond important. 
Une coupe du télescope est indiquée sur la figure 3 
— ainsi que les éléments géométriques essentiels. 


FIc:0 


bis, — Le télescope partiellement démonté. 
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L’angle solide est de : 2 10—2 stéradians et pour 
une distance cible radiateur de 10 em adoptée dans 
les expériences l’ouverture angulaire maximum est 
de 13032: l’angle moyen de 59, pour une orientation à 
00. Le télescope tourne autour de deux pivots dont 
l’axe passe par le centre du radiateur. Une roue 
porte écran comporte cinq logements sur lesquels 
sont fixés une source de plutonium, un radiateur de 
polyéthylène ou de deutéro-paraffine, une cible 
épaisse de paraffine ou de deutéro paraffine (voir 
$ 3), le radiateur en étude et un support de 
radiateur nu pour la mesure du bruit de fond. La 
roue peut-être tournée de l’extérieur par l’inter- 
médiaire d’un passage étanche ; son position- 
nement précis est assurée par un billage (fig. 3 bis). 

Le cristal d’iodure de césium (épaisseur 1,5 mm) 
est couplé à un photomultiplicateur K 1382 par 
l’intermédiaire d’un guide de lumière en plastique 
qui assure l’étanchéité. Pour diminuer le bruit de 
fond il y a intérêt à garnir les zones des télescopes 
pouvant donner lieu à des triples coïncidences 
réelles avec un métal Dur. C’est pourquoi la face 
avant et l’ouverture dans la paroi du compteur 
sont recouvertes de plomb. Le boîtier est garni de 
plomb à l’avant et d'aluminium dans les zones ne 
pouvant produire de triple coïncidence. L’expé- 
rience a prouvé que le plomb ne donnait pas de 
bon résultats, étant donné, l’impossibilité de le 
dégazer. Il sera remplacé par de l’or. De même le 
dépôt d’or sur la partie avant du tube de champ 
est insuffisant. Le cristal est métallisé par écrase- 
ment sur la surface d’une feuille d'aluminium de 
200 mg/em? pour éliminer les photons lumineux 
produits au voisinage du fil. Le fil du compteur est 
tendu par un ressort et facilement démontable. 


ITT. Compteurs proportionnels. — Pour déter- 
miner la géométrie des compteurs proportionnels la 
nature et la pression du gaz de remplissage, 1l faut 
tenir compte d’un certain nombre de facteurs. 

a) La nature du gaz joue sur la vitesse de dépla- 
cement des électrons ; comme ceux-ci sont formés 
à des distances différentes du fil, il s’ensuit une 
dispersion dans le temps des impulsions ce qui 
impose une limite à la résolution du circuit de 
coïncidence. 

b) La résolution en énergie comme le temps de 
collection varie considérablement avec la présence 
d’impuretés. 

c) La perte d’énergie totale de la particule le 
long de sa trajectoire modifie le spectre d’énergie 
enregistré par le compteur ce qui réagit sur la 
réalisation du système de discrimination des parti- 
cules. 

d) Les réactions produites dans le gaz, l'efficacité 
de détection des photoélectrons arrachés aux parois 
jouent un rôle important dans la production du 
bruit de fond. 


I est évident qu’il y a un optimum dans le choix 
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des paramètres pour un type de particules d’énergie 
donné. Le but de l'installation étant d’étudier en 
une seule fois une large gamme d’énergie faibles 
conduit à admettre que les particules de fortes 
énergies ne laissent dans les compteurs qu’une 
énergie d’un ordre de grandeur voisin de la perte 
d'énergie maximum. Dans ce cas, la dispersion 


statistique est régie par une courbe de Landau dis- | 


symétrique 7-8, et le maximum de la répartition 
correspondante à la perte d’énergie la plus pro- 
bable est donnée par la formule : 


Cme Me x LCmèe ex 
FE (1 — 2) I?(Z) 


NEC In B? +} 


= 


: masse de la particule, 
: énergie de la particule, 
: épaisseur de l’absorbant, 
: constante fonction de l'absorbant, 
: masse de l’électron, 
: atomique de l’absorbant, 
: potentiel moyen d’excitation de l’absorbant 
: facteur de correction. ; 


5 ©& 
(ex 


Pour une distance radiateur scintillateur donnée, 
on peut choisir d'utiliser la plus grande partie de 
celle-ci pour le compteur proportionnel servant à 
mesurer AE, ce qui permet de disposer d’une‘perte 
d’énergie relativement importante en réservant 
l’espace restant à un compteur proportionnel gros- 
sier uniquement destiné à la définition de la direc- 
tion. On peut aussi utiliser également cette dis- 
tance pour construire deux compteurs propor- 
tionnels identiques de bonne performance, séparés 
l’un de l’autre part une certaine distance ; la perte 
d’énergie mesurée par chacun d’eux est alors rela- 
tivement plus faible et la dispersion statistique 
plus grande. Cependant si l’on met en coïncidence 
les impulsions provenant des deux compteurs pro- 
portionnels dans un circuit de Rossi [9], celui-ci 
sort une impulsion égale à la plus petite des impul- 
sions ce qui réduit la dispersion statistique pour le 
parcours que nous nous étions fixés. L'étude des 
diverses combinaisons possibles au moyen de la 
théorie de Landau montre que du point de vue de 
la résolution en énergie les deux solutions donnent 
des résultats équivalents. 

Nous avons choisi la seconde qui en principe per- 
met de diminuer la résolution du circuit de coïnei- 
dence et devrait réduire les chances de coïncidences 
réelles entre les deux compteurs proportionnels qui 
sont relativement éloignés l’un de l’autre. 

Le mélange utilisé est du krypton avec 5 % de 
CO?, La pression moyenne de fonctionnement est de 
20 em Hg. Au début des expériences, nous utilisions 
un four à calcium solidaire du télescope pour la 
purification du gaz, dont la circulation était assurée 
par une petite turbine. Nous l'avons abandonné, les 
résultats obtenus étant équivalents en prenant la 
précaution d'effectuer le remplissage à travers un 
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mélange réfrigérant. Sur une période d’un mois, on 
observe aucune détérioration appréciable des 
spectres obtenus. 
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Fic. 4. — Spectre de distribution d'énergie à la sortie du 


circuit de Rossi pour des protons de 14,1 MeV et des 
deutons de 12,5 MeV (radiateurs de 15 mg/cm?). Les 
courbes sont les courbes théoriques déduites de la théorie 
de Landau. 


La figure 4 montre, comparés avec les courbes 
théoriques, les résultats obtenus avec des neutrons 
de 14,1 MeV, un radiateur de paraffine de 
15 mg/cm? et un radiateur de deutéro paraffine 
de 15 mg/cm?. L’élargissement à la base du pic de 
deutons est dû à la présence de protons provenant 
de la réaction de rupture du deuton. Ces résultats 
sont comparables à ceux obtenus par Kuo et al. [14] 
qui utilisent un seul compteur de définition de 
perte d’énergie. Par contre nous pouvons utiliser 
une résolution en temps plus faible, 

On voit qu’à première vue il n’est pas possible de 
séparer complètement les protons des deutons à 
cette énergie, Pour un réglage parfait 1l restera 
quelques %, de protons dans un spectre de deutons. 


III. Méthode de sélection des particules. — Que 
Pon ait à utiliser la formule de Bethe ou la formule 
de Landau le produit ESE [ôx est à peu près cons- 
tant et proportionnel à M. Une solution au pro- 
blème de la sélection est donc d’effictuer électro- 
niquement le produit. Toutefois aux basses éner- 
gies, il devient nécessaire de rajouter à Æ un 
terme kAE, k n'étant ailleurs pas constant ce qui 
complique les réglages. C’est pourquoi nous avons 
préféré utiliser la méthode consistant à tracer sur 
un oscilloscope la courbe Æ = f(3E 5x). Chaque 
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impulsion en coïncidence produisant un point ayant 
pour coordonnées l’amplitude des impulsions £ et 
SEJdx ; les protons et les deutons se placent sur 
des courbes voisines ; si l’on veut les séparer, il 
suffit de masquer l’une des courbes et de regarder 
l'écran de l’oscilloscope avec un photomultipli- 
cateur dont l’impulsion mise en forme servira à dé- 
bloquer l’analyseur multicanaux. Le montage que 
nous avons adopté permet de former le cache par 
le phénomène lui-même, tient compte automati- 
quement de la non linéarité apparente de l’énergie 
due aux pertes dans le compteur, des résolutions 
expérimentales et permet d’apprécier directement 
la qualité du réglage. 


Formation du cache 


1 
le) 20 40 


Canaux 60 


Fic, 5. — Spectre de cible épaisse. 


Fic. 6. — Montage du photomultiplicateur sur la caméra. 
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On réalise un spectre continu de protons en irra- 
diant à 00 un radiateur épais de paraffine ; l’épais- 
seur étant limitée par des considérations méca- 
niques la moitié de ce radiateur est recouverte 
d’une feuille d'aluminium d’environ 4 MeV d’épais- 
seur, Ce qui provoque la superposition d’un spectre 
décalé vers les basses énergies. On obtient ainsi le 
spectre de la figure 5 qui peut être corrigée en 1rra- 
diant un temps supplémentaire à un angle conve- 
nablement choisi. Au cours de l’irradiation on pho- 
tographie la courbe sur film transparent au moyen 
d’une caméra polaroïde Dumont, la chambre à 
film est ensuite retirée et remplacée par la chambre 
à plaque de cette caméra sur laquelle on.vient fixer 
un montage où le film obtenu précédemment 
occupe rigoureusement la même position que lors 
de la photographie (fig. 6, 7). Un photomultipli- 


Fic. 7. — Porte-cache. Le cache visible sur la photo occupe 
le plan image de la caméra. 


cateur regarde l’écran de loscilloscope à travers le 
cache et le système de lentilles de la caméra. Le fait 
de pouvoir disposer à tout intant de l’obturateur 
permet de contrôler fréquemment les réglages de 
l’oscillographe en cours d’expérience. De faibles 
dérives du zéro peuvent en effet provoquer des 
erreurs importantes. 

Le même procédé utilisé avec une cible épaisse de 
deutéro paraffine permet de tracer simultanément 
une courbe de deutons et une courbe de protons pro- 
venant de la réaction n(d, 2n)p et d'apprécier direc- 
tement le degré de séparation obtenu. 

Deux clichés caractéristiques sont montrés 
(fig. 8-9). I apparaît immédiatement sur ces clichés 
que les coïncidences fortuites sont répartis unifor- 
mément sur la surface et spécialement dans une 
zone voisine de l’origine, le noircissement de cette 
zone les élimine complètement. 

De plus la méthode de production photogra- 
phique du cache présente de gros avantages. On 
peut par reproduction Jouer sur les contrastes, rem- 
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p'acer une courbe transparente par une courbe 
opaque. Comme il n’est pas possible de disposer 
d’une source de deutons dépourvue de protons, la 
sélection se fait en supprimant les protons avec un 
cache inversé. 


Fic. 9. 
La courbe la plus intense correspond aux deutons diffu- 


Cache obtenu en irradiant de là paraffine lourde. 


sés, l’autre aux protons de la réaction n(D, 2n)p (les 
(amplifications sont différentes sur les 2 clichés). 


Enfin, la densité optique varie rapidement d’une 
façon continue sur le pourtour du cache traduisant 
les dispersions d’énergie des spectres d’impulsions 
(fig. 10). (L’obtention d’une telle répartition néces- 
site de grandes précautions dans le réglage de la 
luminosité et la manipulation des clichés, des ré- 


Spectre des impulsions lumineuses 
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F16G. 10. — Spectre des impulsions lumineuses 
vues à travers un cache. 
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flexions des tracés de doigts ou des rayures la per- 
turbent très facilement. 
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FiG. 11. — Élimination des protors. 


En diseriminant l’amplitude des impulsions de 
photomultiplicateur, on peut en quelque sorte 
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« border le contour », par exemple pour éliminer 
complètement les deutons dans les zones de recou- 
vrement en acceptant une certaine perte de pro- 
tons, on établit alors très facilement une courbe de 
correction de pertes en enregistrement simulta- 
nément un spectre continu de protons dans les deux 
sélecteurs multicanaux avec et sans sélection de 
particules. 

Selon une idée de M. Bonnel, nous étudions la 
possibilité de réaliser un cache ou les zones corres- 
pondant aux deutons et aux protons seraient de 
densité optique différentes. Il suffirait de sélec- 
tionner la taille des impulsions lumineuses pour 
pouvoir enregistrer simultanément des spectres sé- 
parés de deutons et de protons. 


IV. Électronique. — Le bloc diagramme de l’ins- 
tallation est représenté figure 12. 

La mise en route et l’arrêt de l’ensemble de l’ins- 
tallation est commandé par l’échelle prédéterminée 
reliée au moniteur. L’intégrateur permet la con- 
duite de l’accélérateur en se référant au flux de 


Fic. 12. — Bloc diagramme de lélectronique. 


PM : Photomultiplicateur. — CP : Compteur proportionnel, — PA : Préamplificateur. 


— À: Amplificateur. 


dence. — C : Coïncidence. — T : Trigger. — RO : Rossi. 


D : Discriminateur. — R: Retard. — TC : Triple coïnci- 


S : Sélecteur. — AL : 


Allongeur. — Mc : Sélecteur multicanal. — E : Échelle. — EP : Échelle prédéter- 


minée. — I : Intégrateur. 


Les alimentations des détecteurs ne sont pas représentées. Les éléments de référence 
de l’alimentation stabilisée commune aux compteur proportionnels sont thermostatés. 


neutrons. Les échelles E, — E, — E, permettent 
d'effectuer les corrections de pertes et de coïnei- 
dences fortuites, de plus leurs indications per- 
mettent de s’apercevoir rapidement d’un déréglage 


du moniteur ou de l’amplificateur des différentes 


chaînes. 


L2s amplificateurs sont de types non saturables 


et les discriminateurs d’un type présentant peu de 
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variations de retard au déclenchement, en fonction 
de l’amplitude des impulsions (Triggers déclanchés 
par la partie négative de l'impulsion différentiée). 

Le circuit du triple coïncidence comprenant trois 
boîtes à retard incorporés est du type à somme qui 
a l’avantage de pouvoir fonctionner en double coïn- 
cidence par simple réglage d’un seuil de discrimi- 
nation. La possibilité d’observer la sortie du cir- 
cuit somme permet de contrôler simultanément 
en cours d'expérience le fonctionnement des trois 
voies de discrimination. La largeur des impulsions 
est réglable de 0,3 à 2 us. Avant l’attaque du circuit 
de Rossi, les impulsions sont légèrement intégrées 
pour diminuer la dispersion d'amplitude produite 
par les fluctuations des impulsions dans le temps. 

Les allongeurs permettent de maintenir les am- 
plitudes correctes sur les plaques de l’oscilloscope 
jusqu’à l’arrivée de l’impulsion d’illumination du 
Trigger T ; cette dernière est régiable en amplitude 
et en durée de 1 à 10 us. Sous réserve de faire les 
corrections de taux de pertes, la limite du taux de 
comptage acceptable n’est pas imposée par l’élec- 
tronique mais par les compteurs proportionnels. 
On fonctionne sans distorsions jusqu’à 16 000 c/s 
en sortie des discriminateurs D, — D, (seuil corres- 
pondant à 4 keV). 

Après chaque nouveau remp'issage ou chaque 
modification de seuil, il est nécessaire de tracer un 
jeu de courbes des résolutions pour différentes 
bandes d’énergie en utilisant une cible épaisse ; il 
faut en effet tenir compte des variations possible du 
temps de collection et des légères variations de 
retard en fonction de l’amplitude qui peuvent être 
introduites par l’électronique. Les courbes repro- 


0,-0, = 0,7 ps 


Nombre 
de coincidences 


F+-250 


AE 
1,2 1,4 


0:-0:. 
F16. 13. — Courbes de résolution en temps. 
Les 3 courbes correspondent à de$ bandes d’énergie d’envi- 
ron 5 MeV chacune (de 0 à 14 MeV). 


duites figure 13, obtenues avec des largeurs d’im- 
pulsions de 1,2 us dans le triple coïncidence indi- 
quent une fluctuation de temps de collection de 
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0,4 us (déduite de la largeur du flanc de la courbe) 


et présentent un plateau montrant la réduction 
possible de la résolution en temps. 


V. Résultats expérimentaux. — Pour chaque 
série de courbe on effectue les réglages suivants : 

— Ajustage des gains très précis des deux 
chaînes dÆ/dx et contrôle de la largeur de la 
courbe de résolution en énergie avec un radiateur 
mince de polyéthylène. 

— Ajustage de gain de la chaîne PM avec pic de 
protons de 14 MeV. 

— Contrôle du réglage du cache avec des 
spectres de cible épaisse y compris l’enregistrement 
du spectre d’impulsions lumineuses pour le réglage 
de l’intensité du spot qui réagit aussi sur les retards 
de l’impulsion de déblocage. 

— Contrôle du moniteur en enregistrant le 
spectre « sur un sélecteur débloqué ou non par 
impulsion issue du sélecteur de bande. Ces tarages 
sont répétés toutes les heures. Pour déterminer 
l’angle 09 vrai on trace la courbe du taux de comp- 
tage en fonction de l’angle avec un radiateur de 
paraffine lourde, la section efficace différentielle 
(n, d) ayant un maximum très prononcé à 00. 


A) COURBES DE RÉFÉRENCES. COMPARAISON DU 
MONITEUR ET DU TÉLESCOPE COMME MONITEURS. 
ABSOLUS. — La figure 14 montre une courbe de 
protons de 14,1 MeV (radiateur d’épaisseur 
13,5 mg/cm?) avec et sans déduction de bruit de 
fond. On remarque l’existence d’une légère traînée 
due aux neutrons diffusés au voisinage de la cible. 


À 
Radiateur Polyélylène 
+ 2000 0° 6,5 mg/cm° 
r-1000 
Ô 20 Li 40 Canaux 
F1G. 14. — KSpectre'de protons à 09, La courbe inférieure 


est la courbe après déduction du bruit de fond. 


L'efficacité du télescope avec un radiateur de 
polyéthylène de 13,5 mg/em? est de 1,044.10—7 
coups par neutrons par stéradian, + 3,5 % à la 


a 
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source (déterminée en utilisant les tables de Bame 
[16]. L'erreur pourra être réduite par une mesure plus 
précise des éléments géométriques et du poids de 
polyéthylène. | 
.. Une mesure simultanée donne un nombre de 
_neutrons par stéradian de 1,54 1011 E 5 % (dont 
1,5 % d'erreur statistique) avec le télescope et de 
1,554 1011 EL 2,5 % avec le moniteur. 
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B) Érunx p£ LA RÉAcTION l°F(n, d)$0. SéLec- 
TION DES DEUTONS. — On utilise un radiateur de 
téflon. 

La figure 14 montre, dans les conditions du ré- 
glage, élimination des protons d’une cible épaisse 
de paraffine légère. 

La figure 15 montre une mesure simultanée avec 
et sans sélecteur des particules dans le cas d’un 


10 20 30 


#40 50 Canaux 


Fic. 15. — Élimination des protons. 


Comparaison des spectres obtenus en irradiant du téflon avec et sans sélection de 
particules en présence d’un bruit de fond de protons important produit par des 


impuretés sur le support de radiateur. 


2  mb/MeV/st 


0 10 20 


Canaux 80 
LEA 


F1G 16. — Spectre des deutons de la réaction 2F{n, d)'°0. 

Les flèches indiquent, compte tenu de correction de pertes 
d'énergie dans le compteur et de la non linéarité de Pélec- 
tronique, les positions} du niveau fondamental et des 
niveaux excités à 1,99 MeV, 2,45, 3,5 et 3,93 MeV. 


bruit de fond de protons important produit par des 
impuretés sur le support de radiateur. L'analyse 


d’une série de mesures sans sélection de particules 
a permis de vérifier que la sélection au moyen de 
l’oscilloscope n’introduisait aucune perte de deu- 
tons dans le domaine d’énergie étudié. 

La figure 16 montre un résultat définitif obtenu 
au bout de 10 heures d'irradiation (bruit de fond 
compris). On voit nettement les pics correspondant 
au niveau fondamental (8,33 MeV) un pic corres- 
pond au premier état d’excitation de 1,99 MeV. — 
Un groupe plus large correspond aux deux niveaux 
suivant d’excitation 3,5 et 3,93 MeV —. La forme 
du premier pic à gauche du fondamental semble 
indiquer qu'il est constitué de la superposition de 
deux pics, l’un correspondant au premier état 
excité (1,99 MeV) l’autre au niveau de 2,45 MeV 
signalé par Holmgren [15]. 

Si l’on admet pour le fondamental, la répartition 
angulaire donnée par Ribes et correspondant à 
L — 0, on trouve (compte tenu de l’angle moyen du 
télescope) une section efficace de 18,8 mb/st + 
8,9 0) alors que la valeur trouvée par Ribes [11] 
serait dans les mêmes conditions de 19,5 mbyst. 


N. B. — Ces niveaux excités ne sont pas mis en 
évidence par Ward. 
og (19 + 29 niveaux excités) 115,6 

o (valeur fondamental) ‘"? 
valeur largement supérieure à la valeur de Ribes et 
o (39 + 49 niveaux excités) ay 

o (niveau fondamental) L 


Le rapport 


le rapport 


33 
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C) SÉLECTION DES PROTONS. RÉACTIONS DE RUP- partie d’une étude en cours sur la réaction n(D 2n)p. 
TURE DU DEUTON. — Les courbes présentées font Sur la figure 17 on a porté les résultats obtenus. 


mb/MeV/st 
8 


ONF A) 
— FRE 


se | 


JOURNAL DE PHYSIQUE _ N°8-9 


l 


mb/ MeV/st / 


je / © + Sans éhrminat! des deutons 


e Avec ” ” ” 


Fic. 17. — Élimination des deutons. 


La figure montre superposés les spectres obtenus avec une cible de 6,5 mg/cm° de para- 
ffine lourde avec et sans sélection des deutons. Le télescope étant orienté à 150, 1 
avec ce réglage il reste 5 % de deutons. 


avec et sans élimination des deutons, le télescope 
orienté à 15,59 (angle non corrigé). L’analyse de la. 
courbe montre un résidu de 5 % de deutons. Le 
petit pic à 13 MeV est produit par la faible propor- 
tion de paraffine légère présenté dans la paraffine : 
lourde. 

Les figures 18, 19 montrent nos résultats aux. | 
positions 09 et 159 comparées aux résultats de 
Seagraves [13] et aux prévisions de Frank et 
Gammel. Ils correspondent à des irradiations de 
4h 30 bruit de fond compris. Là encore, la méthode 
E (MeV) permet de gagner en rapidité et en précision. Ces 
premiers résultats confirment et précisent les dé- 
saccords avec la théorie notamment la proportion 
plus élevée de protons de haute énergie. 

On trouve 


dofdQ{i) — 58,2 + 3,3 dofdA(1605) = 45 + 1,7 mb/st 


Les marges d’erreur pourront être réduites après 
une pesée plus précise du radiateur et une analyse 
de la valeur C/D de la paraffine lourde, | 

Nous tenons à remercier M. Suzor, Ex-Sous- 
Directeur du Laboratoire de Synthèse Atomique et 
M. le P* Magnan, actuel Directeur, pour toutes les 
facilités qu’ils nous ont accordées. 

L'un de nous (P.L.) tient à rappeler tout ce qu'il 


£ (MeV) 


F1G.18 et 19. — Résultats comparés avec ceux de Seagraves doit au regretté P*' Joliot-Curie qui lui avait suggéré 
à 5° et 1605. La courbe en trait plein est la courbe théo- ce travail. 
rique de Franck et Gammel. 

Les points correspondent à nos résultats, les croix à Manuscrit reçu le 9 mars 1961, 
ceux de Seagraves. 


| 
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REVUE DES LIVRES 


Einstein (A.), Signification de la Relativité. Compléments. 
(1 vol. broché de 57 p., 16 X 24 em, traduit de l’anglais 
par M. Solovine et M. A. Tonnelat, Gauthier-Villars, 
Paris, 1960, prix 15 NF.) 

Il s’agit des annexes de la dernière édition (1954) du 
livre « Signification de la Relatioité ». Ce sont donc les 
dernières publications du grand savant. L’une est consacrée 
au problème cosmologique et aux difficultés qu’on ren- 
contre, notamment une durée trop courte pour la vie de 
l'univers. L’autre à la théorie dissymétrique du champ de 
gravitation. Einstein donne un formalisme cohérent pour 
une telle théorie et l’examine en détail. 

Dans des réflexions finales, il écarte l'hypothèse de zones 
de singularité (susceptibles de représenter les corpuscules). 
Il reste encore, dit-il, à préciser des postulats pour les 
conditions aux limites. Enfin, c’est à une insufflsance de 
l’outil mathématique qu’il attribue l’impossibilité actuelle 
de donner l’expiication de la mécanique quantique. 

J. WINTER. 


Thermodynamique des processus irréversibles. (Comptes 
rendus de l’école internationale de physique Enrico-Fermi, 
réunion de Varenna sul lago di Como, juin 1959), 1 vol. 
de 320 p., 17 X 25 cm, Zanichelli, Bologne, 1960, prix 
3500 Lire. (tous articles en anglais). 

Ensemble de leçons traitant les unes des applications, 
les autres des fondements théoriques des relations de 
Onsager : 1. Kzæin (M. J.), Les lois de la thermodynamique. 
2, Munster (A.), Théorie des fluctuations. 3. DE Groor, Ther- 
modynamique du déséquilibre. (Ces trois études couvrent 
plus de la moitié de l’ouvrage et sont consacrées aux fon- 
dements théoriques.) 4. Frescui (R.), Relations d’Onsager 
pour les phénomènes vectoriels. 5. Mazur (P.), Sur les 
bases statistiques de la thermodynamique du déséquilibre. 
6. KzeiN (M. J.), Le principe de la produetion minimale 
d’entropie. 7. KirkworD (J. G.), Mécanique statistique des 
phénomènes de transport. 8 MonrroLL (E. W.), Même 
sujet. 9. Henri (F.), Processus irréversibles dans les solides. 
10. NaAKkAJIMA (S.), Sur la théorie quantique des processus 
irréversibles. 41. Mazur (P.), Sur la théorie statistique 
mécanique du comportement irréversible. 12. ONSAGER et 
Dupino (M.), Propriétés électriques de la glace. (D’après 
le modèle cristallin de Pernal et Fowler.) 13. CINI CAsTA- 
exoziet Ricet, Diffusion dans les liquides simples (annoncé) 

J. WINTER. 


Photographie corpuseulaire (Premier colloque inter- 
national). 1 vol. relié, 451 pages, 250 figures, Publications 
du CN. R.5.,1958,53 NF. 


Le premier colloque international de photographie cor- 
pusculaire s’est tenu à Strasbourg en juillet 1957. Son but 
était une mise au point de nos connaissances sur les méca- 
nismes de l’action photographique des particules chargées 
et l’exposé des résultats obtenus ces dernières années. 

Cet ouvrage rassemble les communications faites au 
colloque sous 8 titres de chapitre : 

— L'image latente dans les émulsions ordinaires. 

— L'image lantente dans les ému'sions corpusculaires. 

— Préparation et proptiétés des émulsions. 

— Quelques caractéristiques des émulsions corpuscu- 
laires soviétiques. 

— Traitement des émulsions corpusculaires. 

— Méthodologie. 

— Applications à d’autres domaines. 

En plus des communications assez spécialisées, on trou- 
vera dans ce livre un certain nombre de renseignements 
généraux et relativement détaillés, notamment sur la sensi- 
bilité photographique, la sensibilisation des émulsions et la 
formation de l’image latente. 

N. ZELLER, 


Demers (Pierre), Photographie corpuseulaire (II) (1 vol. 
relié toile bleue, 464 pages, 360 figures.) Les Presses Uni- 
versitaires de Montréal (Canada) 1959, 

Ce livre est le compte rendu du second colloque inter- 
national de photographie corpusculaire (Montréal, 1958). 
Il s’agit donc plus d’un recueil de communications sur un 
certain nombre de points précis que d’une revue d’ensemble. 
Le plan adopté comporte six parties : 1) Processus de la 
lumière et des rayons À. — 6 communications y traitent de 
l’action des radiations corpusculaires sur les émulsions, 
ainsi que de la théorie du développement de l’image latente 
obtenue. — 2) Préparation et sensibilisation des émulsions 
tonographiques. — Partie importante où sont exposées di- 
verses techniques pour la préparation des émulsions à 
grains ultra-fins requises pour l’ionographie. La sensibili- 
sation chimique et l’hypersensibilisation à la triéthanol- 
amine sont abondamment traitées. — 3) Développement. 
Exposition de divers appareillages adaptés au dévelop- 
pement des couches épaisses. Diverses communications 
envisagent l’aspect théorique du gonflement et de la péné- 
tration du révélateur dans ces couches. — 4) Image latente. 
Propriétés relatives aux trajectoires. — T1 s’agit des relations 
entre les trajectoires et les radiations corpusculaires qui les 
ont engendrées, en particulier discrimination des particules. 
— 5) Appareils spécifiques des mesures sur les traces, —- 
6) Applications représentatives. 
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OBTENTION DE LIQUIDES POLAIRES IONIQUEMENT PURS 


Par GEorGes BRIERE, 


Laboratoire d’Électrostatique et de Physique du Métal, 
Faculté des Sciences de Grenoble. 


Résumé. — Il est possible d'appliquer aux liquides polaires les méthodes, aujourd’hui très pra- 
tiquées, de désionisation des solutions électrolytiques aqueuses par l’emploi de résines synthétiques 
échangeuses d’ions fortement acides et basiques utilisées conjointement [1]. 

L’éthanol, liquide très apparenté à l’eau par son grand pouvoir ionisant a permis de montrer 
comment ce procédé de purification ionique convient aux solvants polaires. 

Il semble bien que nous ayions atteint un état voisin de l’ionisation naturelle ; les mesures 
électrométriques, déjà anciennes, viennent en apporter une confirmation encourageante. 

Après avoir décrit les manipulations qui nous ont conduit à mesurer la conductibilité électrique 
de l’éthanol ioniquement pur, nous tacherons de préciser comment la désionisation peut être envi- 
sagée dans le domaine des liquides organiques en général. 


Abstract. — It is possible to achieve complete deionization of certain pure polar liquids, as is 
usually done with aqueous electrolyte solutions, by column operation on a mixed bed of a strong 
acid anion exchanger and an equivalent quantity of a strongly basic cation exchanger. 

By applying this process to ethanol, a liquid very similar to water as regards its great ionizing 
power, we were able to show how this ionic purification suits equally well to polar solvents. 

It seems that a state close to natural ionization was reached : confirmation in support of this 


is given by alcoholic electrometric measurements. 


After a description of the test set up, an attempt is made to generalize the above method to non 


aqueous media. 


I. Introduction et rappel historique. — L'étude 
des propriétés électriques des liquides est un do- 
maine qui à suscité de nombreux travaux pendant 
ce demi-siècle dernier mais qui, pourtant, n’a pro- 
gressé que très peu et manifeste même une régres- 
sion d'idées dans certaines publications récentes. 

D'une manière générale, la spécificité des résul- 
tats obtenus ne s’est trouvée vérifiée qu’à propos 
de l’eau ; les travaux de Kolhrausch (1895) ont 
permis de déterminer avec précision la constante 
d'équilibre de l’ionisation de l’eau, hypothèse de 
Van’t Hoff vérifiée ensuite par les mesures électro- 
métriques, et de concevoir ainsi la première théorie 
des acides et des bases. 

Dans tous les autres cas, et quelque soit le 
liquide étudié, il est impossible de conclure que les 
valeurs des conductibilités électriques rencontrées 
sont définitives et reflètent le comportement intrin- 
sèque du liquide. 

Distinguons en premier lieu les liquides non 
polaires ou très peu, dont la permittivité n’atteint 
pas 5, des liquides polaires et très polaires, pour 
qui elle atteint et dépasse 20. Pour les premiers, 
hydrocarbures aliphatiques ou benzéniques tels que 
l’hexane et le benzène par exemple, il a été démon- 
tré que par les procédés de purification vulgaires 
(déshydratation chimique, distillation et absorp- 
tion des impuretés par des surfaces actives), il est 
possible d'atteindre des résistivités bien supé- 
rieures à 1014 Qcm et pouvant aller jusqu’à 1020, 


et qu’alors, l’ohmicité disparaissant, la conducti- 
bilité résiduelle est difficilement interprétable (tout 
processus conduisant à une loi de variation expo- 
nentielle du courant avec le champ trouve sa justi- 
fication dans les données expérimentales : émis- … 
sion froide des électrodes, ionisation des molécules 
par le champ électrique intense, etc...). 

Il ne parait pas utile de faire appel à une 1onisa- 
tion spontanée du liquide, indépendante de toute 
action extérieure, étant donnée la petitesse des 
courants minimum limites observés [2]. 

Le rôle contrôlé de l’eau dans ces liquides a été 
précisé récemment et une explication de l’influence 
des pollutions en présence d’eau en a été donnée [3]. 

Il en va bien différemment pour les liquides dié- 
lectriques polaires dont les résistivités obtenues 
dépassent rarement 107 Qcm et varient selon les 
auteurs, et pour un même liquide dans des propor- 
tions qui jouent sur des puissances de dix | 

L'influence des pollutions électrolytiques est 
dans ce cas prépondérante et la diversité des résul- 
tats reflète les différents degrés de purification 
atteints ainsi que l’influence des procédés employés. 
C’est là une manifestation du pouvoir ionisant de 
ces diélectriques, que leur confère une permittivité 
élevée. La plupart des auteurs, après Walden, fai- 
sant l’analogie avec le cas de l’eau admettent une 
ionisation naturelle responsable du courant mesuré 
dans ces liquides prétendus purs. Il a même été 
conclu récemment qu’il est impossible d'attribuer 
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aux semi-conducteurs liquides tels que l’acétone 
et le nitrobenzène, une résistance déterminée [4]. 

Cependant, dès 1914, un expérimentateur soi- 
gneux avait atteint des valeurs particulièrement 
élevées des résistivités de l’acétone et différents 
alcools, valeurs que les travaux postérieurs sem- 
blent ignorer ou mettre en doute, cet auteur fut 
pourtant le seul à considérer ses résultats comme 
erronés par défaut ! Il a fallu un demi-siècle pour 
pouvoir lui donner raison ; en effet, un nouveau 
moyen de désionisation, bien plus puissant qu’au- 
cun des précédents, est apparu et a révolutionné les 
techniques de purification et de séparation, tant 
au laboratoire que dans l’industrie : il s’agit des 
résines synthétiques échangeuses d’ions qui per- 
mettent de produire quotidiennement de l’eau 
dans un état de pureté comparable à celui obtenu 
jadis par Kohlrausch après de très nombreuses 
distillations fractionnées dans un appareillage en 
silice pure ! 

Il est aujourd’hui aisé de retrouver une résis- 
tivité voisine de 25 MQcm (à 18 °C) par élution 
d’une eau ordinaire sur colonne de résines appro- 
priées. L’emploi de ces mêmes résines peut en êtra 
fait afin de désioniser des solutions salines dans 
d’autres liquides, l’échange s’effectue d'autant 
mieux que le liquide est plus polaire [5]. Il doit 
montrer, évidemment, une inertie chimique au 
contact de l’échangeur et celui-ci doit être inso- 
luble dans le liquide d’élution. 

Les vitesses d'échange sont cependant moins 
importantes que dans l’eau et un plus faible écou- 
lement est nécessaire. Il a été montré en 1959 que 
leur emploi convient aussi à la désionisation de 


liquides organiques chimiquement purs [6] et 


qu’alors des concentrations d’ions très faibles peu- 
vent être éliminées. Le choix d’un liquide tel que 
alcool éthylique nous a permis d’en juger l’effi- 
cacité, nos conclusions étant fortement étayées par 
les mesures électrométriques faites sur ce liquide 
pendant les années 1920-1930. Dans le cas parti- 
culier de l’éthanol, il est remarquable de constater, 
là encore, la grande diversité des résistivités mesu- 
rées par les différents auteurs : l’élimination des 

ollutions était dans tous les cas insuffisante. Dès 
1914, Carvallo obtenait après purifications chimi- 
ques et distillations fractionnées successives, une 
valeur de 1,1.108 Qcm 17 0C [7]: Ensuite, réappa- 
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raissent, citées dans les tables de constantes, des 
valeurs bien plus faibles pouvant descendre jus- 
qu’à 10% Qcm (alors que la résistivité de Peau pure 
est de 25 MQcm à 18 0C | 

Danner et Hildebrand en 1922 trouvent 
7,4:108 Qcm à 25 0C, la meilleure valeur rencontrée 
dans la littérature [8]. Cependant les mesures de 
f. e. m. de piles en milieu alcoolique font prévoir 
par le calcul de la constante d’ionisation, des résis- 
tivités 50 à 100 fois plus élevées. 

Nous sommes parvenus à obtenir une concor- 
dance satisfaisante avec les données de l’électro- 
métrie, montrant ainsi que les échangeurs d’ions 
sont un moyen puissant d'élimination d’impuretés, 
même dans des solvants autres que l’eau et à des 
concentrations extrêmement faibles. 


II. Calcul de la conductibilité électrique limite 
des liquides ionisés : eas de l’éthanol. — Il est pos- 
sible de calculer la conductibilité électrique limite 
de l’éthanol sans qu’il soit nécessaire d’effectuer 
des mesures sur le liquide totalement privé d’ions 
étrangers, pour cela deux données expérimentales 
sont nécessaires : La mobilité w des protons et des 
ions ethylates dans l’alcool ; le produit ionique ÆX 
(ou constante de dissociation) de la réaction équi- 
lbrée, d’ionisation du liquide. 

Ces deux valeurs s’obtiennent par extrapolation 
à des concentrations nulle des séries de mcsures 
faites en présence de quantités notables et connues, 
d’ions introduits et devant lesquelles les pollutions 
peuvent être négligées. 

Nous disposons des données nécessaires au calcul 
de o7, à 25°: 

NS CANONS ENS IQ omeA par mole/ko: 
et 
PTE RAT CRE: AN 1e AD 

Si nous posons &« — A/A,, égal à la fraction de 

moles dissociées ou taux d’ionisation : 

Æ = «?.yono—-.ya*t #o?, 
la conductibilité limite de éthanol est alors donnée 
par : 


4 000 4 000 
1/64 = # — = Pl, 
GRAN VAT OS 
Soit 
MOMIE oorones MIO ECM 
ct 29.400 < ecx,om < 62.500 MQ cm. 


TABLEAU 1 


Fo =1 opt .cmr! e MQ x cm : 
Ag 259 QTT.cm AUTOPROTOLYSE CALCULÉE CALCULÉE e (à 25°) 
[2 = le 5 48. EE? 
po a on Zoo = 548 1,04.10—14 Fr S 18,25 MQ em 78,5 
P 0 = 240 rt) 
= 52 s 25.10774 
“es ï, “ He 149 Art LAURE Le 1 280 MQ cm 31,5 
0 SFA 
. É — &R9 16 (O—1 — 1 
C H5 OH A6 CG H3 OH = 83 +5 Q°.cm 4,2 à 15,10— 20 29,4 à 625.10 MO cm 94,3 


par mole/kg. 
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Notons qu’un caleul analogue peut être fait à 
propos de l’alcool méthylique, liquide apparenté 
aussi bien à l’eau qu’à l’alcool éthylique ; à titre 
comparatif nous donnons dans le tableau 1, p. 501 
les valeurs qui nous intéressent ici pour ces trois 
liquides : 


IIT. Obtention d’éthanol ioniquement pur, par 
l’emploi d’échangeurs d'ions. — Nous allons décrire 
maintenant les moyens qui nous ont permis 
d'atteindre un état de pureté ionique, dans l’étha- 
nol, tel qu’une résistivité effectivement mesurable 
a été trouvée, en accord remarquable avec le cal- 
cul précédent. Pour alléger la suite de cet exposé, 
nous insisterons plus sur les manipulations qui 
permettent dès à présent de reproduire un tel état, 
plutôt que sur l’historique de nos expériences. 


1) Déshydratation de l’éthanol. — Nous avons 
utilisé au cours de notre travail, de l’alcool « abso- 
lu », Prolabo R. P., chimiquement pur et dont la 
teneur en eau est variable et oscille entre 500 et 
1 500 10-$ (p. p. m.) à la livraison. 

Si l’on tente une évaluation de la conductibilité 
due à cette eau parasite [9], en appliquant la loi 
d’action de masse à la dissociation de l’eau dans 
l’alcool, on obtient en fonction de la concentration 
d’eau, la résistivité du mélange privé d'ions : 


00 % 1% 100.10—$ 1.10—$ 


107 108 169 1010 Q cm 


D'une manière plus précise on trouve pour 
1 000 p. p. m. d’eau à 259, unc résistivité maximum 
de : 
10+7 x 1 000 


VAG A4 3 Q 
0 x VAU #3 200 MAcm, 


P# 
en admettant que la mobilité de l’ion OH- 
est du même ordre que celle des autres ions mono- 
valents dans l’éthanol (A,, H, O # 100). En fait, 
au cours des premiers essais de désionisation sur 
résine, une telle résistivité a pu être dépassée 
avec de l’alcool contenant environ 5 % d’eau ! Ce 
qui montre bien, d’une part le rô'e du solvant dans 
l’ionisation des corps dissous, d’autre part, que les 
mesures antérieures devaient être faussées par les 
pollutions résiduelles que contenait l’éthanol après 
déshydratation chimique et distillation, et non par 
l'humidité restante. 

L'emploi de résines synthétiques échangeuses 
d'ions a nécessité cependant un dispositif de déshy- 
dratation énergique (et consommant peu d’alcool), 
car l’échangeur lui-même contient environ 50 % 
d’eau de constitution qu’il ne cède que lentement 
au liquide efflucnt. 

De l’alcool « absolu » (à — 1 000 p. p. m. d’eau) 
recyclé constamment sur résine se charge d’une 
teneur maximum limite de quelques % lorsque la 


N°89 


résine est neuvé, puis par suite de la compétition 
entre les affinités réciproques de la résine et de 
l’alcool à la fixation d’eau, cette teneur maximum 


décroit elle-même progressivement, sans devenir 


nulle pour autant, lorsque l’on renouvelle par de 
l'alcool sec, (les mesures sont faites par la méthode 
de Karl Fischer). 

Il nous a été nécessaire de disposer de 30 litres 
environ, d'alcool absolu, pour déshydrater 500 g 
de résine humide, l’alcool d’élution ne contenant 
plus à la fin que 1 500 p. p. m. d’eau, après un recy- 
clage prolongé. Les tamis moléculaires sont un 
moyen indispensable au séchage en courant liquide 
continu. Leur efficacité permet d’atteindre, dans 
l'alcool éthylique des teneurs en eau de l’ordre de 
10 p. p. m., et pour des quantités notables de - 
liquide. Nous avons employé des tamis « Linde », 
type 4 À en batonnets de 1/16 inch, qui possèdent 
un réseau cubique d’alumino-silicate de sodium ou 
de calcium absorbant préférentiellement les molé- 
cules dissymétriques, petites et polaires comme 
celles de l’eau. On peut dessécher théoriquement 
8 1 d’éthanol (contenant 500 p. p. m. d’eau) à 
moins de 50 p. p. m., si l’on dispose d’un kg de 
tamis mis en colonne verticale. Le liquide doit cir- 
culer très lentement dans le sens ascendant avec 
un débit de quelques cm? par minute. 

Les colonnes employées étaient en Pyrex, lon- 
gueur : 130 cm, diamètre intérieur : 30 mm, elles 
contenaient ainsi environ 800 g d’absorbant ; un 
dispositif de régénération constitué par une gaine 
électrique chauffante, porte le tamis à une tempé- 
rature comprise entre 200 et 300 ©C, sous courant 
d’azote sec. 

Le schéma général de l'installation est donné plus 
loin à propos de la désionisation. Notons seulement 
que la circulation est assurée par la pompe utilisée 
au recyclage sur résine, son débit étant réglable et 
compris entre 0 et 12 litres-heure. 

L'alcool, à la sortie du tamis, se trouve chargé 
d’une poussière colloïdale de silicate que l’on éli- 
mine par distillation, afin de ne passaturer par la 
suite les groupes actifs de l’échangeur. 

Dans notre montage nous avons obtenu au cours 
d’un seul cycle, après regénération totale du des- 
séchant, 5 à 8 1 d’alcool contenant 200 p. p. m. 
d’eau en moyenne, en partant d’alcool « absolu » 
à 1 000 p. p. m. environ. Un recyclage de cet 
alcool sur tamis renouvelé, pourrait conduire à une 
teneur en eau inférieure à 50 p. p. m. En fait, la 
limitation est venue de la résine qui cède son eau 
constamment à l’alcool effluent, et de moins en 
moins facilement. 


20 Désionisation en milieu alcoolique. — La 
désionisation dans les solvants organiques polaires 
avait fait l’objet de travaux expérimentaux, uni- 
quement dans le cas de fortes concentrations d’élec-‘ 
trolytes introduits dans un liquide donné, et dont 
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on dosait la quantité disparue après élution sur 
résines à différentes vitesses d'écoulement (10). 
L'efficacité des résines était done, mesurée par 
simple analyse volumétrique. Dans le domaine qui 
nous intéresse, seules les mesures électriques sont 
sensibles à des concentrations d'ions aussi faibles 
que 10-10 jons-grammes/kg, nous avons suivi le 
degré de purification du liquide traité par conduc- 
timétrie, dans la cellule même où la conductivité 
limite devait être mesurée. 

Notons un fait essentiel qui permet de juger 
combien est prépondérant le rôle des pollutions, 
toute modification des raccords entre la sortie de la 
résine et la cellule exige ensuite une longue période 
de lavage (parfois plusieurs journées). Le courant 
liquide entraîne peu à peu les impuretés situées sur 
les surfaces qu’il baigne. A l’arrêt, le liquide se 
pollue par simple solubilisation du verre ; pour les 
mêmes raisons, un essai de mesure avec des élec- 
trodes de platine platiné à montré que le temps de 
lavage devient très long (15 jours ne suffisent pas l), 
les impuretés ioniques occluses dans le platinage 
ne sont libérées que très lentement dans le courant 
liquide et faussent toute mesure. 


a) DÉSIONISATION SUR RÉSINES, À LA TEMPÉRA- 

RATURE ORDINAIRE. — La désionisation, ainsi qu’il 
a été dit précédemment, se fait par élution du liqui- 
de sur un lit de résine constitué d’un mélange en 
parties égales d’échangeur de cations (fortement 
acide) et d’échangeur d’anions (fortement basi- 
que). En effet, pour une meilleure utilisation des 
groupes actifs, le liquide à purifier doit être en 
contact avec les deux variétés acides et basiques 
afin que le PH moyen reste voisin de la neutralité. 
De plus il est certain que les radicaux fonctionnels 
de l'échangeur doivent posséder une fonction 
(acide ou base) forte si l’on veut aboutir à une 
désionisation pratiquement totale. Nous avons 
employé indifféremment les résines américaines 
Amberlite IR. 120 et IR. 400 et les résines alle- 
mandes Merck I et IIT. Les variétés acides sont 
des dérivées polymérisés du polystyrol (ou poly- 
styrène), porteurs de groupes actifs sulfoniques 
substitués, les variétés basiques sont des poly- 
amines substituées en ammonium quaternaire. On 
admet que l’échange s’effectue, dans l’eau, selon 
l'équation de principe : 


R SO, H + CINa 2 RSO; Na + CH, 
(Rh NOH + CINa <= (Rha NCI + NaOH. 


Le radical R étant hautement insoluble dans la 
plupart des solvants ordinaires. Le groupe sulfo- 
nique libère des protons en absorbant les cations 
et le groupe aminé libère des ions hydroxyle en 
absorbant les anions : au total il y a bien échange 
équimoléculaire entre l’impureté saline (supposée 
ici mono-monovalente) et l’eau dans la solution. 
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L’équation d'équilibre citée plus haut laisse 
supposer que léchange s’effectue selon la loi 
d'action de masse mais l’hypothèse d’une absorp- 
tion (pouvant être irréversible à de très faibles 
concentrations) a été proposée par certains au- 
teurs (11). 

1) Appareillage. — Nous allons décrire briève- 
ment le montage effectué. 


Schéma. — L’échange est opéré en circuit fermé, 
construit entièrement en Pyrex : le liquide à désio- 
niser passe sur le lit de résine contenu dans une 
colonne (longueur : 80 cm, diamètre : 35 mm) 
munie à sa partie supérieure d’une réserve de 
600 cm, il traverse ensuite un serpentin situé dans 
un réchauffeur dont la température est thermos- 
tatée, avant de pénétrer dans la cellule de mesure, 
il est enfin recyclé vers la réserve de la colonne. 

Dans nos dernières manipulations, la totalité de 
la colonne de résine était à double paroi et pou- 
vait être portée à la température désirée par une 
circulation auxiliaire de liquide. Nous justifierons 
plus loin ce montage. 

La pompe qui assure la circulation continue de 
l’alcool, à vitesse controlable, doit apporter le 
moins possible de pollutions dans le liquide : un 
modèle, type Vary Pump à membrane, dont le 
corps de pompe et les clapets sont entièrement en 
Téflon, s’est montré très satisfaisant. Des cellules 
de contrôle ont été réparties le long du circuit 
(avant et après la pompe en particulier) pour en 
juger. Le liquide est donc à l’abri de l’atmosphère 
sur la totalité du cycle, il ne prend contact qu'avec 
le pyrex, verre dur, hautement insoluble, le téflon, 
polymère inerte à l’action des solvants, et le métal 
des électrodes. 

Le schéma donne l’allure des opérations succes- 
sives. Tous les raccords entre les différentes parties 
du montage sont assurés par des rodages sphé- 
riques en pyrex, les connections vers la pompe 
étant en tube téflon. 


2) Résultats expérimentaux. — De l’alcool éthy- 
lique ordinaire, contenant plus de 5 % d’eau, après 
circulation sur un lit de résine haut de 50 em, et 
recyclé à une vitesse de 6 litres-heure, passe d’une 
résistivité initiale voisine de 2 MQcm à une résis- 
tivité de 3 500 MQcm, en une semaine, courbe 
1 — Log po —#f(t). La résine est, au préalable, 
mise à macérer dans de l’alcool 969, afin de lui oter 
la majeure partie de son eau. La colonne de résine 
et la cellule de mesure sont alors à température 
ordinaire (25 0C environ). La loi d’Ohm se trouve 
vérifiée pour des champs électriques compris entre 
0 et 2 000 V/cm, même dans une cellule conduc- 
timétrique ayant des tiges de métal pour élec- 
trodes. Des électrodes planes et parallèles donnent 
une linéarité aux courbes 7? — f(V) jusqu’à des 
champs de 4 000 V/cem. 
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SCHÉMA. — 1 : courant d’azote sec. — 2 : Alcoolsec chargé de tamis. — 3 : colonne à distiller. — 4 : Réfrigérant à 
eau. — 5 : Réserve d’alcool déshydraté. — 6 : entrée d’air. — 7 : dessicateur à silicagel. — 8 : prise d’alcool déshy- 
draté. — 9 : réserve d’alcool (960) refroidi venant de la double paroi. — 10 : Réserve d’alcool pur désionisé., — 
11 : chauffe-ballon. — 12 : Colonne de résines échangeuses d’ions {entièrement calorifugée). — 13 : Pompe centri- 


fuge. — 14 : Serpentin de cuivre. — 15 : Vase Dewar (acétone + CO, sol.). — 16 : Cellule de contrôle. — 17 : 
Colonnes de tamis moléculaires. — 18 : Dispositif de régénération à 300 °C. — 19 : Pompe en téflon « Vary Pump ». 


— 20 : Prise d’alcool pur. — 21 : Serpentin réchauffeur (en pyrex). — 22 : 


Thermostat. — 23 : Circulation d’eau. 


— 24 : Cellule conductimétrique. — 25 : Alcool hydraté par l’échangeur. 


Notons d’ores et déjà, que la diminution du 
courant mesuré, ne se fait que très lentement tant 
que le circuit n’est pas entièrement nettoyé et 
qu’un arrêt de la circulation en provoque une 
remontée au bout de quelques minutes au plus, 
parfois quelques dizaines de secondes. Il y a là, 
une preuve de l’importance relative des pollutions 
venant des parois et des surfaces en contact avec 
le liquide pur sortant de la résine. 

— La mise en service d’une colonne plus longue 
(80 cm) et de résine neuve met en évidence la libé- 
ration progressive d’eau dans l’alcool contenant au 
départ 5 % d’eau ; on remplace l’alcool hydraté 


par l’azéotrope obtenu après simple distillation, 
lorsque le courant tend vers une limite ; chaque 
fois la limite du courant est abaïissée (courbe 2 
log p — j(t)). Lors du renouvellement on constate 
une brusque diminution. Pour chaque palier, on 
note la résistivité limite atteinte, elle est fonction 
(décroissante) de la densité de Palcool sortant. 


b) ExXCLUSION D’IONS DANS L'ÉTHANOL ET RÔLE 
DU FROID DANS LA DÉSIONISATION EN MILIEU 
ALCOOLIQUE. — Notre montage a permis de prou- 
ver le rô'e efficace du froid sur l’activité de la 
résine, À cet effet, nous avons fait circuler dans 
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la double paroi de la colonne d'échange, de l’alcool 
refroidi après passage dans un serpentin de cuivre, 
lui-même plongé dans un vase Dewar contenant 
un mélange eutectique acétone-neige carbonique 
(0 — — 80 0C). La température du liquide réfri- 
gérant, après l'échange thermique, permet de fixer 
une valeur maximale de la température du lit de 
résine. La circulation (assurée par une pompe cen- 


. trifuge, Prolabo, à fort débit : 100 à 500 litres/ 


heure) a été effectuée de bas en haut dans la 
colonne, afin que le liquide à désioniser soit en 
contact avec la partie la plus froide du lit d’échan- 
geur avant de pénétrer dans le réchauffeur. 


CE MAQ, cm 


du flacon 
(5% eau) 


0 36 72 
COURBE 1. — p2500 — f{t) sous E — 1 000 V/cm 


# (heures) 


L'ensemble de l’appareillage (cellule de mesure 
exceptée) étant à la température ambiaute 
(8 — 18 0C) la résistivité de l’éthanol à 25 0, 
n’a jamais dépassé 15 000 MQcm, quelque soit le 
degré de déshydratation de l’alcool effluent (pour- 
tant aussi faible que 500 p. p. m.). Par contre, et 
pour une teneur en eau constante, la résistivité 
croit avec l’abaissement de température de l’échan- 
geur, on relève ainsi pour différentes températures 
les valeurs p, 25 °C = f (0 résine), courbe 3. L’allure 
asymptotique des courbes obtenues semble montrer 
que l’éthanol se trouve dans un état ioniquement 
pur et que la résistivité mesurée alors est celle due 
à l’ionisation spontanée du liquide. 

La valeur de la résistivité naturelle limite de 
l’éthanol, serait, en extrapolant à une concentration 
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en eau nulle, voisine de 45 000 MQcm à 25 0C, 
Notons que la plus proche valeur, effectivement 
mesurée est de 43 500 MOcm, l'alcool ne conte- 


12 72 132 156 168 
€ (heures) 
COURBE 2. — gço5oc — f{t) sous E — 1 000 V/cem. 


COURBE 3. 


— p250C — f(0) 


Gén Es 600 P- P. m. 


nant plus alors que 500 p. p. m. d’eau. Cette résis- 
tivité se conserve, indépendamment de la nature 
du métal des électrodes, sous des champs faibles 
allant de 50 à 1000 V/cm. 

Pour des champs plus intenses, le courant croit 
plus vite que la tension mais une étude de cette 
variation nécessiterait l'emploi d’électrodes épaisses 
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à bords arrondis, de manière à éviter les effets de 
pointe sur les angles aigus. Il est remarquable 
cependant de constater que la linéarité des courbes 
I, 25 0C — f (E) disparait pour des champs d’au- 
tant plus faibles que le liquide est plus résistant. 
Les variations de la conductibilité limite, ont été 
mesurées en fonction de la température de l’étha- 
nol lui-même (la résine étant toujours maintenue 
dans les conditions d’efficacité optimum), entre 
20 et 30 °C, la conductibilité est multipliée par un 
facteur 1,6 environ. Dans le cas de l’eau pure, il est 
de 1,84, ce qui correspond bien à une énergie d’ioni- 
sation de l’alcool, très supérieure à celle de l’eau. 


Conclusion. — En résumé, la désionisation en 
milieu alcoolique, sur échangeurs synthétiques, est 
efficace quelque soit le degré d’hydratation de 
l’alcool employé. Il est aisé d’atteindre une résis- 
tivité de quelques milliers de mégohms-cm, avec 
de l’alcool à 969 et une résine neuve n’ayant subi 
aucun traitement préalable. 

Sur une résine longuement déshydratée, et à 
température ordinaire, l’alcool « absolu » du com- 
merce conduit à une résistivité comprise entre 
10 000 et 15 000 MQcm. 

La résistivité naturelle, voisine de 50 000 MQcm, 
ne peut être approchée avec les résines mention- 
nées; qu’à des températures d'échange, inférieures 
à — 50 0€ ;' s’il est possible que d’autres variétés de 
résines plus actives puissent conduire plus aisé- 
ment à de telles valeurs, l’action du froid n’en 
reste pas moins utile, voire indispensable à de tels 
résultats. 
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Il apparait donc que la distillation est un moyen 
de purification ionique insuflisant dans le cas de 
liquides très polaires, peu ionisés naturellement. 
Deux causes viennent en limiter l’action : d’une 
part, la solubilité des parois de l’appareil (rappe- 
lons que les concentrations ioniques détectées sont 
de l’ordre de 10-19 ion.g/kg dans l’alcool ioni- 
quemeñt pur), d'autre part l’entrainement, lors de. 
la distillation de gouttelettes liquides, dispersées 
dans la phase vapeur et contenant des impuretés 
ioniques dissoutes. L’échange sur résines assure 
une désionisation pratiquement totale et de plus, 
l’écoulement du liquide pur permet de lutter contre 
la diffusion des particules solubles issues des parois : 
il s’agit d’une purification dynamique. Notons 
avant de conclure, ce fait remarquable : l’alcool 
éthylique ayant été purifié ainsi qu’il a été dit voit 
sa résistivité décroitre lentement de 50000 à 
500 Mfcm environ, après quelques heures 
d’arrêt total de la circulation, or, c’est précisément 
une telle valeur qui n’a jamais pu être dépassée 
par les procédés habituels. 

D'une manière plus générale, il est permis de 
penser que les échangeurs d’ions doivent conduire 
à des résultats semblables avec la plupart des ” 
liquides organiques polaires. Des résistivités éle- . 
vées ont déjà été obtenues avec l’éther, l’acétone 
et le nitrobenzène. Le rôle de l’eau est, dans tous les 
cas, prépondérant et les modalités de l’échange en 
milieu alcoolique nous seront d’un enseignement 
précieux pour une étude plus étendue des pro- 
priétés des diélectriques liquides. 

Manuscrit reçu le 15 avril 1961. 
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EFFETS A GRANDE DISTANCE 
DANS LA STRUCTURE ÉLECTRONIQUE DES IMPURETÉS DANS LES MÉTAUX 


Par ANDRÉ BLANDIN, 
Service de Physique des solides (BP. 11), Orsay (S.-et-O.). 


Résumé. — L'étude de la structure électronique des impuretés dans les métaux met en évidence 
l'existence d’oscillations de densité électronique se comportant à grande distance en cos 
(2kr r + ©) fr, dans l’approximation des électrons libres L’introduction du réseau périodique 
permet d'étendre ces résultats au cas d'électrons de Bloch. Enfin, on discute l'influence de l’inté- 
raction coulombienne sur les propriétés de l’écran et les oscillations de dersité électronique. 


Abstract. — Studying the electronic structure of impuretes in metals, we point out the existence 
of electronic density oscillations behaving at large distance as cos 2 Apr fr (in the free electron ap- : 
proximation). These results are generalized for Bloch electrorns. We study then the influence of 
Coulomb interactions on the screening and the charge oscillations, 


Introduction. — L'objet de cet article est une 
étude détaillée des propriétés à grande distance de 
l’écran, autour d’une impureté de substitution dans 
un métal. 

En dérivant la règle de Friedel [1] (pour des 
électrons libres) on peut voir qu’il existe, en dehors 
de l’écran électronique concentré autour de l’impu- 
reté, de petites variations de densité électronique à 
grande distance oscillant avec un vecteur d’onde 
2kr et décroissant en r—#. Ces effets à grande dis- 
tance ont été introduits avec succès dans l’expli- 
cation de certaines observations expérimentales : 


a) Le couplage indirect entre spins nucléaires : 
Pour expliquer les largeurs de raies observées dans 
la résonance magnétique nucléaire des métaux purs, 
Ruderman et Kittel [2] ont montré l’existence d’un 
couplage d’échange entre spins nucléaires par l’in- 
termédiaire des électrons de conduction de la 
forme : 

A; A f(2krr)l.b 


où A, et À, sont les constantes de couplage hyper- 
fin des noyaux 1 et 2 et 


æ COS Æ — sin æ 
x 


fe) — 


Ce résultat peut être interprété de la manière 
suivante : le spin nucléaire [, polarise la bande de 
conduction (ici, il s’agit d’une polarisation de spin). 
Le spin nucléaire /, interagit par interaction de 
contact avec la polarisation créée, qui se comporte 
à grande distance comme cos 2kr rr°. 


b) La variation du déplacement de Knight dans les 


. alliages : La variation du déplacement de Knight 
dans les alliages dilués est expliquée par le chan- 


gement de densité électronique Apf(r) à la dis- 
tance r et au niveau de Fermi[3][4]. Pour grand r: 


9 
AeF{r) & sin PEL. 


(La décroissance en r—? est due au fait que, dans 


‘ce cas, on ne s'intéresse pas à la variation totale de 


densité électronique, mais uniquement à la varia- 
tion créée par les électrons au niveau de Fermi.) 


c) Les efjets quadrupolaires dans les alliages 
dilués : Les fortes variations de l’intensité des raies 
de résonance magnétique nucléaire dans les alliages 
dilués ont été interprétées [5], [6] par Pinteraction 
du moment quadrupolaire © des noyaux de la ma- 
trice et des gradients de champ créés par les impu- 
retés à grande distance. L?s gradients de champ 
sont essentiellement dus aux oscillations de charge, 
mais, dans ce cas, l’influence de la matrice (et du 
potentiel périodique) est importante. 

Les trois exemples ci-dessus montrent l'intérêt 
d’une étude systématique des effets à grande dis- 
tance autour d’une impureté dans un métal. Nous 
nous placerons dans le cas d’une impureté de substi- 
tution et ne nous intéresserons qu'aux impuretés 
normales (Zn dans Cu par exemple), en excluant les 
impuretés pour lesquelles les phénomènes 
d'échange prédominent (impuretés de transi- 
tion Mn ou Fe dans Cu par exemple). 

Dans la première section, nous simplifions le pro- 
blème au cas d'électrons libres (c’est-à-dire non sou- 
mis à un potentiel périodique) perturbés par une 
charge ponctuelle Z. L’approximation de Hartree 
des électrons indépendants permet de décrire le 
système par un potentiel à un électron représentant 
l’ensemble des effets de la charge centrale Z et de 
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son écran sur le mouvement d’un électron. Nous en 
déduisons l’existence d’effets à grande distance. La 
deuxième section étend ces résultats au cas d’élec- 
trons de Bloch, soumis à un potentiel périodique. 
Enfin, la troisième partie justifie ces résultats : nous 
montrons comment les caractéristiques de l’écran 
sont liées aux propriétés de la constante diélec- 
trique du milieu. L'influence de l’interaction cou- 
lombienne entre électrons ne modifie pas sensi- 
blement les oscillations à grande distance. 


À. La méthode de Hartree (pour des électrons 
libres). — Cette approximation décrit le mouve- 
ment d’électrons indépendants. L’effet global d’une 
charge ponctuelle Z et de son écran est alors repré- 
senté par un potentiel V(r) ne possédant pas de 
terme coulombien à grande distance. 
connaissons pas en détail les propriétés de V{(r), 
mais Friedel [1] a démontré que le potentiel V(r) 
doit satisfaire une condition de self-consistance : 


2-25 +) nn (1) 


où les m(Ær) sont les déphasages créés par le poten- 
tiel V(r) et décrivant la diffusion au niveau de 
Fermi. Une seconde condition doit être imposée 
à V(r) : que le rayon moyen d’action soit de l’ordre 
de grandeur du paramètre de Thomas-Fermi À tel 
que À? = 4krfrT. 

Ces deux conditions semblent donner un large 
choix de V(r), mais pratiquement seuls les dépha- 
sages 10(Æ) n1(E) sont différents de zéro de manière 
appréciable. En imposant deux conditions, No Eb Ma 
sont déterminés et les propriétés à assez grande dis- 
tance de V{r), ne dépendant que des déphasages, 
sont peu sensibles au choix de V{r). 

Regardons maintenant la variation de densité 
électronique à grande distance de l’origine (les ré- 
sultats présentés ci-dessous ne sont pas nouveaux ; 
nous en présentons une démonstration par ondes 
progressives qui nous sera utile par la suite — 
cf. électrons de Bloch —). 

La fonction d’onde décrivant la diffusion vers 
l'avant d’un électron du vecteur d’onde k et d’éner- 
gie E — k?/2. 

Vé = Vifetker + gkfr)]. 
La variation de densité électronique par unité 


d’énergie est alors (en tenant compte des deux di- 
rections de spins) : 


2k : 
D) EEK) e—iker ce + legtr)l] àQ (2) 


Aprlr) + (27 


et le changement total de densité électronique est 
obtenu en sommant (2) dans la bande jusqu’au 
niveau de Fermi : 


Acalr) k dk, (3) 
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Nous ne. 


N° 8-9 ; 

On voit sur les expressions (2) et (3) que le calcul 
de la densité à l’origine dépend du calcul détaillé 
des fonctions d’onde. Daniel [4] a montré que la | 
densité totale, ou par unité d’énergie, était finie à 
l’origine. 

À grande distance de l’origine, nous pouvons cals | 
culer la forme asymptotique de (2) et (3) : 


gklr) = f& (0) eikr/r 


où fk (0) est l'amplitude de diffusion vers l’avant. 
Alors Aok(r) prend la forme après intégration par 
parties et en conservant le terme principal : 


+ ec 


2r f& (0) 
kr? 


Aoglr) = STE 


LT kr? 


At x) 


— 0 +2 fu (el aa] 


que nous pouvons simplifier en utilisant le théo- 
rème optique [7] : 


fuir an = 7m (0). 


Il vient : 


fe (r) evikr + ce 


Afx'r) = DE ir? 


Cette expression ne dépend que de l’amplitude de 


. diffusion vers l’arrière et peut être écrite en utili-. 


sant la relation donnant l’amplitude de diffusion en 
fonction des déphasages : 


fé (0) — HE (27 + 1) [ein — 11 P, (cos 0) 
Aogfr) = + pa sin (2kr + ©) (5) 


où « et o sont déterminés par : 


&« sin ® = X (—) (27 + 1) sin? m(Æ) (6) 


& cos @ = à © (—} (2 + 1) sin 2m(E). 


Le terme principal du changement total de den- 
sité électronique est donné par intégration par 
parties : 


k $ 
Apr) = — 5 [fk (r) ete + co] 


œ 
es C0 (2kr r + gr) (7 


où ax et r sont pris maintenant au niveau de 
Fermi. 

Nous avons ainsi une description de l’écran qui a 
les propriétés suivantes : 

1. La charge électronique d’écran est égale à Z 
(en satisfaisant la règle de Friedel) ce qui assure 
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x 


l'absence de terme coulombien en 1/r à 
distance. 

2. La densité à l’origine est finie. 

3. Il existe des oscillations de densité éootro 
nique à grande distance, le changement de densité 
totale ayant un comportement oscillant avec le vec- 
teur d'onde 24r et décroissant en 7 #. 

On peut se demander dans quelle mesure l’inter- 
action coulombienne entre électrons modifie ces 


grande 


résultats, en particulier si elle n’a pas pour effet 


d’écranter les oscillations à grande distance. C’est 
ce que nous étudierons dans la section €, mais nous 
allons auparavant étendre ces résultats au cas 
d'électrons de Bloch. 


B. Généralisation au cas d’électrons de Bloch. — 
Dans ce deuxième paragraphe, nous généralisons 
les résultats obtenus dans le premier paragraphe 
pour les électrons libres : nous nous plaçons donc 
dans l’approximation de Hartree, la structure pé- 
riodique introduit quelques nouveaux concepts in- 
téressants. Nous rappelons d’abord sommairement 
les résultats obtenus par Roth [8] dans la diffusion 
d'électrons de Bloch par un potentiel. Nous dis- 
cutons ensuite l’existence et la forme des oscil- 
lations de densité électronique à grande distance et 
montrons qu’elles ne dépendent que de « l’ampli- 
tude de diffusion vers l’arrière » au niveau de Fermi 
(nous préciserons plus loin la signification de la 
« diffusion vers l’arrière »). 


1. DIFFUSION DES ÉLECTRONS DE BLOCH PAR UN 
POTENTIEL V(r). — Les fonctions d'onde électro- 
niques dans un potentiel périodique ont la forme : 


Onkr) = eik.r un!r, k) 


où z,(r, k) est une fonction périodique. Leur énergie 
est Æ,(k). Si l’on introduit un potentiel pertur- 
bateur V(r), on peut montrer que les fonctions 


| d’ondes diffusées correspondant à la diffusion vers 


l'avant prennent la forme : 


Paglr) — onk/r) )—f Go ir, r, E) Vir’) Pa k(r’) dèr 
où Gÿ(r, r’, E) est la fonction de Green (*) pour 
l'énergie E : 
7 1 ” Enk r) enr) 

“ NGC 5 3 . 
RDELE IE Lee LE 

Alors, Roth [8] a montré que la forme asymp- 
totique à grande distance de Wi(r) est : 


Gô (r, 


on k, ki) eiD: 


D = 
PS 


Enk r) rdet %;| /2 Pmki(r). (8) 
AEL 


(*) Nous laissons de côté le problème de l’analyticité des 
fonctions Æn(k) et w,k(r) et supposerons que les expres- 
sions [E — ÆEA(k) + 1e]! etc... ont un sens. Ce problème 
délicat n’a pu, jusqu’à présent, être traité qu’à une dimen- 
sion [9], [10 
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fin(k, ki) est la généralisation dé l’amplitude de 
diffusion : 


1 À Vlr 
finlk, ki) = — 57 f'osulr) PU 


La sommation sur met k; est étendue à toutes les 
bandes (**) et points de l’espace réciproque où 


Wa (r') dr. (9) 


et 
NS 0: 


Vk Emtki) = àr 
est le tenseur de l’inverse des masses effectives : 
| a —= Vk Vk Enlki) 
et ®; un facteur de phase valant 
0 si & est défini positif, 


x si 0; est défini négatif, 


[es 


x [2 si & est non défini. 


On peut avoir une idée physique assez simple de 
ce résultat : Considérons par exemple une structure 
de Bloch à deux dimensions à maille élémentaire 
rectangulaire où la surface E(k) = E est comprise 
dans une seule bande. Dans la direction r indiquée 
sur la figure 1, l’onde sortante comprend deux 
termes correspondant aux deux points de la courbe 
E = E(k) ayant un gradient parallèle et de même 
sens que r. L’onde diffusée dans la direction r a sa 
vitesse de groupe Vk E(k) parallèle à r. 


Fic. 


4. — Diffusion dans la direction r 
dans une structure de Bloch. 


Nous pouvons maintenant discuter les propriétés 
de l’écran électronique accumulé autour de la per- 
turbation, en particulier la charge totale d’écran. 

Dans ce cas, on ne peut pas diagonaliser simple- 
ment la matrice de déphasage n, comme il est 
possible de le faire pour des particules libres diffu- 
sées par un potentiel central. On peut faire un 
calcul direct de la charge produite par la pertur- 
bation, en considérant l’amplitude de diffusion 
f*(k, k”) définie en (8) comme l’élément de matrice 


(**) Par la suite, nous régligerons les indices m ou n, 
sachant que dans toutes les sommes sur k;, nous avons à 
tenir compte de toutes les bandes. 
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de l'opérateur e” sin n entre les états k et k’. En 
fait, la règle de somme de Friedel est une règle de 
trace et elle peut être généralisée à tout système 
(cf. Appendice A). La charge de l’écran pour les 
électrons d’énergie inférieure à Æ est égale à : 


À traces " (10) 
où n est l'opérateur déphasage, diagonal en énergie, 
et la trace est prise sur les fonctions propres d’éner- 
gie Æ. La condition de self-consistance pour «écran- 
ter » une charge ponctuelle Z est que le potentiel 
V(r) dans la structure périodique donne lieu à une 
matrice n telle que 


Z = (2/7) tracerr n. 


2. CALCUL DES VARIATIONS DE DENSITÉ DE 
CHARGE À GRANDE DISTANCE — La variation de 
densité électronique créée par la diffusion de l’élec- 
tron de fonction d’onde de vecteur k est : 


Ak(r) = [FPE — |çKl? = >: à Ofr) œklr) + k, k;) + cc 
+ ct 0 ps glor) fT ke ki JR k) 
kiki ? 
cri 


« 2 2 2 É 

où l’on a noté aifr SPORE 
qui restera constante dans toute la suite puisque 
nous nous intéressons à une valeur déterminée 
de r, définissant ainsi des vecteurs k; et o; et œ& 
ne dépendant que de Æ. La variation de densité 
par unité d’énergie, pour l’énergie Æ a pour valeur 
asymptotique (en tenant compte des deux direc- 
tions de spin) : 


Aeg'r) = (11) 


See Fk'r) Nero 


Dans cette a il y a deux . termes: 


a 
i Sk 
A3 — te { fs 
a) 4 r méur) fau 1° \k, ki) Ven 
qui peut s’écrire 
a *# { : * ; cik,r dSk 
RU D Ve 
r og (r) ufr L 87 ( L i) IVKkE ki] 


où l’opérateur gradient opère toujours à droite et 
alors la fonction u(r, — 1y) n’opère pas sur la 
partie spatiale incluse dans la fonction u. Le terme 
principal pour r grand de l'intégrale 


fr dSk 
IVk Ek)| 


a été calculé par Roth [8]. L’argument de phase 
stationnaire mentionné ci-dessus montre que le 
terme principal de cette intégrale pour grand r 
correspond aux points k; de la surface Æ — cte 
pour lesquels yKk E(k;) est parallèle à r. Le calcul 
exact donne alors : 


cik-rAg = E 2m (e) 
Kj 


Var 5 kg. ce) (12) 
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de telle sorte que : 


Ài CH EET Fa pËlast + (k;, ki) œk;(r r) 
*% 
— af f+ (—k, ki) ?—kilr)]: 
b) Un autre genre de termes est : 


SES dk 
2 où ch) [A5 (6 Ka FFE Ki ge 


17 intégrale est calculée dans l’appendice B et est. 
égale à : 


*# 
“— 


En groupant ces résultats, toute une série de 
termes s’annulent et on obtient pour Apg(r) (en 
utilisant le fait que œk(r) = g#k(r)) : 


Aez(r) fe di a @kalr) okslr) + ki, ki) — ce]. 


_ 27? ir? 


On peut alors déduire la variation totale de den- 
sité en sommant jusqu’au niveau de Fermi Er : 


Er 
setri= fAestr; dE. 


Nous gardons dans cette intégrale le terme prin- 


cipal en intégrant par parties, et en remarquant. 


que ki, k; sont fonctions de l’énergie E pour r 
donné. La primitive de eKi+kÿ-r est : 


etki+kil.r 


7 


arr | IVk al 


et il vient pour le premier terme de Ab décroissant 


enr: 


Tr) FC kiki) 


+ ce 


i 
[Ve Elki | ‘ [Vx Eïk;)| 


où maintenant les k; sont pris au niveau de Fermi. 


3. CONGLUSION. 
oscillations de charge est assez compliquée, mais on 
peut remarquer : 


a) La décroissance en r—? à grande distance. 


b) La seule présence des amplitudes de diffusion 
vers l’arrière f+(ki, — k;). Dans le cas général, on 
peut avoir plusieurs termes s’il existe plus d’un 
vecteur k. Par exemple, dans le cas de la figure 1, 
pour une onde incidente de vecteur K;, il existe 
deux amplitudes de diffusion vers 
1 (ki a 1) je (ki ka). 


c) (14) est la somme de termes oscillants avec les 


vecteurs d'onde k; + k;, modulés par les fonctions : 


périodiques wki(r) wk;(r) prises au niveau de Fermi. 
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— La formule (14) donnant les À 


l’arrière 
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Dans le cas d’une surface de Fermi convexe, (14) 


se simplifie : 


Aer) = — 1 UK Elko)| 


kT° r$ 
Latlek/r)} f+(— ko, ko) + ce} (15) 


Nous voyons donc que les résultats obtenus avec 
les électrons libres se généralisent sans grand chan- 
gement aux électrons de Bloch. Mais nous avons 
précisé en particulier la dépendance périodique 
de Ap(r), superposée au facteur oscillant et décrois- 
sant en 7. Ce fait est important, car il peut modi- 
fier profondément l’analyse en harmoniques sphé- 
riques de la répartition Ab(r) et les effets électro- 
statiques dans une cellule atomique différente de la 
cellule centrale [5], [6]. 

À titre d'exemple, nous allons appliquer les ré- 
sultats (14) et (15) au cas d’une zone de Brillouin 
presque pleine et retrouver pour un gaz de trous 
des résultats analogues au cas des électrons. Soit un 
réseau cubique simple (cf. fig. 2) : on peut alors 
amener l’origine au point A(K), coin de la zone de 
Brillouin, et poser k, —K + kr pour kr petit : 


E = Elk) — (K3[2m*) 


ko r) — cKr) cikre-r. 


ky 


F1G. 2. — Variations de densité pour une sphère de trous. 


On peut remarquer que 


DE fr 


pre) = PAG) = PÉersr) = PElr), 


rétant le vecteur AB de la figure 2. Alors œk(r) est 
réel et l’on a : 


PT Lek)IE UF (= kos Ko) etérr + cc]. (16) 


En comparant (16) avec l'expression (7) obtenue 
pour des électrons libres, il ne faut pas s'étonner de 
la présence de la msse effective m*. Ceci vient de la 
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différence de définition entre f+ (k, k'} et fkr(0) pour 
des électrons libres. En fait : 


m* + ke, kr) = f(x). 


On obtient alors exactement l’expression (7) pour 
les électrons libres, modulée par la fonction 
Cœk(r)I : 


Apfr) M 


pe ar eee RE LS A IEEE 


Il faut remarquer dans (17) que Ao(r) est une 
densité d'électrons. On trouve donc pour une sphère 
d’électrons et de trous la même variation de densité 
de charge à grande distance. Ayant affaire à des 
trous, on s’attendrait à un signe opposé. En fait, ce 
signe est contenu dans la définition de l'amplitude 
de diffusion, où l’on a pris le potentiel V(r) pour 
des électrons alors que si l’on traitait un gaz de 
trous de manière analogue à un gaz d’électrons, 1l 
faudrait considérer le potentiel — V{(r). 


C. Infiuence de l’interaction coulombienne sur 


les propriétés de l’écran. — Dans les sections précé- 


dentes, nous avons décrit de manière phénoméno- 
logique l’influence de l’interaction coulombienne en 
introduisant un potentiel V(r) à un électron soumis 
à certaines conditions : ne pas avoir de terme cou- 
lombien à grande distance, satisfaire la règle de 
Friedel, avoir un rayon d’action déterminé. 

Nous voulons discuter plus en détail, dans cette 
section, la nature et les propriétés de l’écran et 
essayer de tenir compte de manière plus détaillée 
de l’interaction coulombienne. Mais le problème 
d’un gaz d'électrons en interaction est compliqué et 
ne peut être résolu que moyennant certaines simpli- 
fications. L'introduction d’une charge localisée n(r) 
compliquant encore la situation, nous nous borne- 
rons à étudier le cas de petites charges localisées, 
c’est-à-dire à ne considérer, dans les modifications 
apportées au gaz d'électrons, que les termes li- 
néaires en A(r) (nous nous limiterons au cas d’élec- 
trons libres). 

Le premier paragraphe pose le problème et le 
relie aux propriétés de la constante diélectrique. 
Dans le second paragraphe, nous calculons la cons- 
tante diélectrique et la répartition électronique 
dans l’écran dans divers approximations. Enfin le 
troisième paragraphe compare les résultats obtenus 
et discute en particulier la forme et l’amplitude des 
oscillations à grande distance. 


1. RÉPONSE A UNE CHARGE LOCALISÉE. LA CONS- 
TANTE DIÉLECTRIQUE. — Soit dans un gaz d’élec- 
trons une petite charge localisée extérieure n(r). La 
présence de cette charge induit dans le système une 
répartition électronique Ap(r). [Pour un système 
quantique Ap(r) est la valeur moyenne dans l’état 
fondamental perturbé de l’opérateur X d[r—ri]]. 

+ 
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Le potentiel électrostatique V(r) dans le système 
vérifie donc : 


AV{r) = 4nfn(r) — Apfr)]. (18) 


Comme nous nous intéressons aux effets linéaires, 
il est intéressant de prendre la transformée de 
Fourier de (18) : 


— K? V(K) = 4x[n(K) — As K° (19) 
etintroduire, ainsi que le font Nozières et Pines [11] 


la constante diélectrique statique pour le vecteur 
d’onde K, e(K, 0) : 


— e{K, 0) K? V(K) = 4nnK). (20) 
On tire alors de (19) et (20) : 
JE ET es APR}, 
= (K, 0) n{K) 


La charge induite est donc : 


Aetr) (x) LeK, 0) 


= ff nlK) eik-r dk. (21) 
En particulier, pour une charge extérieure 
ZÔ(r) : 


AE elle 0) 1] A Re 


L’étude des propriétés de l’écran revient donc à 
l’étude de la constante diélectrique du gaz per- 
turbé. En particulier : 

— la propriété d’écran, c’est-à-dire la non-cxis- 
tence de terme en 7-1 dans le potentiel électrosta- 
tique global V(r) du système, est liée au compor- 
tement de e{K, 0) pour K petit ; 

_— l’existence d’oscillations à grande distance est 
liée aux singularités de la constante diélectrique. 


2. CALCUL DE LA CONSTANTE DIÉLECTRIQUE STA- 
TIQUE ET DES PROPRIÉTÉS DE L'ÉCRAN DANS DI- 
VERSES APPROXIMATIONS. — La première approxi- 
mation utilisée a été l’approximation de Thomas- 
Fermi linéarisée [1], [12] : elle consiste à écrire que 
la densité électronique perturbée b(E%, r) est : 


(Er, r) = polEr — Vir), r]. 
eL£r, r] étant la densité non perturbée, 
Po — kr [3r?, 
d’où en linéarisant : 
Apr) = — (kr{x?) Vir), 
d’où, avec (19), il vient : 


AT 
K? + x2 


n'K) avec À? = 4kpfr 


et 
ETF(K, 0) — 4 + (4kp/rK?), (23) 


et dans l’espace réel, pour l’écran d’une charge 


ponctuelle Z : 


V{r) = — (Z}r) ex 
APTE {r) = (à? ZJ4nr) e—". 


Sur les résultats (24), nous voyons que l’approxi- 


mation de Thomas-Fermi donne un bon écran à 


grande distance (pas de terme coulombienne dans … 


V(r) à grande distance), mais d’une densité infinie 
à l’origine et pas d’oscillations à grande distance. 


La méthode de Hartree self-consistante considère le. 


mouvement d’électrons indépendants dans le po- 
tentiel V(r). Alors la densité électronique per- 
turbée Ap(r) peut être calculée en fonction de V{r). 
En particulier, dans l’approximation linéaire, on 
peut calculer Ab(K) dans l’approximation de Born 
(cf. Appendice C), et l’on trouve : 


Ap(K) = — (kr/27°?) V(K) 8(K), 
où la fonction g(K) est : 


K? 


ke K + 2kr 
{ — Se DrES 0° 
£g(K) 1 + k (1 2e Log 


K—2kr 


: (25) 


en utihsant (1 ),il vient pour la solution self-consis- 


tante de V{(K) : 


LT k 
= = 7 (A 
KE Gr in) 8K) 


d’où la constante diélectrique : 


VKy= (K}, 


SEK, 0j = 1 + (24r/rK°) glK) (26) 


et la densité électronique d’écran pour une charge 
ponctuelle Z : 


Z g(K) 
AÇSFE r) — è iK-r dŒK. (27) 
FU Er. RTE FKES 


Les résultats (26) et (27) appellent quelques re- 
marques : a) On obtient à partir de (26) la cons- 
tante diélectrique de Thomas-Fermi, en rempla- 
çant g(K) par sa valeur pour K = 0, g(0) = 2. 
L’approximation de Thomas-Fermi, bonne pour les 
faibles valeurs de ÆX, ne peut donc pas bien décrire 
les propriétés d’écran pour les plus grandes valeurs 
de Æ. b) Le résultat (27) avait été obtenu direc- 
tement [13] dans l’approximation de la phase aléa- 
toire, valable à haute densité. Il semble donc que 
Papproximation self-consistante soit identique à 
approximation de la phase aléatoire. On peut mon- 
trer, en effet, que ce résultat est valable non seule- 
ment pour les phénomènes statiques, mais aussi 
pour ceux dépendant du temps. Ce point est discuté 
plus en détail dans l’appendice D. c) L’existence 
d’un point critique logarithmique pour K — 2% 
dans la fonction g(K) : 


8K) = 1 + (K —2kr/2kr) Log |À — 2kr| pour À æ 2kr 


donne lieu, pour Ap(r), à des oscillations de densité 
électronique en cos 2Æyr fr. 
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) On peut essayer de trouver une approximation 

plus générale que l’approximation de Hartree self- 
consistante (c’est-à-dire de la phase aléatoire) 
valable à forte densité. Hubbard [14] par exemple 
a proposé de rendre compte des termes négligés en 
remplaçant la fonction g(K) par 


£g'K) + SE tie 
= Ghrires Ge UN = REG 


Cette expression rend bien compte de l’écran 
pour K petit. Pour les grandes valeurs de K, le 
facteur /(K) vaut 1/2 et représente l'influence des 
termes d'échange. L’interpolation entre ces deux 
extrêmes est empirique. On obtient alors pour la 
densité électronique de l’écran d’une charge Z : 


f TRE À La BG LEURS li MR TR iK 
RL 1 | 
28 


(28) a essentiellement les mêmes propriétés que 
dans le cas de approximation de la phase aléatoire. 

Il est difficile d’aller plus loin dans le calcul de la 
constante diélectrique, mais nous pouvons dégager 
les propriétés importantes. 

a) le comportement à petites valeurs de K, qui 
décrit les propriétés d’écran et l’absence de terme 
en 7-1 dans le potentiel d’une charge localisée ; 

b) le comportement à grandes valeurs de K 
(g(K) «K—?) donnant une densité électronique finie 
à l’origine ; 

c) l'existence d’un point critique logarithmique 
pour À — 2k-F donnant naissance à des oscillations 
de densité électronique à grande distance. Ce fait 
est essentiel pour les problèmes qui nous intéres- 
sent. Il est lié à l’existence d’une surface de Fermi 
et au caractère métallique des premiers états 
excités. 


3. COMPARAISON DES RÉSULTATS OBTENUS. — 
Nous pouvons maintenant comparer les résultats 
obtenus. En tenant compte de l'interaction cou- 
lombienne, l’écran d’une charge ponctuelle Z a les 
mêmes propriétés que dans l’approximation de 
Hartree décrite dans la première section : charge 
totale d’écran égale à Z, densité électronique finie 
à l’origine, oscillations à grande distance en 
cos 2kr rfr°. Un tel écran est représenté sur la 
figure 3 comparé avec le résultat de la méthode de 
Thomas-Fermi linéarisée. Ces propriétés ne sont 
que le reflet dans l’espace réel, du comportement 
de la constante diélectrique e(K, 0). La figure 4 
montre la variation de e(K, 0) avec X dans l’appro- 
ximation de Hartree self-consistante, et dans 
l’approximation de Thomas-Fermi. 

Le point essentiel de cette étude est la mise en 
évidence des oscillations de densité électronique. 
Pour les calculs pratiques, il est intéressant de 
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comparer les amplitudes du terme principal à 
grande distance. En l’écrivant sous la forme : 
Z cos 2krr 
RU algue Meme 
on obtient le tableau I des amplitudes relatives « 
pour diverses valeurs du vecteur d'onde de Fermi, 


FiG. 3. — Variations de densité électronique autour d’une 
charge fixe. 
Méthode de Hartree. 
— — — — — Méthode de Thomas-Fermi. 


F1G. 4. — Constante édiélectrique e(K, 0). 
Approximation de Hartree self-consis- 
tante. 
— — — — Approximation de Thomas-Fermi. 


dans les approximations étudiées (les densités mé- 
talliques usuelles correspondent à kr compris entre 
0,5 et 1 u. a). Dans ce tableau, l’approximation a) 
correspond à la méthode de Hartree avec un poten- 
tiel du type de Thomas-Fermi — Ze*/r (qui 
satisfait à la règle de Friedel) : 

a = 1/(2kr + 2/x) 

b) donne l’amplitude des oscillations dans la 
méthode de Hartree self-consistante (ou phase aléa- 
toire) : 

op = 27° Ar/(27r Xp + 1)?; 
c) donne l’amplitude dans lapproximation de 
Hubbard : 
oc = 27° kr[(27x kp + 3/5)?2. 
34 
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TABLEAU I 


AMPLITUDES RELATIVES DES OSCILLIATIONS 
A GRANDE DISTANCE 


a) approximation de Hartree avec un potentiel-ZeAr/r 
b) approximation de Hartree self-consistante 
c) approximation de Hubbard. 


Ky 0,4 0,6 0,8 1 
en u.' a. 

a 0,70 0,55 0,45 0,38 

b 0,632 701532": Uk 0.98 

0,81 0,62 0,50 . 0,42 


Les valeurs calculées concordent à 15 % près. 
L’approximation de Hartree (approximation à un 
électron, avec un potentiel satisfaisant la règle de 
Friedel) rend donc bien compte de l’existence et de 
l’amplitude des oscillations de densité électronique ; 
elle représente une description satisfaisante des 
propriétés de l’écran formé autour d’une petite 
charge localisée dans un gaz d'électrons libres. 

Nous terminerons ce paragraphe par quelques 
remarques : 

a) Pour les grandes valeurs de r, on peut déve- 
lopper la vraie valeur de Ap(r) : le terme principal 
a la forme cos 2kr r/r%. Dans l’approximation de 
Hartree avec un potentiel régulier, le terme suivant 
est en sin 24r r/r#, alors que les descriptions plus 
correctes de l’effet d’écran (phase aléatoire ou 
Hubbard) donnent un terme en 

Log 2kr r (sin 2krr/r°). 
L'interprétation de ce terme n’est pas claire, mais 
il doit exister. Dans un calcul d’électrons indé- 
pendants, le seul terme ayant une signification phy- 
sique claire est donc le terme en cos 2#r r /r°. 

b) On peut se demander quelle est la validité des 
résultats obtenus dans la section A. — Quand la 
charge centrale n’est pas petite. Il s’introduit alors 
une phase ®#r et les oscillations ont la forme 
cos (2krr + œor)/rè. Il est pratiquement impos- 
sible d’aborder le calcul réel des oscillations et de 
donner un sens à la phase or dans une approxi- 
mation autre que celle de Hartree. Mais la discus- 
sion précédente justifie au moins l’ordre de gran- 
deur de l’amplitude des oscillations. Les résultats 
d’un calcul numérique sont donc certainement 
valables s’ils ne dépendent pas crucialement de la 
phase. 

c) Les résultats obtenus dans cette section pour- 
raient être étendus au cas des électrons de Bloch. 
On obtiendrait des conclusions analogues, après des 
calculs rendus plus complexes par la présence du 

potentiel périodique. Pour avoir une idée sur les 
modifications apportées à l’amplitude des oscil- 
lations, nous pouvons comparer les résultats du 
tableau I à la valeur de l’amplitude relative « cal- 
culée pour un puits carré de rayon atomique rs 
(cette approximation a été utilisée et justifiée dans 
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le cas de métaux réels [15, [4] ; elle revient à écrire | 
d’onde d’un électron sous la forme dr) = q(r) ufr), 
u(r) est une fonction périodique sauf dans la cellule 
de l’impureté, o(r) peut alors être écrite comme la 
fonction d’onde d’une particule libre diffusée par 
un puits carré de potentiel. L’amplitude relative & |! 
est alors, pour un métal monovalent par exemple : | 


= 5 (sin Xp rs — 2kp rs cos 2kr rs) © 0,185. 


Cette valeur est de l’ordre de 2 à 3 fois plus 
faible que celles du tableau I, ce qui donne un 
ordre de grandeur du changement de l’amplitude 
relative des oscillations en présence d’un potentiel 
périodique. 


D. Conclusion. — Nous avons donc montré 
l'existence d’oscillations de densité électronique 
autour d’une impureté. Pour un gaz d’électrons 
libres, ces oscillations sont du type cos (2k4rr + o)/r3. | 


La généralisation au cas des électrons de Bloch met” 


en évidence l’importance du potentiel périodique 
(dans les effets quadrupolaires par exemple [5], [6]). 
Cette étude nous a permis de préciser l’origine des 
oscillations, liée à l’existence d’une singularité du” 
type (À — 2kx) Log |X — 2kr| dans la constante” 
diélectrique. Ces propriétés sont caractéristiques 
d’une substance métallique, c’est-à-dire de l’exis- 
tence d’une surface de Fermi. 

L'ensemble de ces résultats justifie l’interpré- 
tation donnée au couplage indirect par échange [2], ! 
au déplacement de Knight dans les alliages [3], [4]. ! 
et aux effets quadrupolaires dans les alliages [5], | 
[6]. Ils sont à rapprocher de la remarque suivante ! 
de Kohn [16]: la constante diélectrique statique ! 
d’un gaz d’électrons présente une singularité pour | 
K = 2e (dérivée infinie). Si l’on considère que les 
phonons sont créés par le déplacement de charges 
(ions), la courbe de dispersion des phonons doit pré-" 
senter une singularité (rapide changement de fré- 
quence) près de À = 2kr. En effet, on peut écrire 
grossièrement que la fréquence des phonons «œ(K) 
est : 


OK) = PJeK, w), 


Q étant la fréquence de plasma des ions. © est 
petit et présente donc une singularité pour K = 24 
L’argument, généralisé à une surface de Fermi non’ 
sphérique, permet de définir une image de la sur- 
face de Fermi dans la surface de dispersion des 
phonons. 


Appendice A. 
Généralisation de la règle de somme de Friedel. 
F Soit Z1, l’hamiltonien décrivant le mouvement 


d'électrons indépendants, V une perturbation loca- 
lhisée. H = H, + V. 


. N°89 


On définit les opérateurs fonctions de Green non 
perturbées et perturbées (*) : 
1 1 
FÉES AMEN Les Er 
e petit > 0. 
Soient R*(E) et X(E) les opérateurs décrivant la 


diffusion pour une onde sortante [R+(E)]|, ren- 
trante [R(E)] ou stationnaire [ X(E)] : 


F&(E) = V + VGX(E) F*(E) 


{ 
PT. 
(F signifiant partie principale) 


(AA 


K(E) = V + VS 


K(E) 


n(E) et n(E) sont les nombres d’états d’énergie 
inférieure à Æ pour les systèmes non perturbé et 
perturbé : 


E 
no(E) = [ dE’ trace S(E’ — H) 


E 
et n(E) — f dE’ trace S(E’ — H) (A2) 
5 Ô(E— 1H) =. 1(G+F(E)— Gr(E))] 
por 
GE(E”) = GÉ(E”) + GE(E) RE(E”) GE(E), 
d’où 
1 lE 
n(E)—n(f) = f dE 


trace [G+(E’) Rt(£) G+(E’) — G-(E') R—(E) G-(EN. 


Comme on considère une trace, on peut permuter 
_cycliquement les opérateurs : 


“ GE2(E') F*(E') 
= GE’) {1 +... + [GE(E’) 7} +...) 
Mais on a : 
d : 5 
Trace TE [GE(E) Vi? = — p trace GÈ(E) [GE(E) Vpr, 


d’où par intégration : 


pe nolE) = À trace { — DE Gt (E)LG#(E) Vi? + cc} 
= = trace { log [1 — Gt (E)V] — log [1 — G(E)] } 


(A.3) 


| mais on peut écrire : 

Trace [log (À + :B)—log (A—1B)] = 2itrace tg—1!(BA—!) 
et 

D OGHE) —1—5 


1 
E—H 
— K(E) V—1 + inô(E — Hi) V 


V + inô(E — Ho) V 
0 


(*) Nous employons ici le langage de la théorie formelle 
de la diffusion. Voir par exemple Messiah [17]. 
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d’où l’on tire de (A.3) : 


nn (Eee - trace te—1[r8(E— Ho) K(E). (A) 
On peut introduire dans ce formalisme l’opéra- 

teur déphasage n, diagonal en énergie dont les élé- 

ments de matrice entre les états |£,0, > et [£su >> 

(x représente l’ensemble des indices autres que 

l’énergie, définissant l’ensemble complet de vec- 

teurs [Ex >) : 


< Es cultg nlE2@ > TT 
(A.5) 


= — 7 < E,«lK(E:)lEat > S(E, — Fi). 
Alors (A.4) devient : 
n(E) — no(E) = (1/7) 2 Eu|n|Ex > 


que l’on peut noter tracez n c’est-à-dire trace sur 
l’ensemble des vecteurs d’énergie £. 


n(E) — nolE) = (1[x) tracez n. -(A.6) 


S'il y a des états liés d’énergie inférieure à Æ 
dans le système non perturbé ou perturbé au 
nombre de o, et o, il suffit, dans les expressions 
(A2), de restreindre les intégrales au continuum, 
on trouve alors : 


n(E) — no(E) = 6 — 69 + (A.7) 


On en déduit la généralisation de la règle de 
Friedel [13]. La charge d’écran Z(E) localisée par 
les électrons d’énergie inférieure à Æ est (en tenant 
compte des deux directions de spin) : 


Z{E) = 2{(6 — 60) + 


(1/7) traceg n. 


(AS) 


En particulier, si V représente un potentiel self- 
consistant décrivant l’écran d’une charge totale Z, 
on doit avoir : 


(2/7) traceg n. 


VAT 


Densité d’états et déplacement des niveaux. — 
Pour un ensemble de centres perturbateurs loca- 
lisés, on a [18]: 


F#(E) = 2 RY(E) 


Z  FÈ(E) GERE(E).. 


L 
2 hi His üu R£ (E) 


Pour une fable concentration x par unité de 
volume de centres perturbateurs, le changement de 
densité d’état Ab(E) est : 

A6(E) = p(E) — polE) = 


traceg 1. (A.9) 


æ d 
Tr AE 


Soit S(ÆE) le déplacement d’un niveau de Æ, dans 
l’état non perturbé à Æ dans l’état perturbé : 


E — Ey + S(E) et p(E) dE — p0(E0) dE 


d’où 
QE) = SOL Or) — IS(E) eolE) 
PE 5E (E) 


| (Ceci est valable pour petit æ et en dehors des 
valeurs de Æ pour lesquelles do,/dE est infini), 
d’où : 


S{E) polE) = — (x/x) tracer n. (A.10) 


Approximation des bandes rigides. — Dans l’ap- 
proximation de Born, (A.4) devient : 
n(E) — no(E) = — trace S(E — H;) V, 
d’où le déplacement d’énergie : 
S(E) p0(E) = x trace d(E — H;) V. 


En particulier, pour des électrons libres, les élé- 
ments de matrice diagonaux de V ne dépendent pas 


de E et l’on a : 
F)-= 2 [rt dr, 


ce qui est exactement l’approximation des bandes 
rigides. 


(A1) 


Appendice B. 
Le théorème optique pour des électrons de Bloch. 


L’amplitude de diffusion est définie par : 


f+{ ki) = — 5 f over 


PK(r) est la fonction d’onde décrivant la diffu- 
sion vers l’avant d’un électron de Bloch de fonction 
d’onde œk(r). 

Calculons : 


)FHir) &r.. (BA) 


* : ds 
Bikes, k5) = 54 (ke ki) PR ko gente 


1 LA 14 
tes ir) Pr) ki) sfr) 


diSk 
Lonellr ee) a 
Jet KT TX ER] 
(en utilisant le fait que VK(r) = TÉk(r)). 
Effectuant l’intégration sur la surface d’énergie 
constante £ — E(k), 1l vient : 


as _ disk 
ic (KT) pe ET k)| 


— — in[G+ir,r, E) — G-1r,r’, E)] 


L_2 
où Gt(r, r', E) sont les fonctions de Green décri- 
vant la diffusion vers l’avant ou vers l’arrière du 
système perturbé. Mais : 


gate) = Pé(n + ef Vlr) Gti, r', E) cut) dr, 
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d’où il vient : 
B(ki, k) = — à [à r Ve 
É Lot) te og) — opte) Lo) — Fk ll} 
TT (Uk, ko — F4 ko Ki] (B.2) 


En particulier, on tire de (B.2) la relation : 


2 ASE 
JR gere = 


Un calcul tout à fait analogue, en utilisant les 
fonctions de Green non perturbées, donne de même : 


: RU: mft\ (Ke k’). 


dSk 
Ciks kj) D ls ki FH) gr 
= a Æ {+ (ki, ki) — JE (, Kj] = B*(k, kj). (B.4) 
En particulier : 
dSk’ x 
fur RE YEN — à — Imf+(k,W.  (B.5) : 


Les relations (B.2) et (B.4) d’une part, et (B.3) 
et (B.5) d'autre part, sont en fait les mêmes rela- 
tions et se déduisent facilement les unes des autres 
en remarquant que l’amplitude de diffusion vérifie 
la relation aisément démontrable (*) : 


FR RTE FFERKEESRE 


Ces relations ne sont qu’une généralisation du … 


théorème optique dela diffusion des particules 


N° 8-9. 


(B.3): 


libres, reliant la section efficace totale à l’ampli- … : 


tude de duffusion vers l’avant [7] : 


5 — [+ (0)1? 80 = Im 10) 


Appendice C. 


Caleul des variations de densité électronique 
dans l’approximation de Hartree 
pour de faibles perturbations 


Soit un gaz d'électrons au niveau de Fermi 
Ex = kK$/2 soumis à un potentiel V(r) petit. Dans 
l’approximation de Born, la fonction d’onde décri- 
vant la diffusion d’un électron de vecteur d’onde k 
est : 


Pr) = &k r — 2 si Gx(r, r”) Pr’) eik-r” dr’ 


(*) En remarquant que : 


& HV Se SPP 
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où 


eik. (r—r : : 
= Er nf a A Jar dk (n petit > 0) 
est la fonction de Green pour une particule libre. 
On en tire la variation totale de densité de charge 
(en tenant compte des deux directions de spin) : 


Gi\ (Ti tE) 


2 \12 
et) = os f KUPE (I — 1 
d’où 
AREA TE 
<Ue (2x)° pes HS 


gitk'—k) .(r—r") 


8 3 r') ; 
Jfakarr Ne Re er UE 
= mn eik.r ds K 
27) 


{ f _… ef ds k | 
E<kr À? — 4 el 2 k2— (k— K)?— in 


K —k’— ket V(K) est la transformée de Fourier 
de V{(r). L'intégrale sur k est aisément évaluée et 


-. donne : 


she [1 + (1 — el] 


ik pr) L°el 2kr 


que nous noterons r#r g(K). On en tire : 


Appendice D. 


Approximations de Hartree self-consistante 
et de la phase aléatoire. 


Introduisons dans un gaz d’électrons une petite 


charge extérieure n(r, t). Cette charge induit une 


variation de densité électronique Ao(r, t). Le poten- 
tiel électrostatique est donc : 


; Ap'r',t)— n'r’,t) # 
ve, = [= SR LA 


(DA) 


Dans la méthode self-consistante de Hartree, 
Ab(r, t) est calculé en supposant les électrons indé- 
pendants soumis au potentiel V(r, t). Dans l’appro- 
ximation linéaire, 


Le fe-uttK.NPIK, ©) f(K, &) d' Kdw, 


Aer, t) = D 97° 


V(K, «) est la transformée de Fourier spatiale et 
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temporelle de V{r, t) et {(K, w) une fonction cal- 
culée par Hubbard [14] : 


et 


h{z,y) = 0 six > x(x + 2) 
ou six > 2 et y < x(z — 2) 


Le Ts 1 {y he 
= fifa) [sie >» 


et (x — 2) < y < x(x + 2) 

ou,siæ <2 et x(2— x) < y < x(x + 2) 

fe net si æ < 2et 0 < y < x(2—x) 
Rx, — y) = fx, y). 


En prenant la transformée de Fourier de (D.1), 
il vient alors : 


— 47 


KE 4e E f(K, uw 


n'K, w) 


V 


et la constante diélectrique pour la fréquence w 
et le vecteur d’onde K : 


2: 


ECO) nl ie = eo (FD.2) 


(D.2) est exactement le résultat obtenu dans 
l’approximation de la phase aléatoire [11]. La simi- 
litude des deux approximations a été discutée par 
Ehrenreich et Cohen [19]. De la connaissance de la 
constante diélectrique, on peut en effet obtenir 
(cf. Nozières et Pines [11]) : l’énergie de l’état fon- 
damental du système d’électrons non perturbé ; 
les propriétés des plasmas ; la fonction de corréla- 
tion dans le gaz d'électrons ; l’énergie des premiers 
états excités. 

En particulier, l’existence d’ondes de densité 
électronique dans lapproximation de Hartree, qui 
a été beaucoup discutée dernièrément [20], [21], 
est liée aux zéros de la constante diélectrique sta- 
tique e(K, 0). De l’expression ci-dessus pour 
e(K, 0), il ressort qu’il n’y a pas de zéros et qu'il 
ne peut donc pas exister de petites ondes de densité 
électronique. 


Manuscrit reçu le 22 mars 1961. 
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REVUE DES LIVRES 


Problèmes actuels en théorie de la Relativité. (Réunions 
Louis de Broglie de 1958.) 1 vol. de 124 p., 16 X 24 cm, 
Paris, édition de la Revue d’Optique, 1960, prix 12 NF. 
Recueil varié de conférences, expérimentales ou théo- 

riques. Signalons deux très beaux exposés de M. LIGHNERO- 

wicz, sur les ondes électromagnétiques et gravifiques et de 

Mme TonneLaT sur les diverses théories unitaires. Que les 

lecteurs seront heureux de les lire sans avoir à fournir le 

gros effort que nécessitent habituellement les écrits sur la 

Relativité générale ! — (Mesure du temps en Relativité, 

par O. Costa DE BEaAurEGARD. Chronométrie hertzienne, 

par B. Decaux. Décalage vers le rouge des raies spectrales, 
par E. ScHATzMANN. Décalage vers le rouge et vieillis- 
sement de la lumière, par S. MavrinËs. Astronomie de 
position et relativité générale, par J. Lévy. Induction élec- 

+ tromagnétique en relativité générale, par PHAM MAN QuaN. 

Conditions d’isothermie en relativité générale, par Pnau 

Tan Hoanc. Choix et interprétation des système de coor- 

données en relativité générale, par S. KIGHENASSANUY. 

Électricité et relativité générale, par C. LaAncos. Quantifi- 

cation du mouvement interne des masses fluides relativistes, 

par P. Hizzionet J. P. Vicier et les deux conférences déjà 
mentionnées, le tout intéressant et récent.) 
J. WINTER. 


Daniel BLancet Georges AmBrosiNo, Éléments de Physique 
Nucléaire. (1 vol., 16 X 24 cm, 328 pages, Masson et Cie, 
éditeurs, Paris, 1960). 

Comme le fait remarquer M. l’Inspecteur général Baïssas 
dans la préface, en France il n’y a pratiquement pas de 
traité de Physique nucléaire, mises à part les traductions 
des bons ouvrages étrangers. Nous devons donc remercier 
les auteurs d’avoir entrepris d’écrire un tel ouvrage et 
surtout d’y avoir pleinement réussi. 

Il a été écrit pour les étudiants des Facultés des Sciences 
et aussi pour les ingénieurs qui désirent avoir une vue 
d’ensemble précise et utilisable sur la Physique nucléaire. 
Sa lecture demande de connaître déjà la Physique, il ne 
s’agit pas là d’un ouvrage de vulgarisation. Les dévelop- 
pements théoriques sont limités à l’essentiel ; les résultats 
sont toujours indiqués de façon à ce que le lecteur ait la 
possibilité d'effectuer des applications quantitatives. Des 
références bibliographiques nombreuses renvoient aux ou- 
vrages dans lesquels on retrouve le détail des raisonnements 
qui conduisent aux résultats cités et qui auraient alourdi 
beaucoup le texte. 

Les auteurs ont insisté sur l’aspect physique des phéno- 
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mènes, qui est essentiel pour une parfaite compréhension 
des faits, 

} J’ouvrage comporte neuf chapitres : Évolution des 
recherches sur la structure de l’atome, Les réactions nu- 
cléaires et l’énergie atomique, Structure et stabilité de 
l'atome, Les noyaux, Les désintégrations radioactives, Les 
réactions nucléaires, Le neutron, Fission nucléaire, Réac- 
tions en chaîne, Structure et cohésion nucléaires, Les réac- 
tions thermonucléaires. 


Dans un appendice on trouve les tableaux, courbes et : 


données numériques utiles-aux applications. 
Dupuy. 


Duquesne (M.), GRÉGoiRe (R.) et LerorT (M.), Travaux 


pratiques de Physique Nucléaire et de Radiochimie. | 


(1 vol.,16 X 24 cm, 324 pages, Masson et Cie, Paris, 1960) 

Dans cet ouvrage les auteurs présentent les travaux pra 
tiques qu’ils font effectuer aux étudiants, à l’Institut du. 
Radium, pour le Certificat de Physique nucléaire et Radio- 
activité. Il est évidemment très précieux pour ces étudiants, 
mais aussi pour tous ceux qui, dans l’industrie ou dans”les 
laboratoires, ont à s’initier à ces techniques particulières. 

Les bons livres de Physique nucléaire abondent, mais ils. 
ne sufflsent pas pour former un bon expérimentateur, c’est 
pourquoi Frédéric Joliot avait demandé aux auteurs de 
rédiger ce recueil qui concourt fort heureusement à la 
compréhension de cette branche de la Physique. 

D'ailleurs, les rappels théoriques sont nombreux avant 
chaque manipulation ; dans la description du mode opéra- 
toire, les détails techniques sont nombreux et précis ; les 
auteurs ont insisté sur les ordres de grandeur des phéno- 
mènes misen cause, sans la connaissance desquels la réali- 
sation d’une bonne expérience et aléatoire. 

Les travaux pratiques sont répartis en deux groupes ; le 
premier est relatif aux méthodes de détection et de mesure 
des rayonnements nucléaires, le second aux méthodes. 
d’études de la structure nucléaire et des phénomènes chi- 
miques associés à ces réactions. Dans un appendice le 
Dr'Devorer,spécialiste en radioprotection, traite des effets 
pathogènes des rayonnements, des normes de protection, 
de la protection contre les divers rayonnements et des 
unités utilisées en ce domaine. Dans un autre appendice 
sont rassemblés les constantes physiques utiles et des 
énoncés de problèmes. | 

Duruy. 
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LUMINESCENCE DES CRISTAUX D’AgCl PURS ET DOPÉS AUX BASSES TEMPÉRATURES 


Par K. VACEK (*) et J. RINGEISSEN, 
Laboratoire de Spectroscopie et d’Optique du Corps Solide, Institut de Physique, Université de Strasbourg. 


Résumé. — La luminescence des cristaux de AgCI, AgCI + 0,01 mol % CdCl, et AgCl + 0,5 
mol % TICI a été étudiée aux températures de l’azote et de l’hélium liquides. La bande d’émission 


à environ # 800 


est décalée à 77 °K en fonction de l’ion dopé. L'émission dans la région spectrale 


rouge à la température de l’hélium liquide disparaît et la bande d’émission principale est décalée 


vers les grandes longueurs d’onde. 


Abstract. — The luminescence of AgCI, AgCI + 0.01 mol % CdCl, and AgCI + 0.5 mol %, TICI 


crystals was measured at liquid nitrogen and helium temperatures. 
was observed at 77 0K in function of the doped ion. 


band at about 4 800 À 


À displacement of the emission 
The emission 


in the red region of the luminescence spectrum disappears at 4 0K and the fundamental emission 
band is displaced towards lower energies when the temperature decreases. 


Introduction. — La luminescence des cristaux 


. d’AgCI purs a été mesurée par plusieurs auteurs [1], 


[2], [3], [4], essentiellement à la température de 
l’azote liquide. Les résultats des expériences ont 
montré que le spectre d'émission se compose de 


. deux bandes, l’une à environ 4 700-4 800 À 


et l’autre à environ 6 400 À. La luminescence des 
cristaux d’AgCI avec 0,5 mol % TICI a été mesurée 


. par Berger [5]. A la température de — 180 °C, 


il a trouvé deux bandes d'émission à environ 4 700 
et 6100 À. L'interprétation de la nature de la 
bande d'émission près de 4 800 À a été donnée par 
Matyas [6]. Il a supposé que cette bande d’émission 
apparaît par la recombinaison des électrons de la 


bande de conduction avec les centres formés par un 


 Jhélium liquides. 


trou positif lié, soit à une lacune cationique (a), 
soit à une lacune cationique sur un cran de la dis- 
location (b). Dans ce travail, nous avons étudié les 
spectres de la luminescence des cristaux purs et 
dopés de AgCI aux températures de l’azote et de 


> 


LL 

Conditions expérimentales. — Les cristaux uti- 
lisés ont été préparés par recristallisation lente de 
la substance fondue entre deux lamelles de verre, 
par la méthode décrite ailleurs [7]. Trois types de 
cristaux ont été étudiés : AgCI pur, AgCI + 0,01 
mol % CdCl, et AgCl + 0,5 mol % TICI. Les 
spectres d'émission ont été observés à l’aide d’un 
spectrographe « Huet AIT ». Les plaques utilisées 
étaient des plaques « Super Panchro Lumière ». La 
densité optique des spectres obtenus a été mesurée 


- avec un densitomètre enregistreur. La position du 


maximum des bandes était déterminée à moins 
de 5 À près. Les échantillons étaient plongés direc- 


* tement dans l’azote ou l’hélium liquides. Les 


(*) Institut de Physique, Université Charles IV, Prague. 


cryostats étaient ceux utilisés au Laboratoire de 
Spectroscopie et d’Optique du Corps Solide [8]. La 
source pour l’excitation des cristaux était une 
lampe à vapeur de mercure « Philora » qui ne laisse 
passer que la lumière ultraviolette à partir de 
4 078 À. Les conditions expérimentales pour les 
différents échantillons étaient les mêmes. 


Résultats obtenus. 

a) TE = TT KES 

Les courbes de la luminescence des cristaux étu- 
diés à la température de l’azote liquide sont repro- 
duites sur la figure 1. On voit que la position du 
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F1G. 1. — Le spectre de la luminescence (1) à la tempé- 
rature de l'azote liquide. Courbe 1 AC 2 


AgCl + 0,01 mol % CAC, : 3 : AgCI + 0,5 mol %TICL 


maximum dépend de la nature de l’ion dopé. En 
même temps, dans le cas des cristaux de 
AgCI + CdCl,, l'intensité de l’émission est beau- 
coup plus faible que dans le cas des cristaux purs. 
Les cristaux de AgCI + TICI ont, dans la région 
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spectrale rouge, une forte bande d’émission qui 
. pourrait correspondre [5] à la superposition de deux 
bandes, l’une correspondant à la présence des ions 
T1+ (6 388 À) et l’autre à la présence des ions excé- 
dentaires Ag+ (— 6 200 À). 

b} T==24°9K. 

Aux très basses températures, la position des 
maxima des bandes principales de tous les cristaux 
reste la même, mais est décalée vers les grandes 
longueurs d’onde par rapport à leur position spec- 
trale à la température de l’azote liquide (fig. 2). 


22 2 KV 


GG 


unilés arbitraires 


N 


8000 à enlA] 


F1G. 2. — Le spectre de la luminescence (I) à la tempé- 
rature de l’hélium liquide. Courbe 1 : AgCl; 2: 
AgCI + 0,01 mol. % CdCl, ;3 : AgCI + 0,5 mol. % TICL. 


Pour le cristal pur de AgCI, ce décalage correspond 
à environ 0,08 eV. L’émission rouge des cristaux 
d’AgCI et AgCI + Ti+ disparaît à cette température. 
Dans le tableau I, on a indiqué la position spectrale 
des bandes. 
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Discussion. — Les résultats obtenus montrent, 
d’une part que la bande principale d'émission est 
décalée vers les grandes longueurs d’onde quand la 
température s’abaisse à 4 0K, d’autre part que la. 
forte bande d'émission des cristaux d’AgCI + TI* 
dans la région spectrale rouge disparaît. Le premier 
effet peut s'expliquer en supposant l’existence de 
deux types de centres de luminescence (a, b).. 
Dans [9], on a indiqué que chacun de ces deux | 
centres a une série de niveaux excités et que la 
bande d'émission observée est finalement formée 
par la superposition des différentes transitions ra- 
diatives des électrons de la bande de conduction à 
ces niveaux. On pourrait se demander si ce dépla- 
cement signifierait qu'à 4 ©K, les transitions de 
luminescence seraient plus favorables pour les 
niveaux d’énergie plus élevés et que, par contre, à 
77 OK, elles seraient plus favorables pour les 
niveaux d'énergie plus bas (*) (par rapport à la 
bande de conduction). La disparition de la bande » 


d'émission rouge des cristaux de AgCI + Ti+ ob-. | 


servée à 4 °K est cependant difficile à expliquer. 
L'influence de l’impureté sur la bande d’émission 
principale est observée seulement à la température … 


de l’azote liquide quand les effets de polarisation du 


réseau par l’impureté peuvent influencer les tran- 
sitions radiatives. Le fort abaissement de l’intensité u 
d'émission à 77 CK dans le cas des cristaux de 
AgCI + Cd++ est probablement dû à la décrois- 
sance de la concentration des lacunes cationiques 
libres sous l’influence des ions Cd*+ qui peuvent 
constituer des complexes du type Agf Cdt+ 
(lacune cationique + ion Cd++). 

Ce travail a été effectué au Laboratoire de Spec- 
troscopie et d’Optique du Corps Solide à l’Univer- 
sité de Strasbourg. Il a été rendu possible grâce à 


TABLEAU TI une bourse du Gouvernement Français auquel l’un 
Crisrac TEMPÉRATURE Are panne 2e paxre de nous (K. V.) exprime sa reconnaissance. Nous 
(°K) (À) (À) tenons à remercier M. le PTS. Nikitine, Directeur 
La ae, ue nu du Laboratoire de Spectroscopie et d’Optique du 
AgCl 77 a 809 4 850 6 200 Corps Solide pour ses précieux conseils ainsi que … 
AgCl + Cdt+ 77 L 880 4 954 ee Mme Kiradijeffova pour la préparation des cristaux, 
AgCl + TI+ 97 PE RSR 6 388 Manuscrit reçu le 21 avril 1961. 
AgCI À 4 954 5 026 re NU, À | 
ne Fes ette tentative d'explication est basée sur la théorie 
AgCI + se 2954 5 026 de Matyas [6]. D’autres explications ont également été 
AgCI + TI+ 4 954 5 026 — envisagées. Elles sont encore à l’étude. 
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RAYONNEMENT DE FREINAGE DANS LES EXPÉRIENCES DE DIFFUSION MUON-PROTON 


Par JEAN-MaARIE ABILLON et Pauz KESSLER, 


Laboratoire de Physique Atomique et Moléculaire du Collège de France. 


Résumé, — En vue d'expériences de diffusion muon-proton dans une région où le transfert 
d’impulsion est de quelques centaines de MeV, nous avons calculé à l’aide de la méthode des pro- 
cessus quasi réels la contribution du rayonnement de freinage à la section efficace de diffusion. Le 
résultat obtenu doit être additionné avec la correction de Schwinger qui inclut à la fois les effets 
radiatifs virtuels et l’émission de photons réels de très basse énergie. On obtient ainsi la correction 
radiative totale, qui est fonction de l’énergie initiale du muon, de l’angle de diffusion, et de la lar- 
geur de résolution pour l’énergie finale. Une application numérique est donnée pour une énergie 
initiale de 2 GeV et différentes valeurs de l’angle de diffusion et de la largeur de résolution. 


Abstract. — In relation with muon-proton scattering experiments in a region where the mo- 
mentum transfer is several hundred MeV, we have calculated by the method of quasi-real pro- 
cesses the contribution of bremsstrahlung to the scattering cross-section. Our result has to be 
combined with the Schwinger correction, which includes both virtual radiative effects and emission 
of real photons with verylowenergy. One thus gets the total radiative correction as a function of 
initial energy of the muon, scattering angle, and resolution width for the final energy. A nume- 
rical application is given for an initial energy of 2 GeV and different values of the scattering angle 
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and the resolution width. 


Introduction. — Des expériences de diffusion 
de mésons y sont actuellement en cours ou en pré- 
paration dans divers laboratoires. Le transfert 


* d’impulsion (du muon au noyau-cible) se situe 
. dans ces expériences entre 10? et 10% MeV. Leur 
but est de déceler dans ce domaine une anomalie 


éventuelle dans les sections efficaces de diffusion 
élastique : anomalie qui serait due soit à une struc- 
ture du méson y, soit à une interaction non élec- 
tromagnétique entre le muon et le noyau. 

En raison de la faible intensité des faisceaux 
de mésons y susceptibles d’être obtenus (à partir 
de faisceaux de mésons rx se désintégrant en vol), 
il est pratiquement impossible dans ces expériences 
d’avoir une bonne résolution sur l’énergie finale 
au voisinage du pic de diffusion élastique. Les 
sections efficaces mesurées incluront donc forcé- 
ment une part d'effets inélastiques. Ces effets, 
dans le cas particulier de la diffusion muon-proton 
(que nous considérons ici) se réduisent au nombre 
de deux : rayonnement de freinage, et production 
de mésons x. Il est donc nécessaire de calculer 
ces effets (compte tenu des valeurs prévisibles 
pour la largeur de résolution sur l’énergie finale), 
pour pouvoir les soustraire des sections efficaces 
expérimentales et obtenir ainsi les valeurs correctes 
des sections efficaces de’diffusion’élastique.t Dans 
le présent travail, c’est au calcul du premier des 
deux effets mentionnés que nous nous attaquons. 

Notons que nos calculs seraient également (et 
même a fortiori, puisque nous utiliserons une appro- 
ximation relativiste) valables pour la diffusion 
électron-proton dans le même domaine d’énergie, 
seule la masse de la particule diffusée variant d’un 


problème à l’autre. Mais du point de vue expéri- 
mental, nos résultats sont beaucoup moins utiles 
dans le cas de l’électron, pour lequel il est facile 
d’obtenir une bonne résolution, et en fait d'isoler 
le pic élastique. Signalons néanmoins que nous 
avons fait précédemment un calcul.semblable pour 
la diffusion électron-noyau dans un domaine 
d’énergie où la structure nucléaire ne se manifeste 
pas encore [1]. 


Caleul du rayonnement de freinage par la mé- 
thode des processus quasi réels. —— Nous utilise- 
rons pour ce calcul la méthode d’approximation 
dite des processus quasi réels [1]. Cette méthode 
consiste à associer à un diagramme de Feynman 
élémentaire tel que celui représenté sur la figure 1 


un spectre donné par : 
CAR NE 
2} k E® Es)’ 
soit encore par l'expression équivalente : 


i dE’ E?+E" 
D PRPE TERRE 


€ 


3 
P(E, à) dk =? [in NE 
T m 


; œ E 
P(E, E) dE’ © (in — Emax 


Cette formule est valable à condition que le 

fermion entrant et sortant soit relativiste 
(E, ES m). 

Elle permet dans de nombreux cas en électro- 
dynamique de ramener un processus complexe à 
un processus plus simple en « amputant » le dia- 
gramme de Feynman correspondant. 

Revenons à notre problème. Soit m et M les 
masses respectives du méson u et du proton, 
E et E’ respectivement l’énergie initiale et l’énergie 
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finale du méson w, 0 l’angle de diffusion, et A la 
largeur de résolution pour £’. Appelons encore 
E* l'énergie finale correspondant au pic élastique, 
g le transfert d’impulsion pour la diffusion élastique, 
et enfin Æ6(— E* — A) la limite inférieure sur 
l’énergie finale (coupure expérimentale du spectre). 
Nous prendrons un système d’unités où c — 1, 


FAGAT 


L’allure générale du spectre différentiel qui nous 
intéresse (pour un angle 0 déterminé) est indiquée 
par la figure 2. Ce spectre contient le pic élastique 


plus la queue due au rayonnement de ‘freinage ; 
il se confond avec le spectre mesuré, une fois que 
l’on a déduit de celui-ci l’effet de production de 
mésons. Notons qu’il est impossible de séparer 
strictement le rayonnement de freinage du pic 
élastique, étant donné que dans le premier phé- 
nomène le photon émis peut avoir une énergie 
aussi voisine de zéro que l’on veut. Nous pouvons 
néanmoins opérer une telle séparation en intro- 
duisant un paramètre artificiel e qui définit une 
limite inférieure pour l’énergie du photon émis. 

La section efficace de rayonnement de freinage, 
intégrée sur le spectre d’énergie, dépendra de 
quatre paramètres : d’une part, les limites À et 
e ; d'autre part, deux paramètres que l’on est libre 
de choisir parmi les trois quantités Æ, q et 0 (liées 
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“1R à 


e € À, et d’autre part nous postulons pour les 
paramètres q, 0 et À les conditions expérimentales 
suivantes : 


m<& q < M ; 8 < 1 radian ; À < Ef2. 


On vérifie aisémént qu’avee ces conditions nous 
sommes dans le domaine de validité de la méthode 
des processus quasi réels qui nous permettra ici 
de relier de façon directe, par la formule donnée 


N°89. 


entre elles par une relation). Nous choisissons. 


plus haut, le rayonnement de freimage à la diffu-. 


sion élastique. 
En effet, nous avons à considérer les deux dia- 
grammes (a) et (b) de la figure 3 ; chacun de ces 


a) 


ice 


diagrammes contient un diagramme élémentaire 
du type de la figure 1, « branché » sur un processus 
de diffusion élastique. Les conditions posées en- 


traînent le fait que les angles d et y indiqués sur 


les deux diagrammes de la figure 3 sont essentiel- 
lement petits par rapport à l’angle total de diffu- 
sion 0. Il en découle deux simplifications impor- 
tantes : 

19 Pour chacun des diagrammes (a) et (b), 
l’angle de la diffusion élastique (qui intervient au 
vertex marqué d’une croix) peut être sensiblement 
confondu avec 0. 

20 Nous pouvons négliger l'effet d’interférence 
entre les deux diagrammes, l’état final étant es- 
sentiellement différent dans les deux cas : dans 
le processus (a) le photon est émis à peu près dans 
la direction de la particule sortante, dans le pro- 
cessus (b) à peu près dans la direction incidente. 

Soit ca, 6m et on les sections efficaces respectives 
de la diffusion élastique et du rayonnement de 
freinage correspondant au processus (a) ou (b). 
Pour la simplicité de l’écriture, nous définirons ces 
sections efficaces comme fonctions de l’énergie seu- 
lement, considérant les autres paramètres (6, À, €) 
comme fixes ; d'autre part, nous écrirons o pour 


la section efficace différentielle do /d Q. Enfin, nous 


À 


.N°8-9 


_ écrirons directement les facteurs de branchement : 


Se — oM(E) JoalE) ; 5 — ch(E) Joa(E). 


On obtient immédiatement : 


E*—£e 
Vus fi ; 
Eo 


fE— 
3b — Je © PE, E*) dE* ca(£'*) Jou(E), 
0 


P(E*, E’) dE’, 


où P(E*, E') et P(E, E'*) sont données par la for- 
mule des processus quasi réels. Les énergies in- 
termédiaires Æ* et E'* se déduisent de Æ et E' 


par les relations bien connues de la cinématique. 


de la diffusion élastique : 
2E QE Ye 0\=1 
Le — ed . er / PE psEpe Er IE 
E = Ef1 +7 sin?) LE, = E'(1 mr Si ;) ; 
d’où également : 
\ ; ” 2 HÉSRC Nr 
EeSE, (1 or sin’) 3 
Le calcul de 5° donne de façon immédiate : 


2 E* 1 se 
D — ER. En 
ë T (in m Paz :) 


MS 
2 


(in E*— FE (E* — . 
£ 


LEY2 


Dans l’argument du premier facteur logarith” 


mique, il est raisonnable de prendre, d’après les 
hypothèses faites plus haut : Pmax & 0 ; d'autre 
part, on peut prendre E* = E et noter que £0 = q. 


D'où : 
_2& qe 1 A A A? 
(nes) (ne + rpm) 
Le calcul de 3? est considérablement plus diffi- 
cile. Il faut considérer en détail la section efficace 
de la diffusion élastique d’un fermion (électron ou 


muon) relativiste par un proton (2) : 
m.?7 008 *.0 F2 
” 4E%sin4 6/2 


Fr + 


— 1 
aulE) (1 Le 2 FAT 12) 


a 


Ke 


ml + AP:)*te" 0/2 + K° Fi] | 


Dans cette formule, X représente le moment ma- 
gnétique anormal du proton (ÆX — 1,79) ; F, et 
F, sont respectivement les facteurs de forme élec- 


. trique et magnétique du proton, qui dépendent du 


quadrivecteur q (énergie-impulsion transférée dans 
dla diffusion). Il est difficile d’avoir, à la lumière 
des données expérimentales (sur la diffusion élec- 
tron-proton), une expression analytique précise 
de F, et F,. Si nous utilisons pour la distribution 
de charge du proton le modèle dit « exponentiel », 
relativement simple et assez bien confirmé dans 
un large domaine d’énergie, ce modèle nous donne : 


PS (A + j, at 4) 


. RAYONNEMENT DE FREINAGE 
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où l’on prend pour la constante a (« rayon » 
électromagnétique du proton) la valeur de 0,8 
fermis 4 (GeV)-1 ; d’autre part, q? est donné 
par : 


à : 2E . rs 
qg? = 4E? sin? 0/2 (1 + Sr sint0/2) . 


Comme nous avons à considérer dans notre pro- 
blème un rapport de deux sections efficaces rela- 
tives à des énergies différentes, nous ne garderons 
dans l’expression de ca que les facteurs qui varient 
fortement avec l’énergie. Compte tenu des condi- 
tions initiales, cela nous donne : 


celle ETF, 


D’autre part, dans l’expression de F, nous pren- 
drons (compte tenu toujours de nos hypothèses) : 
g? = 4ËE? sin? 0/2, ce qui nous donne: 


Fr (1 +2 E?)-4 


où 


Avec ces simplifications, on a : 


E? (1 [2 E?}* 


E—e 
RE PUESRES) LE n + ppp 


soit encore : 


nee Le 
à AU DRE RL RENTE 


E? SE E’*? (1 + L? Fe 


E +2 U +E PEL 
L'intégration, assez fastidieuse, donne fi- 
nalement : 
PP) :) In E — E'* 2 
DEEE LEE RENE ver * 
D — (in — [ + JUEUES | 


où J(LE, LES*) sera donné par l'expression ci-dessous. 
Pour simplifier l’écriture, nous posons : 
LUE = Le) 
et d’autre part : 
u=1+a,u, = 1 + rè. 
Nous avons : 


DA 
Uo t? 


1. CEE" 
Re A TION DU VE 
f(x) Z In ne + Z (u 4) In id 


de - [(Suÿ + Gu? + 8u + 16) 


— 2 (350 + 30u? + 24u + 16)] x arc tg Re 
Ne 


4 ua 
eg À SU — 20 0) 23 + 216 — xo lu + 6)] e 


524 


x? + 8 


2 
x 


u 
u 


+ 3[8 — x, (ou? + Gu + 8)] 
äl u? 
RE 2[2(2 + 9x$ + 624) — xo x(8 + 252$ + 156)] 73 
t uG 

+ DIEU + 228) — 20 e(5u + 6) (828 + 2)] 


2 
+ 3[8 — 167 — 3x0 (Su? + Gu + 8)]}u 


la8 
— (92 + x? + 311at + 33148 + 105a). 


La somme de &* et 3° nous donne la contribution 
du rayonnement de freinage dans les limites dé- 
finies. Toutefois, si nous voulons déterminer la 
correction radiative totale, et en même temps 
éliminer le paramètre artificiel e, nous devons 
ajouter un troisième terme 9° qui est la correction 
radiative calculée par Schwinger [3], et qui cor- 
respond : 

19 Aux effets virtuels, c’est-à-dire à l’émission 
et à la réabsorption d’un photon virtuel ; 20 A 
l'émission d’un photon réel de très basse énergie 
(Éee). 

Ce terme est donné dans le cas relativiste par : 


Laprq À E 47 
Ko [(in m 5) (in € 42 72 
où®D(6) est une fonction compliquée, qui se réduit 
cependant pour 8 << 1 radian à : 


Te œ(0)], 


Il ne reste plus qu’à additionner les trois termes 
pour obtenir la correction radiative totale : 


20 


D = 


[in 1 :) [Aur, LE, x) 


E? ASS si 


In 5 se 
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9 . 
APPLICATION NUMÉRIQUE. — Une expérience qui 


est actuellement en préparation au C. E. R. N. 
prévoit la diffusion de mésons 1 d’une énergie 
de 2 GeV sur des protons ; l’angle de diffusion 
pourra varier entre 109 et 300 (ce qui correspond 
à un transfert d’impulsion compris entre 340 et 
1 020 MeV), et la résolution sur l’énergie finale 
sera certainement de l’ordre de quelques centaines 


de MeV. 
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Nous donnons dans le tableau ci-après les résul- | 
tats que nous obtenons pour à (en pour-cent) en. 
fonction de différentes valeurs de 0 et A. 


0 JA (GeV) 0,80 0,55 0,30 
300 14,60 5,17 0,50 
200 8,25 AT ECS O0 
100 0,40 US er 0,8E 


Il convient de noter que pour 0 = 300 (9 = M), 
nous sommes déjà à la limite de nos conditions 
initiales ; d’autre part, il semble, d’après des ex- 
périences récentes [4] sur la diffusion électron- 


proton, que le modèle exponentiel ne soit plus va- | 


lable dans ce domaine, la décroissance du facteur 
de forme électrique (en fonction de g?) apparais- 
sant comme moins rapide que ne le prédit ce mo- 
dèle. Les approximations que nous avons faites 
d’une part, le choix du modèle d’autre part, nous 
ont sans doute conduit à une surestimation des 
valeurs de à pour 300 ; les résultats donnés pour 


cet angle peuvent donc être considérés comme des | 


bornes supérieures. 

D'autre part, le fait d’avoir utilisé une approxi- 
mation relativiste introduit pour l’ensemble des. 
résultats une erreur relative qui est au maximum - 
de l’ordre de (In g/m) 1. Cette erreur est évidem- 
ment particulièrement sensible pour 6 — 100 


(q = 3m), mais il apparaît que pour cette valeur n | 


la correction radiative est de toute façon pratique- 
ment négligeable, quelle que soit la résolution. 

Malgré leur caractère approximatif, les chiffres : 
donnés ci-dessus (chiffres qui permettent d’ailleurs 
une interpolation facile pour des valeurs intermé- 
diaires de 0 et A) peuvent être d’une grande utilité 
pour les expérimentateurs. Ils doivent en effet 
leur permettre de déterminer la résolution à ob- 
tenir pour un angle de diffusion donné (autrement 
dit, de décider des moyens de détection à mettre. 
en œuvre pour le faisceau de muons sortant sous : 
cet angle) en fonction de la borne qu’ils entendent 
assigner à la correction radiative. 

Nous désirons remercier MM. L. Goldzahl et 
J. Banaïgs, avec qui nous avons eu de fructueuses 
discussions. Nous tenons d’autre part à signaler 
que ce travail a bénéficié de l’aide du Commissa- ” 
riat à l'Énergie Atomique, ainsi que du Centre 
National de la Recherche Scientifique. 


Manuscrit reçu le 9 mai 1961. 
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… Résumé. — La désintégration de 1%Au a été étudiée à l’aide d’un spectromètre 8 à focalisation 
intermédiaire et des techniques de comptage par scintillations. On a déterminé le rapport : 
r = B (E2; 2°—2)/B(E2 ; 2’ —> 0) par la mesure des intensités des raies de conversion. On 
compare les énergies des trois premiers états excités et les valeurs de > pour les différents isotopes 
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pair-pair du platine. 


Abstract. — The decay of Au has been investigated with an intermediate image 6 ra 


spectrometer and with scintillation techniques. 
has been determined by measurement of conversion electron intensities. 


Ni 
The ratio r = B(E2 ; 2’ — 2)/B(E2 ; 2’ — () 
Comparison of the 


energies for the first three excited levels and the 7 value is given for different even-even Pt 


nuclei. 


Les noyaux pair-pair dans le domaine transitoire 
entre les régions déformées et sphériques ont été 
récemment beaucoup étudiés du point de vue théo- 
rique et expérimental. La théorie de Davydov et 
Filippov [1] partant d’une conception d’un rotor 


asymétrique explique certaines des propriétés de 


ces noyaux. Récemment Jean [2] a montré que les 
noyaux pair-pair peuvent être aussi traités par une 
théorie plus conventionnelle, en introduisant un 
modèle qui effectue une transition graduelle entre 
les oscillations quadrupolaires des noyaux sphé- 
riques et les vibrations et rotations des noyaux dé- 
formés. Les résultats obtenus ainsi sont qualita- 
tivement semblables à ceux de Davydov et Filippov 
au moins en ce qui concerne la variation du 
rapportr = B(E2 ; 2° — 2)/B(E2 ; 2° — 0) avec 
E2"+)/E0-—+). 

Du point de vue théorique il est intéressant de 
pouvoir comparer les états excités de différents 
isotopes du même élément. Il existe de nombreux 
travaux sur les noyaux de 1?Pt, 1%4Pt et 1%6Pt 
[3, 4, 5, 14, 15]. Par contre les transitions de 1°Pt 
jusqu’à ce jour ont été peu étudiées. Ce noyau a été 
identifié mdépendamment par G. Albouy et col. [6] 


. et Anderson et Ringh [7] comme le dernier membre 


de la désintégration 1°°Hg —> 190Au > 190P4, Les 


mesures de G. Albouy ont été effectuées par des 


méthodes de spectrométrie y sur des sources sépa- 
rées isotopiquement. Andersson en utilisant la 
même méthode pour la préparation des sources 
s’est servi d’un spectromètre 8 comme appareil de 
mesure. [l signale deux transitions à 295 keV et 


_ 301 keV dont, pour des raisons d’intensités, la pre- 


mière va au niveau fondamental. 
Nous avons repris l’étude de Pt afin de pou- 
voir déterminer les propriétés de ses premiers ni- 


(*) De l’Institut de Recherches Nucléaires, Varsovie. 


veaux excités et les comparer avec les autres iso- 
topes du platine. Nos sources ont été obtenues par 
irradiation d’or massif (1,7 g) avec des protons de 
94 MeV du synchrocyclotron d'Orsay et ensuite par 
séparation isotopique [8] de °Hg. La séparation 
faite sur cette quantité d’or durait habituellement 
5 minutes, temps nécessaire pour fondre le métal 
dans la source d'ions du séparateur. Ce temps de 
5 minutes a été déterminé par collection de “#Hg 
contenu dans la même masse d’or et par observation 
de l’activité collectée directement sur un compteur 
à scintillations monté au deuxième étage du sépa- 
rateur [8, 9]. Pour la spectroscopie y la collection 
des ions de Hp se faisait sur une feuille d’alumi- 
nium. Les sources pour la spectroscopie 8 ont été 
déposées directement sur un support en acétate de 
cellulose avec un dépôt d'aluminium de 1,3 mm de 
diamètre. Pendant les premières 90 minutes après 
la séparation nous nous sommes servis de ces 
sources pour observer les transitions de #°Au [10]. 
Après ce temps l’activité de He (7,y, — 20,4 min 
devenait négligeable par rapport à celle de 1°Au 
(Pr = 59,9 ma). 

Nous avons mesuré le spectre y à l’aide de cris- 
taux de Nal(T1) de 1 1/2” X 1” et de 3" x 3". Les 
valeurs des énergies des principaux pics observés 
sont en très bon accord avec les résultats de 
G. Albouy et coll. [6]. La seule différence est l’exis- 
tence d’un faible pic d’annihilation de 511 keV ainsi 
qu’un pic photoélectrique de 440 keV. En compa- 
rant l'intensité du pic d’annihilation avec la 
raie XX nous trouvons pour le rapport f*/cap- 
ture X une valeur voisine de 2 %, un peu plus 
grande que la limite supérieure de 1 % assignée 
par [6]. Nous avons étudié aussi le spectre y en 
coïncidence avec la raie la plus importante de 
300 keV. Le spectre résultant montre un fort pic 
apparaissant à l’énergie — 300 keV, ce qui indique 
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qu’il y a deux transitions y d’environ la même 
énergie en coïncidence. D’autre part une partie au 
moins de la raie de 600 keV, déjà mise en évidence 
par [6], est en coïncidence avec un ou les deux y de 
300 keV. 

Pour pouvoir déterminer le rapport r nous avons 
étudié le spectre des électrons de conversion à l’aide 
d’un spectromètre à focalisation intermédiaire type 
Siegbahn-Slatis. La détection des électrons se fai- 
sait soit à l’aide d’un compteur G. M. soit d’un 
cristal d’anthracène. La résolution du spectro- 
mètre était — 1,2 % dans le meilleur cas. 

Nous avons étudié plus particulièrement les deux 
domaines d’énergie des électrons de conversion 
correspondants aux transitions de 300 keV et de 
600 keV. Huit raies d'électrons de conversion ont 
été mises en évidence et leurs énergies, intensités 
relatives et attributions proposées sont données 
dans le tableau 1. 


TABLEAU I 


RAIES DE CONVERSION OBSERVÉES 
DANS LA DÉSINTÉGRATION DE 1%Au 


ÉNERGIE keV INTENSITÉ ATTRIBUTION 
217,0 100 295 — K 
229,6 BASE 301 —K 
240,6 8 + 1 319 — À 
282,7 54 +10 295 — LL 
289,8 DANCE 0 SUR LL 

+ 295 — M 
519,6 7 0,7 598 — K 
544,7 15 + 0,5 623 — X 
580,0 faible 898-— L 


Nous estimons que les valeurs relatives en Bb 
sont exactes à 3.10% près et les valeurs absolues 
à 102. On voit qu'aux erreurs relatives près la 
transition de 598 keV peut être considérée comme 
le cross-over d’une cascade d’un y de 301 keV et 
d’un y de 295 keV. En accord avec Andersson nous 
trouvons que la transition de 295 keV est plus 
intense et va au niveau fondamental. Nous ne pou- 
vons pas décomposer suffisamment bien les deux 
raies de conversion L et M de 295 keV et 301 keV 
pour donner la limite de pourcentage de A1 dans 
cette dernière transition. Étant donné que dans 
les autres noyaux de platine cette limite n’est pas 
supérieure à  %, nous allons considérer par la 
suite la transition 2°+ — 2+ comme Æ2 pure. 
Ainsi, d’après les intensités relatives des raies de 
conversion K, on peut calculer le rapport des pro- 
babilités réduites de cascade cross-over : 


r — 85 + 30. 


L'erreur sur cette valeur de r ne tient compte que 


JOURNAL DE PHYSIQUE 
des erreurs commises sur les intensités des raies de 
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conversion. La présence de mélange de M1 dans 
la transition de 301 keV pourrait diminuer ce rap- | 
port. Ainsi, en supposant qu’il y a 10 % de Hi on 
trouve r — 55 au lieu de 85. £ : 4 
Sur la figure À nous comparons la‘variation de ce 
rapport en fonction du nombre de neutrons W pour 


200 


100 


0 | 125241 114 16 118 

B(E2 2° — 2) 
B(E2°2’— 0) 

en fonction de N pour les noyaux pair-pair du platine. 


F1ic. 1. — Variation de 


les différents isotopes du platine. Les valeurs expé- 
rimentales pour ?Pt et I6P£ sont celles données 
par Van Patter [3]. Pour Æ4P£ nous avons calculé 
ce rapport d’après les données précises d’intensités 
des raies de conversion À mesurées par Backstrôm 


et coll. [4]. Nous avons utilisé les valeurs théoriques | 


de «x calculées par Rose. 

Dans 1%Pt, outre les trois transitions 2’ + — 2 +, 
2'+ + 0+et 2+ — 0+ le spectre des électrons 
de conversion révèle deux autres raies que nous 
avons interprétées comme X de 319 keV et X de 
623 keV. Il est intéressant de remarquer que, si 
cette interprétation est correcte on peut composer, 
aux erreurs relatives près, les deux sommes sui- 
vantes : 

E,(295 keV) + Ey(623 keV) = 918 + 3 keV 
E,(598 keV) + Ey(319 keV) = 917 + 3 keV. 


Aïnsi on pourrait avoir un niveau de 918 keV qui 
serait le troisième niveau excité de 1°0Pt. Les résul- 
tats des coïncidences y-y qui montrent l’existence 
d’une raie de 600 keV en coïncidence avec les y de 
300 keV confirme aussi cette conclusion. Ce niveau 
correspondrait au niveau mis en évidence dans 
H4Pt [4] et aussi dans 192Pt [5, 12, 131. 

La théorie de Davydov et Filipov prévoit l’exis- 
tence d’un niveau 3+ pour lequel l’énergie peut 
être déterminée par la relation : 


EG+) =EQ2+) + L(2'+. 


Les récentes mesures de corrélations angu- 


"3 
ne, 


N° 8-9 


- laires [11] dans 1*?Pt ont montré que le spin du 


niveau de 921 keV dans 1°2Pt est 3. 

D’après [4] le niveau à 923 keV dans 1*Pt à pro- 
bablement le spin 44, mais la possibilité de 3 + ne 
peut pas être exclue. 

Sur la figure 2 nous représentons la variation de 
l'énergie des premiers niveaux excités des isotopes 


Eev 
1000 


_—— 
_— 


— 


x 
Points théor (1) 
pour 3 + 


112 : 114 116 118 


… Fic. 2. — Comparaison de l’énergie des premiers niveaux 


excités dans les noyaux pair-pair de platine. . 


} du platine avec le nombre de neutrons N. Sur cette 


figure nous avons aussi indiqué le niveau de 918 keV 
de Pt et la courbe théorique [1] pour le ni- 


. veau 3+."(Dans la littérature on annonce l’exis- 


tence d’un autre niveau dans ?2Pt à l’énergie 
785 keV. Nous ne l’avons pas marqué sur la figure 2 
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d'autant plus qu’il semble y avoir un désaccord 
entre [5] et [12, 13] en ce qui le concerne. Nous ne 
pouvons rien dire sur un tel niveau dans 1°0P4t,.) 

D’après la figure 2 on voit que l’énergie des pre- 
miers niveaux des noyaux pair-pair de Pt change 
d’une façon monotone avec le nombre de neu- 
trons NV. Le point tombant nettement en dehors des 
lignes droites est le niveau 2°+ dans Æ6Pt, le noyau 
ayant comme autre particularité le rapport r très 
grand. 

Il est intéressant de remarquer que dans la théo- 
rie développée par Davydov et Filippov ainsi que 
dans celle de Jean, le rapport E(2'+)/E(2+) ne 
peut pas être plus petit que 2. Dans les noyaux de 
platine il est plus petit pour le nombre de neutrons 
N —=114, 116 et 118 mais plus grand pour Y —112. 
Ce fait ne semble pas avoir une influence sur les 
propriétés de ces noyaux. 

Nous remercions Mme N. Marty pour l’aide et 
l’intérêt constant qu’elle a portés à ce travail et 
pour les nombreuses discussions sur les résultats 
de nos mesures. 

Ce travail n’a pu être mené à bien que grâce aux 
nombreuses séparations isotopiques effectuées par 
l’équipe du séparateur d’isotopes de M. Bernas et 
plus particulièrement M. Sarrouy et M. Camplan et 
nous les en remercions vivement. 

Nous remercions l’équipe de synchrocyclotron 
pour les irradiations qu’elle a effectuées. 

L'un de nous (J. J.) remercie les autorités fran- 
çaises pour une bourse Louis de Broglie qui lui a 
été attribuée au cours de ce travail. 


- Manuscrit reçu le 29 avril 1961. 
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LETTRE À LA RÉDACTION D 


LE FERROMAGNÉTISME FAIBLE 
A L'APPROXIMATION 
DU CHAMP MOLÉCULAIRE 


Par A. HERPIN, 


Service de Physique du Solide et de Résonance Magnétique, 
Centre d'Études Nucléaires de Saclay. 


Considérons une substance magnétique dont les ions 
peuvent être décomposés en deux sous-réseaux À et B 
également peuplés, la symétrie étant telle qu’on ne 
puisse passer de À à B par inversion autour d’un centre. 
D’après Dzialoshinski et Moriya on peut écrire l’éner- 
gie magnétique, à l’approximation du champ molé- 
culaire, sous la forme 

W = 2M Mg + D M /\ Ms —HMi + M5) 
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M, et M3 étant les aimantations des deux sous- 
réseaux, À une constante positive et D un vecteur que 
l’on suppose dirigé suivant Oz. On en déduit les champs 
moléculaires agissant sur les atomes A et B : 

Ha = H—5Me ED N Ms 

Hg =, H— 1M1 — D AM. 


a) À haute température, nous aurons donc : 


C 

Ms = GI — M8 + D AM] 
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Mg = IH —2M; —D A Mal. 


Si H est parallèle à Oz, la susceptibilité est, comme 
pour un antiferromagnétique normal 
br US 
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où — CA est la température de Curie — 0, 
Si au contraire H est perpendiculaire à Oz, par 


exemple parallèle à Ox, les aimantations M, et Ms 


ont pour composantes : 
Ci—T 
C2? + D?) — T?° 
CD 
C?(? + D!) — T? 
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MX = — ME = C. 
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ce qui montre l’existence de composantes transverse: 
MX et Mÿ. On en déduit pour l’inverse de la susceps 
tibilité ct 
L'ENTOERTE ( Te — <) # 


x C 0—T 
où T, — CAR + D?. Le module de D étant très 
petit devant à, on voit sur la formule précédente que su 


a un comportement analogue à celui d’un antiferro= 
magnétique, sauf au voisinage de 7, température pour, 
laquelle Pinverse de la susceptibilité tend rapidement 
vers zéro, comme le montre la figure 1. # 


b) Au-dessous du point de transition T., les deux. 
sous-réseaux ont une aimantation M, — M, en l’ab4 
sence de champ, les vecteurs M, et M; étant dans le 
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plan perpendiculaire à D et faisant entre eux un angle 


(fig. 2) 4 
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D'où une aimantation spontanée très faible, 
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Pour un champ parallèle à D, l’aimantation a une. 
composante qui est parallèle au champ et lui est pros. 
portionnelle : 
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Pour un champ perpendiculaire à D, l’aimantation 
se présente sous la forme 
D Pa el 
= 2x M\ 2 CES É Ë 
C'est-à-dire que l’aimantation présente à la fois les 


propriétés de PE d’un ferromagnétique { aimantation 

DM, 
spontanée 9x 
tion linéaire avec le champ, la susceptibilité différen- 
tielle étant indépendante de la température). 


) et d’un antiferromagnétique (varia= 
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